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This  dissertation  presents  studies  of  charmed  baryons  the  A^,  Ec,  and  £*.  The 
two  A+!  baryons  are  the  L  =  1  excited  states  of  the  A+,  and  they  decay  to  a 
A+  and  two  pions;  the  Ec  triplet  is  the  spin  1/2  7  =  1  state  corresponding  to 
the  A+;  and  E*  is  the  spin  3/2  excitation  of  the  Ec.  Both  Ec  and  E*  decay  to  a 
A+  and  a  pion.  We  performed  comprehensive  measurements  on  the  masses,  widths, 
fragmentation  functions,  and  the  production  ratios  of  these  particles  from  10  GeV/c2 
e+e~  annihilations. 

A  new  decay  channel  of  the  A+x,  A+7r°7r°,  is  first  observed  and  its  decay  ratio  rela- 
tive to  A+7T+7T-  is  reported.  A  detailed  study  of  A^  decays  through  the  intermediate 
Ec  is  performed. 

The  Ec  intrinsic  widths  are  first  measured  with  consistency  with  theoretical 
predictions.  The  relative  production  ratios  among  the  three  Ec's  are  measured  to 
be  consistent  with  1.  The  first  evidence  of  the  existence  a  new  particle,  the  E*+,  is 
reported. 
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CHAPTER  1 
THEORY  OVERVIEW 

1.1    Quark  Model 


Tremendous  progress  has  been  made  in  high  energy  physics  since  the  beginning 
of  this  century.  We  now  have  a  simple  picture  about  the  ultimate  constituents  of 
matter:  All  matter  is  composed  of  quarks,  which  carry  fractional  electric  charges, 
and  leptons,  such  as  the  electron  and  neutrino.  A  quark  does  not  exist  individually 
and  strong  interactions  hold  quarks  together  in  hadrons,  which  are  built  from  two 
types  of  quark  combination:  a  baryon  consists  of  three  quarks  and  a  meson  consists 
of  a  quark  and  an  antiquark. 

Quarks  come  in  six  flavors  (called  d,u,s,c,b,  and  t),  as  do  leptons  (three  types 
of  charged  and  three  of  neutral).  Each  quark  has  spin  1/2  and  baryon  number  1/3. 
Table  1.1  shows  the  additive  quantum  numbers  (other  than  baryon  number)  of  the 
three  generations  of  quarks.  By  convention  each  quark  is  assigned  positive  parity,  so 
each  antiquark  has  negative  parity.   The  masses  of  quarks  vary  largely  for  different 


Table  1.1:  Additive  quantum  numbers  of  quarks. 
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flavors,*  ranging  from  a  few  MeV/c2  for  u  and  d  quarks  to  about  170  GeV/c2  for 
t  quark.  Since  u  and  d  are  very  light,  the  isospin  SU(2)  symmetry  is  a  very  good 
approximation.  In  the  limit  that  isospin  is  an  exact  symmetry,  the  u  and  d  quarks 
are  considered  identical  except  the  charge  difference. 

Analogous  to  the  photon  in  electromagnetic  interactions,  a  strong  coupling  hap- 
pens via  a  mediating  boson,  call  a  gluon.  In  the  theory  of  quantum  chromodynamics 
(QCD),  there  are  six  types  of  strong  charge,  called  "color,"  representing  an  internal 
degree  of  freedom.  A  quark  can  carry  one  of  the  three  red,  blue,  or  green  colors, 
and  an  antiquark  can  carry  one  of  the  corresponding  anticolors.  The  quark-quark 
force  is  independent  of  quark  color,  so  the  color  symmetry,  SU(3),  is  exact. 

All  baryons  are  three-quark  (qqq)  states,  and  each  such  state  is  an  SU(3)  color 
singlet,  a  completely  antisymmetric  state  of  the  three  possible  colors.  Since  quarks 
are  fermions,  the  state  function  for  any  baryon  must  be  antisymmetric  under  inter- 
change of  any  two  equal-mass  quarks  (up  and  down  quarks  in  the  limit  of  isospin 
symmetry).  The  state  function  of  a  baryon  can  be  written  as 

\qQQ)a  =  \c°l°r)A  x  lsPOce,  spin,  flavor) s, 

where  the  subscripts  S  and  A  indicate  symmetry  or  antisymmetry  under  interchange 
of  any  two  of  the  equal-mass  quarks. 

For  baryons  made  up  of  u,  d,  and  s  quarks,  the  three  flavors  imply  and  approxi- 
mate flavor  SU(3).  This  flavor  SU(3)  symmetry  requires  that  a  corresponding  group 
of  baryons  belong  to  the  multiplets  on  the  right  side  of 

3  ®  3  o  3  =  ios  e  8M  e  8M  e  1.4, 

*The  values  of  quark  masses  depend  on  precisely  how  they  are  denned.  The  quark  masses  we 
referred  to  here  follow  the  meanings  defined  by  the  Particle  Data  Group  [1],  which  are  different 
from  the  constituent  quark  masses. 
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where  the  subscript  M  indicates  mixed-symmetry  states  under  interchange  of  any 
two  quarks.  Taking  spin  symmetry,  SU(2),  into  account,  the  two  8'Ms  will  end  up 
with  one  octet  of  baryon  states  of  \spin,  flavor) s.  This  octet  is  shown  as  the  bottom 
layer  of  Figure  1.1(a)  In  the  ground  state  of  the  multiplet,  the  SU(3)  flavor  singlet  1 
is  forbidden  by  Fermi  statistics.  The  decuplet  formed  by  105,  shown  as  the  bottom 
layer  in  Figure  1.1(b),  contains  the  states  with  the  three  quark  spins  aligned  to  the 
same  direction  (spin  J  =  3/2). 

If  we  add  the  c  quark  to  the  three  light  quarks,  the  flavor  symmetry  is  extended 
to  SU(4).  The  SU(4)  multiples  are  shown  in  Figure  1.1(a)  and  Figure  1.1(b).  All 
the  particles  in  a  given  multiplet  have  the  same  spin  and  parity.  Since  the  mass  of 
the  c  quark  is  large,  this  SU(4)  symmetry  is  very  badly  broken. 

1.2    Heavy-Quark  Effective  Theory 

1.2.1  QCD 

In  QCD,  the  Lagrangian  describing  the  interactions  of  quarks  and  gluons  is 

Lqcd  =  --fijF^+z^l^D,)^  -  £m^>*.  (1.1) 

9  q 

=  %K  -  dvAl  +  gsfabcAlAi, 

a  1 

where  gs  is  the  QCD  coupling  constant  (as  =  g*/4n),  and  the  fabc  are  the  structure 
constants  of  the  SU(3)  algebra.  The  (j)q(x)  are  the  4-component  Dirac  spinors  asso- 
ciated with  each  quark  field  of  color  i  and  flavor  q,  and  the  A*(x)  are  the  Yang-Mills 
gluon  fields. 


Figure  1.1:  SU(4)  multiplet  of  baryons  made  of  u,  d,  s,  and  c  quarks  shown  in 
(IZ,S,C)  coordinates,  (a)  The  20-plet  with  an  SU(3)  decuplet.  (b)  The  20-plet 
with  an  SU(3)  octet. 
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QCD  is  believed  to  be  the  fundamental  theory  of  strong  interactions.  But  to 
date,  there  are  no  practical  means  to  do  full  QCD  calculations  of  hadron  masses 
and  their  decay  widths,  and  a  variety  of  different  approximate  methods  have  been 
introduced.  Potential  models  have  been  well  recognized  to  be  the  most  useful  in 
understanding  the  mass  spectrum  and  the  decay  patterns  of  heavy-flavor  hadrons. 
The  correlation  between  flavor  and  spin  wave  functions  is  outside  the  realm  of  QCD. 
Some  assumptions  can  be  made  using  the  empirical  evidence  of  the  constituent  quark 
model.  Heavy-hadron  perturbation  theory  can  be  used  to  predict  the  decay  widths. 

1.2.2    Potential  Models 

Although  the  details  of  modern  potential  models  have  evolved  greatly  since  the 
discovery  of  J /(f),  the  essential  features  remain  unchanged.  The  potential  between 
two  quarks  is  often  taken  to  be  of  the  form 

Vs  =  --—  +  kr, 
3  r 

where  r  is  the  inter-quark  separation,  and  as  ~  1,  about  two  orders  larger  than 
a,  represents  the  magnitude  of  the  strong  coupling.  The  factor  4/3  holds  for  the 
interaction  between  a  quark  and  an  antiquark.  It  is  2/3  for  the  interaction  between 
two  quarks.  The  first  term  dominating  at  small  distance  r  arises  from  a  single 
gluon  exchange,  and  the  linear  term  is  associated  with  the  confinement  of  quarks 
and  gluons  inside  hadrons  at  large  r.  Because  of  the  linear  term,  attempts  to  free 
a  quark  from  a  hadron  simply  result  in  the  production  of  new  qq  quark-antiquark 
pairs.  The  annihilation  process  e+e"  ->  hadrons  is  viewed  in  terms  of  the  process 
e+e~  ->  qq  followed  by  "fragmentation"  of  the  quark  and  antiquark  into  hadrons. 
One  important  feature  of  potential  models  is  that,  except  for  mass  dependencies, 
the  inter-quark  potential  is  independent  of  quark  flavor. 
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Potential  models  cannot  be  formally  derived  from  QCD,  but  they  can  be  in- 
terpreted by  it.  In  QCD,  the  short-distance  behavior  is  dominated  by  one  gluon 
exchange  described  in  the  lowest-order  perturbation  theory.  The  linear  behavior  at 
large  r  can  be  seen  by  considering  that  a  qq  pair  are  attracted  to  each  other  and 
a  color  flux  tube  forms  between  them.  Lattice  calculations  suggest  that  the  ten- 
sion per  unit  length  in  this  flux  tube  is  a  constant[2,  3],  and  therefore  the  energy 
stored  in  the  tube  is  proportional  to  the  separation  of  the  charges.  Finally,  couplings 
described  in  QCD  are  flavor  independent. 

The  most  well-known  and  successful  potential  model  is  described  Isgur,  et  al. 
[4,  5,  6].  Early  models[5,  6]  are  non-relativistic  and  later  they  are  improved  by 
adding  some  relativistic  terms[4].  The  Hamiltonian  of  a  baryon  is, 

H  =  EK  +  |H  +        +  hibfI  (i.2) 

i=l  *mi  i<j 

3  r 

HliF  =  atij  ■  ^  +  cjttj  ■  fj  +  cr(3^  •  rtj  r-~$i-  fj)  +  cHti  ■  % 
where  i,  j  run  over  all  quarks  inside  the  baryon,  and 

,1         1     4as  lb 
m2k     mim/3r3     2m2kr  J 

1  Aas 

ct  =  r-r, 

miTrij  3  rd 

87r    1    4as  3  . 
3  raj/ft,  3  r 

Hlc3  is  the  confinement  term,  and  the  spin-dependent  term  H%F  is  very  similar  to 
the  Breit-Fermi  Hamiltonian  of  atomic  and  nuclear  physics.  It  is  very  crucial  that 
in  Hc  the  as  and  b  terms  have  opposite  signs,  and  the  two  contributions  to  the 
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spin-orbit  coupling  cancel  almost  exactly.  Since  this  model  does  not  perform  a  full 
relativistic  treatment,  it  is  called  a  "relativized"  quark  potential  model.  Compared 
with  non-relativistic  models,  it  has  slightly  different  values  of  a  and  b. 

In  the  limit  of  mQ  ->  oo,  where  Q  indicates  the  heavy  quark,  the  Hamiltonian 
is  simplified  and  a  new  good  quantum  number  j  the  total  angular  momentum  of 
the  light  degrees  of  freedom,  is  introduced,  j  ,  is  the  direct  analogy  with  the  total 
angular  momentum  of  the  electron  in  a  hydrogen  atom. 

1.2.3    Heavy-Quark  Chiral  Perturbation  Theory 

The  quark  contribution  to  the  QCD  Lagrangian  in  equation  1.1  can  be  considered 
to  have  two  parts.  The  first  comes  from  the  light  degrees  of  freedom,  and  the  second 
part  is  from  heavy  quarks.  Each  of  the  two  parts  has  distinct  symmetry. 

The  light-quark  part  of  the  Lagrangian  has  a  flavor  SU(3)lxSU(3)h  chiral  sym- 
metry in  the  limit  that  the  light-quark  masses  are  set  to  zero.  SU(3)LxSU(3)fi  is 
only  an  approximate  symmetry  since  the  quark  masses  explicitly  break  it.  This  chi- 
ral symmetry  is  spontaneously  broken  and  leads  to  eight  massless  Goldstone  bosons, 
the  7r's,  AT's,  and  77.  Their  couplings  to  hadrons  are  determined  by  PCAC  (partial 
conservation  of  axial-vector  current)  and  current  algebra,  or,  alternatively,  by  the 
nonlinear  chiral  Lagrangians. 

For  the  heavy-quark  part,  in  the  limit  of  infinite  quark  masses,  the  dynamics  of  a 
heavy  quark  in  QCD  depends  on  its  velocity  and  is  independent  of  its  mass  and  spin. 
As  a  consequence,  a  new  flavor  and  spin  symmetry  appear  for  hadrons  containing 
one  heavy  quark.  This  symmetry  is  called  heavy-quark  symmetry.  Many  models 
based  on  different  approaches  have  been  proposed  by  theorists  from  [7,  8,  9,  10,  11, 
12,  13,  14,  15,  16].  This  new  symmetry  implies  that  the  excitation  spectrum  and 
the  transition  form  factors  in  weak  decays  of  heavy  hadrons  are  independent  of  the 
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heavy-quark  species  and  heavy-quark  spins.  This  is  similar  to  a  hydrogen-like  atom 
in  QED  having  an  excitation  spectrum  and  transition  matrix  elements  independent 
of  the  mass  and  spin  of  the  nucleus. 

Since  heavy  hadrons  contain  both  heavy  and  light  quarks,  both  the  chiral  symme- 
try of  the  light  quarks  and  heavy-quark  symmetry  affect  the  low-energy  dynamics 
of  heavy  hadrons  interacting  with  the  Goldstone  bosons.  Experimentally,  strong 
decays  such  as  Ec  -»  A+7r  involving  soft  pions  give  such  examples. 

1.3    Introduction  to  Charmed  Baryons  decaying  to  A.+  (nn) 

1.3.1  Spectroscopy 

The  charmed  baryons  which  are  covered  in  this  work  are  those  containing  a  c 
quark  and  two  u  or  d  quarks.  The  two  light  quarks  are  often  referred  to  as  a  diquark. 
In  the  mc  ->  oo  limit,  the  total  angular  momentum  of  the  light  constituents, 

~fi  =      +  ~t, 

become  a  good  quantum  number.  Here  is  the  total  spin  of  the  two  light  quarks 
and  ~t  is  the  orbital  angular  momentum  between  the  c  quark  and  the  two  light 
quarks.  In  potential  models,  j  l  has  an  analogy  with  the  total  angular  momentum 
of  an  electron  in  a  hydrogen  atom.  Therefore,  the  properties  of  a  charmed  baryon  are 
mostly  dictated  by  its  light-quark  contents.  In  HQET,  the  final  angular  momentum 
of  the  baryon  is  given  as 

where  ^c  is  the  spin  of  the  c  quark.* 

*For  L  =  0,  we  simply  have  7  =  5^,  +  ^c.  But  when  L  is  non-zero,  the  angular  momenta  can 
be  coupled  as  either  7  =  (Sc  +  S,)  +  £  or  7  =  ~SC  +  (3,  +        Although  the  final  physics 
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Table  1.2:  The  s-wave  (both  L  and  L,  are  0)  and  p-wave  (either  L  or  Lt  is  1) 
charmed  baryons  and  their  quantum  numbers:  isospin  (/),  orbital  angular  momen- 
tum between  the  diquark  and  the  c  quark  (L),  orbital  angular  momentum  between 
the  two  light  quarks  (L/),  total  spin  of  the  diquark  (5/),  total  angular  momentum 
of  the  light  degree  of  freedom  (j,),  and  spin-parity  of  the  baryon  (Jp).  5/  and  jt  are 
good  quantum  numbers  in  the  limit  of  mc  — >  oo. 


State 

I 

L 

u 

Si 

./'' 

Ac+ 

0 

0 

0 

0 

0 

1+ 
2 

£c 

1 

0 

0 

1 

1 

1+  3+ 
2  '  2 

Aci 

0 

1 

0 

0 

1 

1-  3- 
2  >  2 

SC0 

1 

1 

0 

1 

0 

1 - 
2 

£C1 

1 

1 

0 

1 

1 

1-  3- 
2  '  2 

2C2 

1 

1 

0 

1 

2 

3-  5- 
2  '  2 

A'c0 

0 

0 

1 

1 

0 

1- 

2 

Ki 

0 

0 

1 

1 

1 

1-  3- 
2  '  2 

0 

0 

1 

1 

2 

3-  5- 
2  '  2 

1 

0 

1 

0 

1 

1-  3- 
2  '  2 

The  properties  of  all  L  =  0  (s-wave)  and  L  =  1  (p-wave)  baryons  of  cim,  cud, 
and  cdef  states  are  listed  in  Table  1.2.  The  lowest-lying  L  =  0  state  is  the  A+,  a  cud 
state  in  which  the  ud  diquark  has  isospin  of  /  =  0.  The  state  is  antisymmetric  under 
the  exchange  of  the  u  and  d  quarks.  Since  it  is  required  that  the  total  \spin,  flavor) 
be  symmetric,  the  spin  state  has  to  be  antisymmetric.  The  the  ud  diquark  must  be 
in  a  Si  =  0  state.  Combining  Si  with  the  spin  of  the  c  quark  gives  a  total  spin  of 
J  =  1/2.  By  convention,  the  A+,  like  all  spatial  ground  state  baryons,  has  positive 
parity. 

For  Ec  states,  which  have  /  =  1,  the  diquark  must  have  5/  =  1  to  satisfy  the 

requirement  of  \spin,  flavor)  being  symmetric.  Adding  S,  =  1  with  the  spin  of  the 

does  not  depend  on  the  order  of  coupling,  the  latter  choice  leads  to  a  larger  amount  of  algebra  in 
the  calculations  of  most  of  the  HQET  models.  The  choice  of  ~f  =  !)c  +  j  ,  coupling  also  helps  us 
derive  the  physical  states  more  easily. 
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c  quark  results  a  doublet  having  Jp  =  1/2+  (Ec)  and  Jp  =  3/2+  (E*).  A+,  Ec,  and 
E*  are  called  s-wave  particles  since  there  is  no  orbital  angular  momentum  involved. 

It  is  worthing  mentioning  that  we  have  the  similar  L  =  0  states  in  the  strange 
baryon  sector,  A,  E,  and  E*.  An  important  fact  is  that 

M(EC*)  -  M(EC)  ms 
M(E*)  -  M(E)  ~  mc' 

for  the  hyperfine  splittings.  Based  on  HQET,  a  hyperfine  splitting  arises  from  the 
color  magnetic  moment  of  the  heavy  quark.  Since  ms/mc  ~  1/3,  the  M (E*)-M(EC) 
is  expected  to  be  about  one-third  of  M(E*)  -  M(E). 

Now  Let  us  consider  the  lowest-lying  L  =  1  (p-wave)  particles.  In  the  /  =  0 
case,  there  must  be  S  =  0  for  the  same  reason  as  for  the  A+  baryon,  so  jt  =  1. 
The  c  quark  brings  hyperfine  splitting  to  give  two  baryon  states  of  Jp  =  l/2~  and 
Jp  =  3/2",  where  the  negative  parity  comes  from  the  one  unit  of  the  orbital  angular 
momentum.  This  doublet  is  named  Acl. 

Same  as  Ec,  for  /  =  1,  we  have  S  =  1.  This  gives  j,  =  0, 1, 2.  Each  of  these  is 
split  by  its  hyperfine  interaction  with  the  c  quark.  This  type  of  excitation  ends  up 
with  five  baryons:  the  Ec0  singlet,  the  Ecl  doublet,  and  the  Ec2  doublet. 

Another  type  of  p- wave  states  come  from  one  unit  of  the  orbital  angular  momen- 
tum between  the  two  light  quarks,  Lt  =  1.  This  implies  that  the  /  =  0  ud  diquark 
must  have  5  =  1.  This  produces  jt  =  0, 1,2.  Then  the  hyperfine  splittings  give  five 
baryon  states:  the       singlet,  the  A'cl  doublet,  and  the  A'c2  doublet. 

Finally,  it  can  be  easily  derived  that  U  =  1  and  1  =  1  require  5  =  0.  Then  the 
hyperfine  splitting  gives  a  doublet  of  E'cl  baryons.  Note  that  the  degeneracy  of  the 
L  =  1  (Lt  =  1)  and  Lt  =  1  (L  =  1)  switches  with  the  change  from  /  =  0  to  /  =  1. 
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1.3.2    Strong  Decays 

All  the  charmed  baryons  we  described  in  the  last  section  will  decay  to  the  A+, 
since  it  is  the  lowest  state  of  it,  d,  c  compositions.  The  purpose  of  these  research  work 
is  to  get  better  understanding  on  how  the  baryons  and  even  p-wave  charmed 
baryons  decay  to  A+. 

To  date,  only  the  Ec  triplet,  two  states  of  the  E*  triplet,  and  the  Aci  doublet  have 
been  found,  and  only  strong  decays  through  one  or  two  pions  have  been  observed. 
Based  on  the  mass  predictions  by  the  relativized  potential  model[4],  other  p-wave 
baryons  may  decay  strongly  to  the  final  states  containing  a  A+  and  one,  two,  or 
three  pions,  but  this  work  failed  see  any  more  p-wave  states  higher  than  the  Acl. 
In  general,  the  higher  states  which  have  the  same  decay  modes  are  expected  to  be 
wide  and  difficult  to  observe.  So  this  dissertation  only  reports  the  studies  of  the 
E^  and  Acl  baryons.  This  subsection  gives  a  brief  review  of  the  strong  decays  of 
these  baryons. 

The  selection  rules  of  these  decay  can  be  easily  derived  from  conservation  of 
parity  and  angular  momentum.  For  a  strong  decay  A  ->■  BC,  where  the  spin- 
parities  of  these  particles  are  denoted  by  JPAA ,  J^B ,  and  jfP,  and  a  possible  orbital 
angular  momentum  between  B  and  C  is  denoted  by  LBC,  it  is  required  that  L  must 
satisfy 

A  =  L  BC  +  ~3  B  +  ~~?Ci 

and 

Pa  =  PbPc{-1)Lbc 

in  order  to  decay  with  a  partial  wave  corresponding  to  LBc-  In  addition,  isospin 
has  to  be  conserved  in  all  strong  decays. 
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Table  1.3:  Allowed  partial  waves  for  strong  decays  of  the  E^  and  Aci  baryons.  The 
A+7T7T  mode  listed  in  the  table  only  represents  the  non-resonance  decay  mode. 


Particle 

Decav  Mode 

Allowed  Partial  Wave 

sc 

A+7T 

P 

E* 

A+w 

P 

Alu/2) 

A+7T7T 

P 

Acl(l/2) 

EC7T 

S 

Aci(3/2) 

A+7T7T 

P 

Acl(3/2) 

EC7T 

D 

Table  1.3  lists  the  all  the  possible  decay  channels  for  strong  decays  of  the 
and  Acl  baryons.  The  E^  baryons  do  not  have  the  A+7T7T  decay  because  there  is  no 
phase  space.  They  should  entirely  decay  through  A+7T.  The  Aci  baryons  do  not  have 
the  A+7r  decay  because  of  isospin  conservation.  Based  on  the  table,  we  should  expect 
that  the  S- wave  Ec7r  would  dominate  the  Aci(l/2)  decay.  The  Acl(3/2)  can  only 
have  D-wave  Ec7r  decay  because  of  the  conservation  of  total  angular  momentum; 
the  Acl(3/2)-»E*7r  decay  would  be  S-wave,  but  is  kinematically  forbidden.  Thus 
we  should  expect  that  P-wave  non-resonance  A+7T7T  decay  dominates  the  Acl(3/2) 
decay. 


1.4    HQET  Predictions  on  £<*)  and  A+x 

Many  quark  models  have  the  ability  to  make  quantitative  predictions  on  heavy- 
baryon  masses,  and  almost  all  of  such  models  use  the  "constituent  quark  masses" 
which  are  derived  from  well-known  hadron  masses.  Predictions  vary  from  model  to 
model,  and  they  all  fall  in  roughly  the  same  range.  Among  all  the  models,  the  quark 
potential  model  developed  by  Isgur  et  al.[4,  6]  is  the  most  complete,  and  it  covers 
all  the  baryons  from  the  proton  to  bottom  baryons. 
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Table  1.4:  The  constituent  quark  mass  assumptions  by  the  potential  models.  All 
numbers  are  in  units  of  MeV/c2. 


Quark 

Non-relativistic  Model 

Relativized  Model 

u 

350 

220 

d 

350 

220 

s 

550 

419 

c 

1500 

1628 

b 

5000 

4977 

Table  1.5:  Some  charmed  baryon  mass  predictions  by  the  potential  models  in  com- 
parison with  average  values  of  experiments.  All  numbers  are  in  units  of  MeV/c2. 
All  uncertainties  of  the  experiments  measured  numbers  are  less  than  3  MeV/c2. 


Quark 

Non-relativistic  Model 

Relativized  Model 

Experiment 

2260 

2265 

2285 

sc 

2440 

2440 

2453 

K 

2510 

2495 

2518 

Ac+i(l/2) 

2510 

2630 

2594 

A+(3/2) 

2590 

2640 

2627 

An  important  part  of  the  potential  models  is  to  specify  the  constituent  masses  of 
quarks.  Typical  values  used  in  both  non-relativistic  models  and  relativized  models 
are  listed  in  Table  1.4.  Although  they  use  different  mass  assumptions,  both  models 
work  extremely  well  for  light  baryons  and  strange  baryons,  and  their  predictions 
on  charmed  baryons  are  also  close.  Table  1.5  lists  some  of  the  predictions[4,  6] 
compared  with  the  average  experimental  values[l]. 

Potential  model  also  predicted  the  masses  of  other  p-wave  charmed  baryons. 
The  masses  of  Ec0  (the  next  state  higher  than  Aci)  and  Ecl  doublet  are  predicted 
to  lie  around  2770  MeV/c2,  within  10  MeV/c2.  Other  nine  states  are  also  close 
to  each  other  in  masses  and  range  from  2780  MeV/c2  to  2900  MeV/c2.  Based  on 
these  predictions,  unless  some  of  these  states  have  very  narrow  width,  it  would  be 
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difficult  to  see  them.  In  our  work,  we  looked  at  these  mass  ranges  for  A+7r±  and 
A+7r+7r" decays,  and  did  not  see  any  obvious  signals. 

The  potential  models  can  make  calculations  on  decay  widths  but  do  not  work 
very  well.  Their  predictions  are  very  rough  and  agree  with  experimental  values  only 
in  their  order  of  magnitude. 

The  most-recognized  theoretical  model  which  can  predict  heavy-hadron  widths 
was  developed  by  Yan  et  al.[16],  who  used  the  nonlinear  chiral  quark  model[17] 
to  calculate  heavy-quark  coupling  to  the  Goldstone  bosons.  In  a  charmed  baryon, 
the  diquark  with  sf'  =  1+  can  be  represented  by  an  axial- vector  field  and  the 
diquark  with  sf'  =  0+  can  be  represented  by  a  Lorentz  scalar  field  0.  The  decays 

~^  ^+7r  anc*  ^  ~^  ^+7r  can  be  described  by  two  independent  coupling  constants 
9\  and  g2,  respectively.  For  the  decay  £c-»A+7r,  the  decay  rate  is  determined  by  g2 
and  masses  of  the  £c,  A+,  and  tt.  Their  mode  predicts 

2 

92  =  -gfl 'A, 

where  gA  is  the  coupling  constant  in  the  single-quark  transition  u  ->  d.  gA  is 
generally  assumed  to  have  a  value  of  1,  but  the  experimental  extracted  value  is 
0.75[17].  g2  values  are  not  very  different  from  other  theoretical  predictions[18,  19, 
20].  Their  Lagrangian  gives  the  following  decay  width: 

where  fw  =  93  MeV/c2  is  the  pion  decay  constant  and  pv  is  the  pion  momentum 
in  the  center  of  mass  frame,  equation  1.3  predicts  the  total  width  of  E°  since  Ec 
decays  100%  to  A+tt.  equation  1.4,  however,  is  invalid  since  the  decay  E*->Ectt  is 
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kinematically  forbidden.  If  we  neglect  the  radiative  decay  width  associated  with  the 
decay  E*— >EC7,*  the  total  width  of  E*  is  given  by 

g2  =  0.612  (with  gA  =  0.75)  will  lead  to 


T(EC)  =  2.45  MeV/c2, 

=  17.6  MeV/c2,   if  Ms.  -  MA+  =  233.5  MeV/c2. 

Equations  1.3  and  1.5  show  that  the  Ec  width  is  directly  related  to  the  E*  width, 

which  has  been  reported  by  CLEO  collaboration  [21].  This  work  includes  the  first 

attempt  to  measure  that  Ec  width.  If  these  widths  can  be  measured  accurately, 

g2  can  be  derived  from  equation  1.3  or  equation  1.4  and  then  further  predicts  the 

widths  of  other  particles,  for  example,  H*. 

Recently  Pirjol  and  Yan[22]  included  all  possible  strong-interaction  couplings 

(S'-wave,  P-wave,  and  D-wave)  among  and  between  s-wave  and  p-wave  baryons  in 

the  chiral  Lagrangian.  This  model  contains  45  independent  coupling  constants  up 

to  and  including  D-wave  interactions.  Besides  gi  and  g2,  coupling  constants  h2, 

/i3,  hi,  ...  are  involved,  which  correspond  to  the  p-wave  baryons.  In  this  theory, 

they  derived  model-independent  sum  rules  which  contain  these  couplings  and  relate 

them  to  properties  of  the  lowest-lying  baryons.    However,  couplings  h2  and  h$, 

which  correspond  to  lowest  p-wave  baryons,  cannot  be  directly  predicted  from  the 

constituent  quark  model.    Using  published        masses[l]  M(A+1(2593)  -  A+)  = 

*A  radiative  decay  can  happen  when  the  strong  decay  is  largely  suppressed.  Since  the  widths 
of  scst->Ec7r  decays  have  been  measured  to  be  around  16  MeV/c2,  there  should  be  almost  no 
radiative  decay. 
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308.6  MeV/c2  and  M(A+!(2625)  -  A+)  =  341.5  MeV/c2  together  with  g2  =  0.57,* 
and  neglecting  radiative  modes  A^j  — ►  A+7,  they  calculated 

r(A+(2593))  =  11.9/i2  +  13.8^  -  0.042h2h8  (MeV/c2), 

r(A+(2625))  =  0.52/12  +  (0.15  x  106)/i2  -  5.2h2h8  (MeV/c2). 

With  r(A+!(2593))  and  r(A+j(2625)),  the  values  of  h2  and  h8  can  be  obtained, 
and  they  can  be  used  to  further  estimate  the  width  of  other  particles.  Based  on  the 
published  values  of  r(A+!(2593))  and  r(A+1(2625))[l],  it  can  be  estimated  that 

h2  =  0.57±°0%, 

hs  <  3.5  x  10~3  (MeV/c2)"1. 

They  further  estimated  that  Ec0  and  Eclhave  widths  larger  than  100  MeV/c2,  If 
these  estimates  are  correct,  it  will  be  very  diffucult  to  observe  signals  of  next-level 
higher  state  baryons. 


•This  value  is  the  average  value  obtained  from  applying  published  masses[21]  M(S*++  -  A+)  = 
234.5  MeV/c2  and  M(S*°  -  A+)  =  232.6  MeV/c2  to  equation  1.5. 


CHAPTER  2 
EXPERIMENTAL  FACILITIES 


2.1  CESR 

The  Cornell  Electron  Storage  Ring  (CESR)  is  an  electron-positron  collider  with 
a  circumference  of  768  meters  located  on  the  campus  of  Cornell  University,  Ithaca, 
NY.  It  can  produce  collisions  between  electrons  and  positrons  with  center-of-mass 
energies  between  9  and  12  GeV/c2.  It  serves  both  CLEO  for  the  study  of  particle 
physics,  and  the  Cornell  High  Energy  Synchrotron  Source  (CHESS)  for  a  variety  of 
biological  and  surface  physics  studies. 

A  diagram  of  CESR's  main  components  is  shown  in  Figure  2.1.  This  accel- 
erator consists  of  three  major  parts:  (1)  the  linear  accelerator  (LINAC),  (2)  the 
synchrotron,  and  (3)  the  storage  ring.  Electrons  are  accelerated  to  150  MeV/c2  by 
the  150-foot  LINAC.  Positrons  are  created  by  50  MeV/c2  electrons  colliding  with  a 
thin  tungsten  target  part-way  through  the  LINAC.  The  electrons  and  positrons  are 
then  boosted  to  the  operating  energy  of  about  5  GeV/c2  by  the  synchrotron  and 
transferred  into  the  storage  ring,  where  they  can  be  maintained  in  the  storage  ring 
for  about  1  hour  in  order  to  achieve  the  highest  integrated  luminosity. 

Electrons  and  positrons  travel  in  opposite  directions  around  the  storage  ring 
in  evenly  spaced  bunches  (currently  there  are  9)  at  390,000  revolutions  per  second. 
Electrostatic  separators  hold  the  electron  beam  and  the  positron  beam  slightly  apart 
from  each  other,  and  two  beams  only  collide  in  one  place  —  the  center  of  the  CLEO- 
II  detector. 
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O   Positron  Bunch  -  Clockwise 


#    Electron  Bunch  -  Counter  Clockwise 


Figure  2.1:  CESR  layout.  The  CESR  lies  in  a  tunnel  50  feet  underground.  CLEO-II 
and  the  LINAC  reside  in  Wilson  Lab. 
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2.2  CLEO-II 

CLEO-II  is  a  multipurpose  high  energy  physics  detector  incorporating  excellent 
charged  and  neutral  particle  detection  and  measurement  used  to  analyze  electron- 
positron  collision  events  generated  by  CESR.  It  is  operated  by  the  CLEO  collab- 
oration of  over  200  physicists  from  many  institutions,  including  the  University  of 
Florida. 

Side  and  end  views  of  the  detector  are  shown  in  Figure  2.2  and  Figure  2.3,  respec- 
tively. The  major  components  (going  outward  from  the  beam  pipe)  are  the  central 
detector  (CD)  which  forms  the  tracking  system,  the  time  of  flight  (TOF)  system,  the 
electromagnetic  calorimeter  (CC),  the  1.5-Tesla  superconducting  magnet,  and  muon 
chambers.  The  central  detector  comprises  the  Precision  Tracking  Layer  (PTL),  the 
Vertex  Detector  (VD),  and  the  Drift  Chamber  (DR).  CLEO-II  is  also  equipped  with 
timing,  trigger,  and  data  acquisition  systems.  Brief  description  of  the  CD,  TOF, 
and  CC  are  given  in  the  following  since  they  are  the  most  important  in  this  analysis. 

2.2.1    Tracking  System 

Charged  particle  momentum  is  measured  by  three  cylindrical,  coaxial  wire  cham- 
bers sharing  a  common  axis  in  the  direction  of  the  beams. 

Figure  2.4  shows  the  structure  of  PTL  and  VD.  The  PTL  is  a  cell-strawtube- 
drift-chamber  which  extends  within  1  cm  of  the  3.5  cm  beam  pipe.  The  device 
consists  of  6  layers  of  mylar  tubes  with  64  axial  wires  per  layer.  It  is  used  to  make 
precise  measurement  of  the  transverse  position  of  the  particle  near  the  interaction 
point,  with  a  resolution  of  90  //m;  and  it  does  not  provide  longitudinal  information. 
The  VD  covers  the  radial  region  from  7.5  to  17.5  cm.  It  has  10  layers  of  hexagonal 
cells  as  shown  in  Fig.  2.4.  arranged  in  10  layers.  The  VD  provides  a  tracking 
resolution  of  150  /zm  in  r  -  (j>  and  0.75  mm  in  z. 


Figure  2.2:  Side  view  of  the  CLEO-II  Detector. 
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Muon  Chambers 


Magnetic  Coil 
Barrel  Crystals 
Time  of  Flight 


Beam  Pipe 
PTLA  and  VD 


Figure  2.3:  End  view  of  the  CLEO-II  detector. 
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Figure  2.4:  The  VD/PTL  wires. 


•  Field  Wires 

o  Stereo  Sense  Wires 

•Axial  Sense  Wires 


DR  inner  Cathodes 
DR  Inner  Shell 
VD  Outer  Shell 


Figure  2.5:  The  DR  wires. 
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Transverse  and  longitudinal  momentum  are  best  measured  by  DR,  shown  in 
Figure  2.5.  This  51-layer  outer  drift  chamber  extends  from  17.5  cm  to  95  cm, 
including  12,240  sense  wires  and  36,240  wires.  One  sense  wire  and  three  field  wires 
make  a  rectangular  cell  structure.  There  are  40  axial  (parallel  to  the  beam  line) 
layers  and  11  stereo  (not  parallel  to  the  beam  line)  layers.  Axial  layers  measure 
transverse  momentum,  the  radial  distance  of  the  closest  approach  to  the  beam  line, 
and  the  azimuthal  direction  of  tracks.  The  stereo  layers  and  2  cathode  layers  measure 
the  polar  angle  and  the  longitudinal  position  of  the  intersection  of  the  track  with 
the  beam  axis.  Sequential  layers  are  offset  in  azimuth  in  order  to  resolve  left-right 
ambiguity  in  drift  distance.  An  r  -  (j)  position  resolution  of  110  /im  for  the  axial 
wires  and  a  z  position  resolution  of  about  3  cm  can  be  obtained  with  this  chamber. 

Particles  can  be  identified  from  the  specific  ionization  energy  loss  (dE/dx)  in 
the  DR.  The  characteristic  bands  of  dE/dx  quantity  as  a  function  of  momentum  for 
different  species  of  particles  are  shown  in  Figure  2.6.  Protons  can  be  well  separated 
from  other  particles  at  momenta  below  1  MeV/c,  and  separation  between  pions  and 
kaons  is  possible  only  below  400  MeV/c. 

Track-reconstruction  programs  are  used  to  find  a  charged-particle  track  in  an 
event.  First,  the  pattern  of  hits,  which  are  the  records  of  the  wirescaused  by  the 
ionized  gas  molecules,  are  recognized;  then  these  hits  are  fitted  to  obtain  the  trajec- 
tory of  the  original  particle. 

The  data  from  all  three  chambers  are  combined  to  measure  the  momentum  vector 
of  charged  particles.  Two  major  factors  limit  the  momentum  resolution  —  multiple 
scattering  and  the  position  resolution  of  the  track.  Multiple  scattering  dominates 
the  resolution  at  low  momentum,  and  the  position  resolution  dominates  at  high 
momentum,  where  the  track  curvature  is  small. 


Figure  2.6:  dE/dx  plotted  against  momentum  for  different  particles. 
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The  CLEO-II  charged-track  momentum  resolution  has  been  measured  to  be 

8p/p  =  yo.0052  +  (0.0015)  V 

where  the  momentum  p  is  in  GeV/c.  The  first  term  represents  the  contribution 
from  multiple  scattering,  and  the  second  term  is  the  contribution  from  the  curvature 
measurement  error.  The  two  terms  give  equal  contributions  at  a  momentum  of  3.3 
GeV/c. 

Charged  particles  with  a  transverse  momentum  greater  than  225  MeV/c  and  a 
polar  angle  less  than  45°  will  reach  the  outer  radius  of  the  DR  and  pass  through  the 
TOF  counters  and  crystal  calorimeter.  These  tracks  are  the  best  measured  in  the 
CLEO-II  detector. 

2.2.2    Time-of-Flight  System 

Besides  dE/dx  measurement  in  DR,  the  TOF  device  is  also  an  important  tool 
for  particle  identification.  The  TOF  system  has  two  major  parts:  a  barrel  system 
covering  polar  angles  from  36°  to  144°,  and  an  endcap  system  covering  polar  angles 
from  15°  to  36°  on  one  side  and  144°  to  165°  on  the  other  side. 

The  barrel  system,  which  is  located  immediately  outside  the  DR,  consists  of 
64  plastic  rectangular  scintillation  counters,  and  an  endcap  system  consists  of  the 
same  but  wedge-shaped  counters.  The  scintillators  are  5  cm  thick  and  2.5  m  long. 
Lucite  light  pipes  are  attached  at  each  end  of  the  scintillators  and  connected  to 
photomultiplier  tubes. 

It  takes  a  500  MeV/c  pion  about  3.5  ns  to  travel  from  the  interaction  region  to 
the  barrel  TOF  counters,  while  it  takes  a  kaon  about  4.5  ns.  The  time  resolution 
for  pions  in  hadronic  events  is  about  150  ps.  Figure  2.7  plots  1/0  {(5  =  v/c)  against 
momentum.  As  we  can  see  in  the  plot,  the  distinction  in  the  bands  separates  K-n 
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up  to  momenta  of  about  1.1  GeV/c,  and  protons  can  be  separated  up  to  momenta 
of  about  1.5  GeV/c2.  Like  the  energy  loss  information,  particles  can  not  be  well 
identified  by  TOF  at  high  momentum. 

2.2.3    Electromagnetic  Calorimeter 

The  Csl  crystal  calorimeter  is  used  to  detect  photons  and  to  measure  the  energies 
of  tracks.  It  is  composed  of  7800  thallium-doped  cesium  iodide  scintillating  crystals, 
6144  in  the  barrel  region  and  828  in  each  endcap,  covering  95%  of  the  solid  angle. 
Each  crystal  is  30  cm  long  and  5x5  cm2  in  cross  section. 

When  a  photon  or  a  charged  particle  interacts  with  the  atoms  in  the  crystal,  an 
electromagnetic  "shower"  is  emitted  by  crystals,  and  it  is  recorded  by  4  Si  photodi- 
odes  at  the  back  of  each  crystal.  The  device  allows  detection  of  photons  with  energy 
of  30  MeV/c2  to  5.3  GeV/c2  with  excellent  resolution. 

Photon  energy  resolution  in  the  barrel  (endcap)  is  1.5%  (2.6%)  at  5  GeV/c2,  and 
3.8%  (5.0%)  at  100  MeV/c2.  The  angular  resolution  in  azimuth  is  3  mrad  (9  mrad) 
at  5  GeV/c2  and  11  mrad  (19  mrad)  at  100  MeV/c2. 

A  typical  tt°  (from  the  decay  )  finding  efficiency  is  about  50%,  and  the  width  of 
the  reconstructed  7r°  peak  is  roughly  5-10  MeV/c2. 


CHAPTER  3 
A+  RECONSTRUCTION 


3.1  Introduction 

In  order  to  study  the  charmed  baryons  which  decay  to  A+'s,  we  first  need  to 
obtain  A+  candidates,  then  add  additional  pions  (one  or  two  in  our  study)  to  recon- 
struct heavier  charmed  baryons. 

Fifteen  hadron  decay  modes  were  used  to  reconstruct  A+  candidates.  They 
are  pl<-7r+,  pR-tt+tt0,  PK°,  pKVtt0,  pKV,  Att+,  Att+tt+tt",  Att+tt0,  AK+K°, 
E+7T+7T-,  E+7T°,  E+K+K-,  E°tt+,  ~-R+tt+,  and  H°K+.  All  these  decay  modes 
have  good  signal-to-noise  ratios  and  reasonable  statistics.  The  decay  pK_7r+  has 
a  very  large  branching  fraction  and  contributes  almost  half  of  all  A+  candidates. 
It  important  to  check  the  consistency  of  measurements  by  comparing  the  results 
obtained  by  using  the  pK~7r+  mode  with  those  obtained  from  other  modes. 


3.2    Data  and  Monte  Carlo  Sample 

CLEO  data  are  taken  on  and  off  the  T(45)  resonance.  The  beam  energy  is  5.187 
to  5.280  GeV/c2  for  off-resonance  and  5.285  to  5.295  GeV/c2  for  on-resonance.  We 
performed  our  analysis  using  both  on-resonance  and  off-resonance  events  in  datasets 
4S1-4SG,*  corresponding  to  a  luminosity  of  4.8  fb~\ 

All  Monte  Carlo  events  generated  were  with  the  A+  mass  2286.7  MeV/c2,  which 
is  our  measured  mean  value  over  all  the  15  decay  modes.   In  order  to  optimize 

•Each  time  CESR  or  CLEO  is  updated,  the  data  are  processed  differently,  and  they  become  a 
new  dataset.  Conventionally,  CLEO  datasets  are  named  4S1,  4S2,  4S9  4SA  4SB  and  so 
on.  "4S"  indicates  the  state  T(45). 
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Table  3.1:  Monte  Carlo  (MC)  generated  branching  fractions  for      and  E<*}  studies. 


UtiLay  IIlOQc 

Average  D/o^prv.  ix  ) 

MC  generated  branching  ratio 

i 
i 

n  i  no 

pK"7r+7r° 

n  fi7+n  f)44-0  1 1 

U.UI  XU.U4XU.  X  X 

U.  1UO 

pis. 

U.4o±U.Uz±U.U4 

A  M  —  O 

0.073 

pK  7T+7T 

0.52±0.04±0.05 

0.082 

0.66±0.05±0.07 

0.104 

Att+ 

0.18±0.03 

0.028 

A7r+7r+7r~ 

0.66±0.10 

0.104 

i           1  A 

A7T+7r0 

0.73±0.09±0.16 

0.115 

AJtS.  IS. 

C\  1  O  1  A  no  l  A  AO 

U.12±0.02±0.02 

0.019 

S+7T+7r- 

0.68±0.09 

0.107 

E+7T0 

0.20±0.03 

0.032 

E+K+K" 

0.07±0.01±0.01 

0.011 

E°7T+ 

0.20±0.04 

0.032 

E-K+7T+ 

0.10±0.02 

0.016 

E°K+ 

0.08±0.01±0.01 

0.013 

All 

1 

the  cuts*  used  to  reconstruct  each  decay  mode,  15  individual  types  of  events  are 

generated  by  sufficient  amounts  (more  than  100,000  events  for  each  mode).  In  our 

later  study  on  the  A+,  or  the  E^  decaying  to  the  A+,  we  need  the  Monte  Carlo  to 

provide  the  correct  over-all  efficiencies.  Thus  the  decay  ratios  of  A+  decay  modes 

were  generated  as  shown  in  Table  3.1.  In  the  table,  the  average  relative  branching 

ratios  for  pR-tt+tt0,  pK°,  pK°7r+7i-,  and  pKV  are  the  values  recently  measured 

by  CLEO[23].  For  other  modes  we  used  the  PDG  average  values[l]. 

*A  cut  is  an  action  to  remove  a  certain  fraction  of  the  events  which  contribute  mostly  to  the 
background.  An  optimized  cut  on  a  measured  physical  quantity  is  the  one  which  gives  the  best 
signal-to-noise  ratio  after  it  is  applied. 
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3.3    Event  Selection 

3.3.1    Track  Selection 

When  a  charged  particle  travels  through  the  tracking  chambers,  a  set  of  hits, 
caused  by  the  passage  of  the  track,  is  generated.  The  tracking  program  will  fit 
each  set  of  hits  forming  a  "track,"  which  is  a  curve  representing  a  possible  particle 
causing  the  hits.  In  CLEO,  this  task  is  accomplished  by  the  track-finding  and  fitting 
program,  DUET. 

Each  track  is  described  by  the  following  five  track  parameters:  (1)  CUCD:  cur- 
vature of  the  track;  (2)  FICD:  0  of  the  direction  of  the  track  at  point  of  closest 
approach  to  the  center  of  the  CLEO-II  detector;  (3)  DACD:  impact  parameter  of 
track  with  respect  to  the  center  of  the  detector;  (4)  CTCD:  cotangent  of  the  polar 
angle,  9  (i.e.  CTCD  =  0  for  tracks  in  the  r  -  0  plane  and  infinite  for  tracks  pointing 
down  the  beam  pipe;)  (5)  ZOCD:  z  coordinate  at  point  of  closest  approach  to  the 
center  of  the  detector.  All  other  information  about  the  track  can  be  calculated  from 
these  five  parameters. 

Before  1996,  the  track  fitter  in  the  CLEO  DUET  processor  ignored  the  energy 
loss  of  a  track  in  materials  and  included  multiple  scattering  approximately.  In  1996, 
the  CLEO  collaboration  reprocessed  the  data  from  4S1  to  4SG,  and  the  new  data 
are  called  the  recompressed  data.  For  the  recompressed  data,  the  old  track  fitter  is 
replaced  by  a  Kalman  Filter.  The  new  Kalman  fitter  optimally  treats  both  multiple 
scattering  and  energy  loss,  so  it  gives  improved  error  matrices  of  track  parameters 
for  all  tracks  and  better-measured  track  parameters  for  low-momentum  tracks.  It 
also  fits  each  track  separately  for  each  particle  hypothesis. 

To  select  a  track  as  a  candidate  for  a  pion,  a  kaon,  or  a  proton,  which  come  from 
the  r  -  0  range  of  e+e~  interaction  points,  we  first  require  it  to  have  not  only  r  -  0 
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also  z  information.  We  reject  all  tracks  which  are  identified  inward  going.  To  ensure 
these  tracks  are  close  to  the  interaction  point,  we  used  the  following  "track  quality" 
cuts: 


where  DBCD  is  the  radial  distance  of  the  point  of  closest  approach  to  the  beam  line, 
P  is  the  momentum  of  the  charged  track;  and 


Note  that  the  DBCD  cut  changes  from  5  mm  at  0  momentum  to  1.2  mm  at  1  GeV/c 
and  then  becomes  flat  from  1.2  mm.  This  is  because  the  low-momentum  tracks  are 
more  poorly  measured,  but  we  still  want  to  keep  those  having  fairly  good  quality. 

Such  selected  tracks  are  used  as  the  charged  track  candidates  which  are  directly 
decayed  from  A+'s.  The  A+  has  a  mean  lifetime  of  about  2  x  10"13s,  so  its  flight 
distance  in  the  CLEO  detector  is  within  the  range  of  interaction  points  and  is 
negligible  in  the  invariant  mass  reconstruction. 

3.3.2    Particle  Identification 

As  we  discussed  in  the  last  chapter,  two  types  of  detectors  are  used  for  the 
identification  of  the  charged  tracks:  the  time-of-flight  counters  and  the  drift  chamber 
for  the  dE/dx  measurement.  The  information  from  each  device,  if  present,  is  quoted 
as  the  number  of  standard  deviations  of  the  measurement  from  the  expected  value 
as  a  given  particle  hypothesis.  For  each  charged  track,  we  combine  the  TOF  and 
the  dE/dx  information  to  define  a  joint 


|DBCD|  <  < 


5  -  3.8|P|  (mm)   if  \P\  <  1, 


1.2  mm 


if|P|>  1, 


Z0CD|  <  5  cm. 


(dE/dx) 


meas 


(dE/dx)exp]2  \{TOF) 


meas 


(TOF) 


exp 


° dE/dx 


Otqf 


where  i  stands  for  pion,  kaon,  and  proton  hypothesis. 

The  x2  ranges  from  zero  to  infinity  and  it  is  not  very  convenient  to  use  when  we 
perform  a  track  identification  cut.  So,  a  equivalent  parameter,  the  \2  probability, 
PROB,  is  calculated  for  each  hypothesis.  PRQB  ranges  from  zero  to  one.*  For  data 
analysis,  it  has  the  advantage  of  being  flat  in  the  zero  to  one  range  for  particles 
which  are  the  right  species  and  peaking  at  zero  for  particles  which  are  the  wrong 
species. 

One  particle  identification  information  will  be  used  if  the  other  device  does  not 
have  the  measurement.  If  neither  TOF  nor  dE/dx  information  is  present,  we  assume 
the  track  is  a  pion  and  assign  it  a  small  PR0B(7r)  value. 

In  this  work,  we  are  more  concerned  about  separating  protons,  kaons  and  pions 
from  each  other,  so  a  normalized  probability,  called  "PROB  level,"  is  used.  To  strongly 
identify  a  proton  or  a  kaon,  we  make  a  cut  on 


PRLEV  =  -  PR0BW 


PROB(p)  +  PROB(A')  +  PROB(tt) 


or, 

PROB(A) 


KALEV  = 


PROB(p)  +  PROB(A)  +  PROB(tt)  ' 

We  do  not  use  a  PROB  level  cut  for  pions,  since  we  take  all  the  tracks  selected  as 
pion  candidates. 

3.3.3    7  and  ir°  Finding 

Although  the  Csl  electromagnetic  calorimeters  can  measure  the  photon  energy 
and  angle,  they  also  interact  with  hadrons.  An  electromagnetic  shower  generally 
has  different  lateral  development  than  a  hadronic  shower.  This  provides  us  a  way 
*Note  it  does  not  really  represent  the  probability  of  the  track  to  be  a  specific  particle  hypothesis. 
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to  distinguish  them  from  each  other.  The  quantity  E9/E25  is  measured  as 

Energy  of  9  central  crystals 


E9/E25  = 


Energy  of  25  central  crystals 


Here  the  9  crystals  in  the  numerator  are  those  in  a  3  x  3  array,  with  the  most 
energetic  crystal  at  the  center.  The  25  crystals  in  the  denominator  are  those  in  a 
5x5  array.  Since  electromagnetic  showers  tend  to  be  more  narrowly  contained 
than  hadronic  showers,  they  mostly  have  E9/E25  near  1.  The  most  widely  used 
cut  on  E9/E25  quantity  in  CLEO  is  to  keep  99%  photons*  and  get  rid  of  most  of 
the  hadronic  showers.  A  E9/E25  cut  also  gets  rid  of  two  photons  which  are  close 
together.  Also,  as  the  high  energy  showers  near  to  each  other  often  deposit  energy 
to  the  same  crystals,  we  can  use  a  program  to  "unfold"  the  nearby  shower  clusters. 
The  unfolding  routine  works  well  to  identify  high  energy  photons  but  become  less 
efficient  when  photon  energy  is  low.  When  we  select  a  photon  in  this  analysis,  we 
mostly  use  unfolded  E9/E25  at  99%  (this  cut  is  called  99%  E9E25U),  and  sometimes 
directly  use  E9/E25  at  99%  (not  unfolded,  called  99%  E90E25)  when  we  are  certain 
that  the  photon  energy  spectrum  is  low. 

7r°'s  are  formed  by  taking  two  photon  combinations  and  calculating  the  di-photon 
invariant  mass.  The  di-photon  combination  must  have  at  least  one  photon  in  the 
barrel  region  of  the  crystal  calorimeters. 

Unless  otherwise  stated,  we  always  selected  7r°  candidates  as  just  described.  In 

addition,  we  require  the  mass  of  each  candidate  to  be  within  2.5  oM  0  of  the  known 

7T°  mass,  which  is  134.976  MeV/c2.  However,  to  suppress  combinational  background, 

we  will  place  cuts  on  7r°  momenta  in  different  situations. 

"This  cut  value  is  actually  obtained  by  a  program  which  returns  E9/E25  cuts  for  1%  efficiency 
loss  in  a  given  region  of  the  detector.  The  1%  is  nominal  only  from  MC,  and  it  is  a  loss  in  an  event 
with  no  other  showers.  Actual  efficiency  losses  are  substantially  more,  e.g.,  5-10%  efficiency  loss 
for  photons  in  hadronic  events  from  a  1%  cut.  Harsh  cuts,  like  5%  or  10%,  lose  too  many  photon 
candidates,  so  they  are  not  often  used. 
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3.3.4  Event  Vertex  Finding 

Finding  the  event  vertex,  which  is  the  e+e~  interaction  point,  can  be  very  useful 
to  improve  the  mass  resolution  of  a  reconstructed  particle.  Most  of  the  charged 
tracks  come  from  the  event  vertex  or  vertices  extremely  close  to  the  event  vertex,  so 
we  can  obtain  the  event  vertex  with  errors  by  performing  a  vertex  fit  using  charged 
tracks.  However,  some  tracks  are  not  from  anywhere  close  to  the  event  vertex, 
and  some  tracks  are  not  measured  correctly.  Using  these  tracks  will  increase  the 
uncertainty  of  the  event  vertex.  Therefore,  to  find  the  event  vertex,  we  start  with 
the  two  tracks  which  are  the  most  consistent  with  the  beam  position  and  also  not 
close  to  back  to  back,  form  a  vertex,  then  add  more  tracks  one  by  one  in  the  order 
of  their  consistencies  with  the  current  found  vertex,  until  further  adding  any  track 
would  make  a  poorer  vertex  fit.  If  the  error  of  the  final  found  vertex  is  not  small 
enough  to  confine  itself  within  the  spread  of  the  beam  position,  or  if  we  do  not  have 
more  than  one  charged  track,  the  beam  position  will  be  used  as  the  event  vertex. 

3.3.5  Vee  Finding 

A  large  fraction  of  A+'s  decay  to  tf°'s  or  A's.  In  CLEO,  the  K°  is  reconstructed 
by  the  decay  Ks  ->•  7r+7i-*  and  the  A  is  reconstructed  by  its  decay  to  jm~ .  Ks's 
and  A's  mostly  fly  a  few  centimeters  away  before  they  decay.  The  Ks  has  cr  =  2.7 
cm  and  the  A  has  cr  =  7.9  cm.  Since  in  each  case  the  shape  of  the  two  daughters' 
paths  looks  like  the  letter  V,  they  are  generally  called  Vees.  The  vertex  of  a  Vee 
formed  by  the  two  charged  tracks  is  called  a  secondary  vertex,  in  distinction  to  the 
primary  vertex,  the  event  vertex. 

To  find  Vee  candidates,  we  first  list  the  tracks  which  satisfy  — DBCD/DDBCD  >  3— 

(DDBCD  is  the  estimated  uncertainty  of  DBCD)  to  exclude  the  tracks  coming  from  the 

*K°  appears  in  the  states  of  Ks  (50%)  or  KL  (50%).  CLEO  detector  cannot  detect  the  KL 
since  it  travels  many  meters  long  (cr  =  15.5  m). 


primary  vertex.  The  tracks  in  the  list,  called  the  secondary  tracks,  are  all  possible 
daughters  of  Vees.  Then  we  take  every  pair  of  tracks  from  the  list  and  fit  them  to 
form  a  vertex.  We  discard  the  pair  if  the  fit  is  very  poor  (x2  <  0.001)  or  if  the 
significance  of  the  three-dimensional  flight  distance  is  smaller  than  3a.  We  require 
the  total  number  of  hits  inside  the  found  vertex  for  both  daughters  to  be  smaller 
than  2.  Also,  the  reconstructed  mass  of  a  Vee  should  be  within  3  aM  of  the  known 
Ks  or  A  mass.  Note  that  the  Kalman  Filter  only  calculates  the  track  parameters 
of  the  secondary  tracks  from  the  origin.  We  corrected  these  track  parameters  at 
the  secondary  vertex,  then  further  recalculated  the  momentum  and  the  position  of 
the  Vees.  Details  of  the  CLEO  secondary  vertex  finding  program  are  described  by 
Prescott[24]. 

In  this  analysis,  we  made  a  few  more  cuts  on  Vee  candidates  to  make  the  signal 
clearer  with  a  tiny  sacrifice  of  losing  events:  a  Vee  is  required  to  come  from  the 
range  of  the  beam  position;  the  separation  in  z  of  the  two  daughters  at  the  r  -  <f> 
intersection  should  not  be  large  (the  significance  is  required  to  be  smaller  than  5); 
for  a  Ks  candidate,  each  daughter  must  satisfy  PROB(tt)  >  0.001;  for  a  A  candidate, 
the  daughter  which  carries  most  of  the  momentum  in  the  decay  in  the  laboratory 
frame  (according  to  special  relativity,  this  almost  guarantees  that  it  is  the  proton  if 
its  parent  is  really  a  A)  must  satisfy  PR0B(p)  >  0.001. 

3.3.6    £  and  H  Finding 

The  decays  of  the  four  £  and  3  are  shown  in  Figure  3.1.  Because  these  particles 
are  supposed  to  be  the  decay  products  from  the  A+,  which  has  very  short  lifetime, 
it  is  very  reasonable  to  assume  that  they  all  come  from  the  event  vertex. 

£+  candidates  are  reconstructed  through  the  decay  £+  ->  p7r°.  The  tt0  selection 
is  described  in  Section  3.3.3.  All  charged  tracks  are  taken  as  proton  candidates. 
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Figure  3.1:  Decay  paths  of  (a)  E+  ->  p7r°  and  tt°  77,  (b)  E°  -»  A7  and  A  ->.  p^-, 
(c)E-  ->  Att"  and  A  -»  p*--,  and  (d)  S°  -»  A7r°,  A  ->  p7r~,  and  tt0  ->  77.  In  each 
figure,  the  solid  curves  indicate  paths  of  charged  particles;  the  dashed  lines  indicate 
paths  of  photons;  and  the  dotted  lines  indicate  paths  of  neutral  hyperons.  All  £ 
and  S  are  assumed  to  be  created  at  the  event  vertex,  which  is  shown  as  •.  The  £+, 
-~»  s°>  and  A  travel  for  a  distance  before  decay.  The  E°  and  7r°  decay  immediately 
after  being  created. 
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A  particle  identification  cut  PRLEV  >  0.5  is  used.  After  adding  the  4-momenta  of 
the  proton  and  the  7r°,  the  flight  distance  of  E+  can  be  calculated  by  constraining 
the  E+  coming  from  the  event  vertex,  and,  further,  the  E+  decay  vertex  can  be 
obtained.  Since  the  E+  has  a  long  lifetime  (cr  =  2.4  cm),  we  require  the  transverse 
flight  distance  of  the  E+  to  be  bigger  than  5  mm  and  DBCD  of  the  proton  to  be 
bigger  than  0.5  mm.  In  addition,  to  reduce  the  combinational  background,  we  make 
a  momentum  cut  of  300  MeV/c2  for  the  tt°  and  500  MeV/c2  for  the  proton.  We 
require  the  mass  of  each  E+  candidate  to  be  20  MeV/c2  within  the  known  E+  mass, 
1189.4  MeV/c2. 

E°  candidates  are  reconstructed  through  the  decay  E°  ->  A7,  where  selections 
of  the  A  and  7  are  described  in  the  previous  sections.  A  50  MeV/c2  photon  energy 
cut  is  used  to  reduced  the  background.  We  require  the  mass  of  each  E°  candidate 
to  be  20  MeV/c2  within  the  known  E°  mass,  1192.6  MeV/c2. 

E~  candidates  are  reconstructed  through  the  decay  H~  Air~.  All  negative- 
charge  tracks  are  used  to  be  the  n~  candidates.  The  E~  has  cr  =  4.9  cm  so  its 
decay  vertex  should  be  distinctly  displaced  from  the  event  vertex.  The  initial  guess 
of  the  H~  vertex  is  at  the  intersection  of  the  projected  A  path  and  the  jt-  track.  The 
combination  is  discarded  if  the  vertex  cannot  be  reasonably  formed  or  the  combined 
direction  does  not  point  back  to  the  event  vertex.  The  invariant  mass  of  Ayr  is 
required  to  be  within  6  MeV/c2  of  the  known  H"  mass,  which  is  1197.4  MeV/c2. 
Then  the  track  parameters  of  the  A  and  the  n~  are  recalculated  at  the  new  vertex 
by  constraining  the  An~  invariant  mass  to  be  1197.4  MeV/c2.  The  H~  momentum 
and  vertex  are  further  obtained  more  accurately.  The  E~  flight  distance  is  required 
to  be  larger  than  5  mm,  and  its  r  -  0  radial  distance  is  required  to  be  larger  than 
1  mm. 

H°  candidates  are  reconstructed  by  a  E°  -»  Att°.  We  initially  build  the  H°  vertex, 
then  combine  the  tt°  angle  with  the  projected  A  path.  The  combined  direction  is 
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required  to  point  back  to  the  event  vertex.  We  correct  the  77°  4-momentum  by  con- 
straining it  to  come  from  the  E°  vertex.  Then  the  H°  4-momentum  is  recalculated. 
The  reconstructed  ~°  mass  is  required  to  be  within  8  MeV/c2  of  the  known  H°  mass, 
which  is  1315  MeV/c2.  The  r  -  <\>  radial  distance  of  the  H°  vertex  from  the  event 
vertex  has  to  be  larger  than  5  mm.  A  100  MeV/c2  ir°  momentum  cut  is  used. 

3.4    A+  Selection 

A  A+  is  reconstructed  by  calculating  its  4-momentum  and  decay  vertex  from 
the  momentum,  mass  assumption,  and  vertex  information  of  its  direct  daughters. 
The  A+  signal  can  be  measured  by  fitting  its  reconstructed  mass  spectrum  us- 
ing Mn_Fit.[25]*  Mn_Fit  is  an  interactive  plotting  and  fitting  package  that  uses 
MINUIT[26]  to  fit  histograms  or  data  read  in  from  a  file  and  to  display  the  fit- 
ting results.  It  is  widely  used  by  high  energy  experimental  groups  to  perform  data 
analysis. 

Although  the  purpose  of  reconstructing  the  A+  baryons  is  to  observe  the  A+x 
and  the  Ec  which  decay  to  the  A+,  our  A+  candidates  are  not  selected  by  direct 
optimizing  the  A+  or  E«  signals.  The  A+  or  Ec  signals  should  not  be  biased  by 
A+  selection.  In  general,  for  each  A+  decay  mode  we  try  to  use  cuts  to  get  the  best 
signal-to-background  ratio  at  xp  >  0.5.  xp  is  the  A+  momentum  normalized  to  the 
maximum  possible  A+  momentum,  defined  as 


Xp  — 


where  E  is  the  energy  of  the  e+(e~)  beam  of  CESR.  Note  that  xp  >  0.5  is  only 
used  to  optimize  the  A+  selection  cuts,  and  not  for  A+  selection.  The  momentum 

*In  this  dissertation,  all  histograms,  fitted  curves  and  direct  statistical  results  are  obtained  from 
Mn-Fit. 
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spectra  of  charmed  hadrons  from  the  e+e~~  annihilation  have  been  shown  to  be  hard 
compared  to  the  combinational  background.  Therefore  a  cut  on  xv  is  a  very  good 
way  to  improve  the  signal-to-background  ratio.  For  most  of  the  decay  modes,  the 
xp  >  0.5  cut  makes  the  signal  be  roughly  same  large  as  the  background.  Later, 
when  we  look  at  or  signals,  we  make  xp  cuts  on  these  heavier  baryons,  and 
there  will  be  no  cut  on  the  xp  of  the  A+. 

Since  a  A+  travels  at  a  distance  of  the  order  of  100  /im  before  it  decays,  its  vertex 
can  be  slightly  displaced  from  the  event  vertex.  We  find  that  constraining  all  direct 
A+  daughters  (except  7r°'s  which  do  not  have  position  information)  to  come  from 
an  unknown  vertex  significantly  improves  the  mass  resolution  of  the  A+.  Thus  for 
each  of  the  15  decay  modes,  the  A+  candidate  is  built  at  its  vertex. 

For  the  decay  modes  containing  a  proton,  strong  particle  identification  cuts  are 
needed.    Since  the  proton  momentum  spectrum  are  different  for  different  decay 
modes  (generally  higher  for  two-body  decays  and  lower  for  three-body  and  four- 
body  decays),  the  cut  value  on  PRLEV  does  not  correspond  to  the  same  efficiency  for 
the  different  modes  of  background  rejection.  The  numbers  we  used  —  0.3  for  pK°, 
0.6  for  pK-7r+,  and  0.8  for  pK-yr+Tr0,  pKVtt0,  and  pKV—  have  roughly  the  same 
efficiencies.  For  the  two  modes  containing  a  proton  and  a  K~,  the  identification  cuts 
are  relatively  softer  since  the  proton  is  already  strongly  identified.  Since  the  signals 
of  E+K+K",  E-K+7T+,  and  ~°K+are  already  very  clean,  very  soft  PWB(K)  cuts  are 
used  to  keep  the  highest  efficiencies.  To  suppress  combinational  background,  Monte 
Carlo  optimized  momentum  cuts  on  charged  tracks  or  7r°'s  are  used.  For  charged 
particles,  we  generally  cut  at  high  values  for  two-body  decays  and  lower  values  for 
multi-body  decays.  All  decay  modes  containing  a  tt0  have  very  high  backgrounds, 
so  we  always  need  restrictive  cuts  on  7r°  momenta. 
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Table  3.2:  Particle  Identification  and  Momentum  cuts  used  in  the  15  decay  modes 
for  the  A+  selection.  The  values  of  physical  parameters  in  the  table  are  required  to 
be  larger  than  the  corresponding  numbers  listed.  The  momenta  PK±,  Pn±,  and  Pno 
are  in  units  of  MeV/c. 


Decay 

PRLEV 

KALEV 

PROB(K) 

1  /i  u 

pK-7r+ 

0.6 

0.2 

100 

pK-TT+TT0 

0.8 

0.5 

200 

500 

pK° 

0  3 

piv  7T  7T 

U.o 

500 

pKV 

0.8 

500 

Att+ 

300 

A7r+7r+7r~ 

200 

Att+tt0 

300 

500 

AK+K° 

0.3 

E+7T+7r_ 

100 

E+7T0 

300 

E+K+K" 

0.02 

100 

E°7r+ 

400 

E-K+7T+ 

0.01 

100 

0.01 

300 

The  cuts  used  to  reconstruct  A+'s  for  each  decay  mode  are  listed  in  Table  3.2. 
Note  that  the  vertex  constraint  routines  are  common  for  all  modes.  The  selection 
criteria  of  all  immediate  A+  daughters  have  been  described  in  previous  sections. 

Figures  3.2  to  3.16  show  the  reconstructed  mass  of  A+  for  the  15  idividual  decay 
modes,  and  Figure  3.17  shows  the  sum  over  all  modes.  For  the  purpose  of  display  a 
xp  >  0.5  cut  is  applied  to  the  A+  momentum. 

The  measurements  on  all  15  modes  and  their  summation  are  tabulated  in  Table 
3.3.  All  parameters  are  obtained  by  fitting  the  signal  with  a  single  Gaussian.  The 
Gaussian  width  of  each  decay  mode  reflects  the  detector  resolution  of  the  A+  mass. 
The  table  shows  that  the  fitted  widths  agree  very  well  with  those  of  the  Monte  Carlo, 
except  some  modes  with  a  tt°  and  "All."  For  the  modes  with  a  tt°,  we  found  that 
the  widths  are  all  wide  and  the  signals  could  not  be  fitted  well  by  a  single  Gaussian. 
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Figure  3.2:  Invariant  mass  of  pK~7r" 
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Figure  3.3:  Invariant  mass  of  pK  n^Tr0. 
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Figure  3.4:  Invariant  mass  ofpK°. 
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Figure  3.5:  Invariant  mass  of  pK°7r07r°. 
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Figure  3.6:  Invariant  mass  of  pK°7r°. 
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Figure  3.7:  Invariant  mass  of  An+. 
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Figure  3.9:  Invariant  mass  of  A7r+7r°. 
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Figure  3.11:  Invariant  mass  of  E+7r+7r~. 
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Figure  3.14:  Invariant  mass  of  E°7r+. 
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Figure  3.15:  Invariant  mass  of  H~K+7r+. 
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Figure  3.16:  Invariant  mass  of  E°K+. 
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Figure  3.17:  Invariant  mass  of  A+  candidates  reconstructed  by  15  decay  modes. 
Although  the  distribution  overlaid  by  15  different  Gaussians  is  not  exactly  in  a 
Gaussian  shape,  we  fit  it  by  a  single  Gaussian  to  roughly  check  the  central  mass 
value  and  the  signal-to-background  ratio. 
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Table  3.3:  Masses,  Monte  Carlo  widths  (aM(MC)),  Fitted  widths  (aM(Fit))  together 
with  the  number  of  A+'s  found  with  xp  >  0.5.  The  mass  for  all  modes  is  the  error- 
weighted  mean  value  over  all  modes  and  the  same  value  is  obtained  when  we  fit 
Figure  3.17.  The  signal  the  background  shown  in  the  table  are  the  fitted  areas  within 
1.6crM(Fit)  around  2286.7  MeV/c2.  Masses  and  widths  are  in  units  of  MeV/c2.  Only 
statistical  errors  are  shown.  The  statistical  errors  on  Monte  Carlo  width  are  all  0  01 
MeV/c2. 


Decay  Mode 

Mass 

aM(MC) 

a  M  (Fit) 

Signal 

Background 

pK-7r+ 

2286.7±0.1 

5.20 

5.26±0.11 

9650 

11660 

pK-7r+7r0 

pK  TY+TC~ 

pkV 

2286.4±0.6 

8.11 

9.04±0.62 

1252 

1928 

2287.1±0.2 

6.51 

6.33±0.23 

1716 

1019 

2287.1±0.3 

3.98 

4.07±0.30 

568 

690 

2285.2±0.6 

8.61 

7.63±0.63 

514 

490 

Att+ 

2287.8±0.3 

7.03 

7.15±0.31 

1368 

658 

A7r+7r+7r_ 

2287.0±0.2 

5.28 

5.37±0.21 

1916 

1812 

Att+tt0 
AK+K° 

2285.0±0.4 

10.38 

11.18±0.56 

2007 

1864 

2286.5±0.4 

2.91 

2.51±0.32 

70 

25 

E+7T+7T~ 

2287.6±0.3 

7.41 

6.96±0.38 

947 

696 

E+7T0 

2284.0±1.7 

10.32 

9.80±0.23 

95 

101 

E+K+K- 

2287.0±0.4 

2.61 

2.54±0.31 

65 

18 

E°7T+ 

2286.5±0.6 

7.54 

7.79±0.59 

707 

896 

H-K+7T+ 

2286.7±0.1 

3.50 

3.27±0.25 

317 

171 

H°K+ 

2286.9±0.5 

6.73 

6.10±0.51 

221 

173 

All 

2286.7±0.1 

20675 

22160 

If  we  fit  them  by  a  double-Gaussian,  the  data  and  the  Monte  Carlo  are  in  better 
agreement  (about  half  way  closer).  Thus  we  can  conclude  that  the  Monte  Carlo  can 
predict  the  mass  resolution  of  the  A+  very  well. 

Note  that  our  measured  mean  value  of  A+  mass,  2286.7  MeV/c2,  is  slightly 
larger  than  the  PDG  mean  value  2284.9±0.6  MeV/c2,  which  is  partially  determined 
by  the  result  of  CLEO-I.5[27].  CLEO-II  has  measured  many  more  decay  modes 
with  much  better  mass  resolutions  and  tiny  statistical  errors,  and  has  constantly 
shown  a  higher  value.  So  in  this  work,  we  decided  to  use  the  CLEO-II  value  as 
the  "known"  A+  mass.  As  we  will  show  in  later  chapters,  as  we  look  at  A+,  or  £(*) 
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mass  by  mass  differences  with  A+,  the  error  on  A+  mass  contributes  little  to  the 
systematic  errors  of  our  measurements.  The  purpose  of  this  work  does  not  include 
the  A+  mass  measurement.  It  should  be  a  separate  research  work  which  will  require 
a  comprehensive  study  of  the  systematic  uncertainties  of  the  mass  scale. 


CHAPTER  4 
STUDIES  OF  A+  BARYONS 


4.1  Introduction 

Two  particles  have  been  found  decaying  to  A+  and  two  charged  pions.  They 
have  been  named  A+(2625)  and  A+(2593)  by  the  Particle  Data  Group[l].  The 
first  observation  of  the  A+(2625)  was  made  by  the  ARGUS  Collaboration[28].  The 
E687[29]  and  CLEO[30]  Collaborations  subsequently  confirmed  this  charmed  state. 
CLEO  also  reported  the  first  observation  of  A+ (2593),  and  E687[31]  and  ARGUS[32] 
published  the  confirmation  of  these  states.  Since  hese  particles  were  seen  only  in 
A+7T+7T-  modes  but  not  in  A+tt°  mode,  they  are  believed  to  be  the  excited  A+ 
baryons  rather  than  higher  states  of  Ec  baryons,  which  decay  through  both  A+?r 
and  A+7T+7T-  channels  but  should  prefer  A+?r  modes.  Proposed  by  the  previous 
experiments  and  strongly  confirmed  in  this  work,  the  two  states  are  a  fine  structure 
doublet  of  the  excited  A+  baryons  in  which  the  light  diquark  ud  has  an  orbital 
angular  momentum  L  =  1  with  respect  to  the  c  quark,  Therefore,  even  though  they 
were  called  A*+  in  previous  experimental  publications  just  referred  above,  we  prefer 
to  call  them  the  A+(2593)  and  A+(2625),*  where  subscript  1  refers  to  the  total  spin 
of  the  light  degrees  of  freedom. 

The  A+j(2593)  decays  largely,  and  perhaps  entirely,  into  £ctt.  As  explained  in 
Chapter  1,  this  is  the  expected  dominant  decay  channel.  The  quark  model  predicts 
that  Ec  has  Jp  =  1/2+  so  almost  certainly  the  A+^2593)  has  Jp  =  1/2".  The 

*They  are  called  A+(l/2)  and  A+ (3/2)  in  Chapter  1  in  order  to  conviently  distinguish  the  spin 
difference  But  the  general  convention  to  name  a  excited  state  particle  is  to  used  the  particle  name 
followed  by  its  first  observed  mass. 
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Jp  of  the  A+j(2625)  is  expected  to  be  3/2".  The  A+i(2593)  and  the  A+^2625)  are 
presumably  the  charm  counterparts  of  the  strange  particles  the  A(1405)  and  the 
A(1520). 

Since  the  last  publication  of  these  states  in  1995,  the  CLEO-II  detector  has 
collected  60%  more  data,  and  the  recompress  made  in  1996  significantly  improved 
the  resolution  of  low-momentum  tracks,  which  is  very  important  in  order  to  observe 
A*i->A+7r+7r-  decays.  We  think  it  is  time  to  confirm  or  update  our  previous  results 
and  look  for  a  new  decay  channel,  A+7r°7r°. 

In  this  chapter,  we  first  describe  Monte  Carlo  simulation  of  A+x  decays.  This  is 
followed  by  the  observations  of  A+tt+tT  and  A+ttV  signals.  Then  we  present  our 
measurements  on  Ec7r  decay  substructures  and  the  search  for  the  decays  of  A+7r° 
and  A+7.  In  later  sections,  we  report  the  results  on  masses,  widths,  fragmentation 
functions,  branching  ratios,  production  ratios,  together  with  a  detailed  discussion 
on  estimating  systematic  errors.  We  also  compare  our  results  with  theoretical  pre- 
dictions. 

4.2    Monte  Carlo  Studies 

In  order  to  determine  the  detector  resolution  of  signals  and  estimate  the  recon- 
struction efficiencies,  large  amounts  (more  than  100,000  events  for  each  decay  mode) 
of  Monte  Carlo  events  were  generated.  Therefore,  statistical  errors  on  all  quanti- 
ties directly  derived  from  Monte  Carlo  can  always  be  neglected  in  the  numerical 
calculations. 

The  A+x(2593)  events  are  first  generated  with  a  mass  of  2593.9  MeV/c2  and 
zero  width  to  obtain  the  mass  resolution  of  the  detector.  After  a  preliminary  mea- 
surement of  the  natural  width  using  the  real  data,  more  A+i(2593)  events  with  the 
measured  natural  widths  were  generated.  We  found  that  our  invariant  mass  recon- 
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struction  program  can  extract  the  Monte  Carlo  generated  natural  width  very  well, 
and  efficiencies  calculated  from  0  width  events  and  events  with  a  finite  width  are 
very  close. 

The  A+1(2625)  events  are  generated  with  a  mass  of  2626.6  MeV/c2  and  zero 
width. 

To  study  the  substructure  of  decays,  we  generated  samples  with  both  Ec7r 
and  non-resonance  decay  modes.  The  non-resonance  modes  are  simply  the  decays 
A+^A+tttt  and  A+^A+ttV.  The  Ectt  modes  are  (1)  the  mixed  50%  A+^E^+tt- 
and  50%  A+->E°7r+,  with  E++->A+7r+  and  E^A+tt",  and  (2)  the  A+-+  E+tt0  with 
E+->A+7T°.  TheA+!-^Ec7r  and  A+^A+tttt  modes  are  used  for  the  A;!i(2593)  and  the 
A+x(2625),  respectively,  to  evaluate  the  corresponding  efficiencies  and  resolutions,  as 
the  data  show  us  that  these  are  the  dominating  decay  modes. 

Determination  of  the  detector  resolution  by  the  Monte  Carlo  is  very  crucial  for 
extracting  the  intrinsic  width  of  a  particle  when  its  intrinsic  width  is  close  or  smaller 
than  the  resolution.  As  stated  in  the  last  chapter,  we  have  checked  for  all  16  A+ 
decay  modes,  and  the  Monte  Carlo  predicted  widths  are  all  in  good  agreement  with 
the  fitted  widths  of  the  data  signal.  However,  the  charged  daughters  of  A+'s  are 
all  at  quite  high  momentum  ranges,  but  the  two  independent  pions  in  A+x->  A+7T7T 
decays  are  both  soft.  It  is  possible  that  the  CLEO-II  official  Monte  Carlo  program 
overestimates  the  performance  of  the  detectors  on  measurements  of  low-momentum- 
track  parameters.  This  can  make  the  Monte  Carlo  predicted  resolutions  smaller  than 
actual  resolutions  of  A^-^  A+nn  signals. 

Measuring  the  signal  of  the  D*+,  which  has  F  <  0.13  MeV/c2  at  90%  confidence 
level[l],  can  be  used  to  check  the  accuracy  of  the  resolution  of  the  low-momentum 
pions  by  the  Monte  Carlo.  To  do  this,  we  looked  at  the  mass  difference  M(D*+)  - 
M(D°)  of  decay  D*+  ->  D°tt+  (and  D°  ->  K~ir+)  in  which  the  momentum  spectrum 
of  the  independent  pion  is  close  to  those  of  the  two  independent  pions  from  A+, 
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decays.  The  fitted  width  of  the  signal  of  M(D*+)  -  M(D°)  distribution  is  (0.550 
±  0.001)  MeV/c2  for  the  Monte  Carlo,  but  about  0.63  MeV/c2  for  the  real  data 
in  the  worst  case,*  about  15%  larger.  We  think  the  actual  detector  resolution  is  in 
the  0  -  15%  range  larger  than  the  Monte  Carlo  predicted  one.  The  D*+  width  is 
largely  affected  by  the  angular  resolution  of  the  soft  pion.  This  is  not  the  case  for 
widths  based  on  our  study  of  comparing  the  Monte  Carlo  widths  with  the  data 
at  different  decay  angle.  So  the  amount  of  the  disagreement  between  real  data  and 
the  Monte  Carlo  for  D*+  may  not  reflect  the  same  amount  of  the  disagreement  for 
Ati  widths.  But  it  appears  that  the  Monte  Carlo  prediction  of  the  resolution  is  too 
optimistic.  This  was  not  noticeable  in  the  previous  research  in  CLEO  since  the  mass 
resolution  itself  was  poor.  For  the  CLEO-II  recompressed  data,  this  effect  should  be 
seriously  considered  as  a  big  source  of  systematic  uncertainties  for  particles  whose 
natural  widths  and  mass  resolutions  are  small. 

Some  studies  are  performed  to  check  the  correctness  of  Monte  Carlo  determined 
efficiencies  of  particle  reconstruction.  We  calculated  the  decay  ratios  of  each  A+ 
decay  mode  relative  to  PK-tt+  based  on  the  actual  yields  from  the  data  and  Monte 
Carlo  predicted  efficiencies,  and  they  all  agree  within  10%  with  the  Monte  Carlo 
generated  ratios.  As  we  measure  all  production  and  branching  fractions  by  ratio 
relative  to  A+,  the  effect  of  the  error  of  A+  reconstruction  efficiency  mostly  cancels 
in  the  ratio. 


For  the  real  data,  the  M(D*+)  -  M(D°)  spectrum  is  difficult  to  fit  due  to  the  uncertain 
background  shape  near  the  peak,  which  is  very  close  to  the  kinematic  threshold.  When  we  fitted 
the  signal  at  a  different  pion  momentum  range,  we  found  that  fitted  width  is  closer  to  the  Monte 
Carlo  predicted  resolution  at  higher  momentum  pions.  By  using  different  fitting  procedures  to  fit 
signals  at  different  pion  momentum  ranges,  0.63  MeV/c2  is  the  largest  fitted  width.  We  believe 
that  15%  is  a  conservative  estimate  of  how  much  the  Monte  Carlo  miscalculates  the  mass-difference 
resolution. 
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4.3    Signals  and  Fits 

To  observe  signals,  we  combine  each  A+  candidate  with  two  pions  (ir+n~ 
or  7r°7r°),  reconstruct  the  invariant  mass  with  the  A+,  then  fit  the  spectrum  of  the 
mass  difference  M(A+7T7t)  -  M(A+).  A+  selection  criteria  are  described  in  the  last 
chapter,  and  the  reconstructed  masses  of  A+  candidates  are  required  to  be  within 
l.6aM  of  2.2867  MeV/c2,  our  fitted  mean  value  of  the  A+mass.  As  A+j  decays 
are  strong  decays,  their  lifetimes  are  very  short.  We  can  assume  that  they  decay 
at  the  interaction  point  of  e+e~  annihilation.  The  event  vertex  finding  procedure 
is  described  in  Chapter  3.  When  we  reconstruct  the  invariant  mass  of  each  A+7T7t 
combination,  we  constrain  the  three  particles  to  be  from  the  event  vertex.  We  find 
that  using  vertex  constraint  improves  the  resolution  of  the  A+1(2625)  by  more  than 
10%  for  the  A+7r+7i-case.  In  order  to  see  very  clean  signals,*  we  require  xp  >  0.7 
for  each  A+irn  combination.  xp  is  defined  as 


yfELm  -  Ml+  ' 


where  pA+  is  the  total  momentum  of  the  A+tttt  combination,  Ebeam  is  the  beam 
energy  of  the  e+(e~),  and  MA+  is  the  mass  of  the  A£. 


4.3.1  A+7T+7T- 


c 


Shown  in  Figure  4.1  is  the  mass-difference  spectrum  of  A+tt+tt"  combinations. 
The  large  peak  at  342  MeV/c2  is  due  to  the  decay  Ac+!(2625)  ^A+tt+tt",  and  the 
smaller  peak  at  approximately  309  MeV/c2  is  from  A+1(2593)^A+7r+7r-. 

To  fit  the  the  mass-difference  distribution,  we  used  a  second-order  Chebyshev 
polynomial  to  fit  the  background  and  two  signal  functions  to  fit  the  peaks.  Each 

*As  we  can  see  in  later  sections,  having  low  background  is  very  important  to  measure  the  E.tt 
substructure. 
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Figure  4.1:  Mass  difference  M(A+7r+7r-)  -  Af(A+) 
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signal  function  is  a  Breit-Wigner  function  with  a  floating  width  convoluted  with  a 
single  Gaussian  function  with  a  fixed  width.  The  Breit-Wigner  widths  reflect  the 
intrinsic  widths  of  the  particles;  while  the  Gaussian  widths,  which  are  1.28  MeV/c2 
for  the  A+(2593)  and  1.62  MeV/c2  for  the  A+^2625)  determined  from  the  Monte 
Carlo  studies,  are  used  to  parameterize  the  detector  resolutions. 

In  order  to  prove  that  signals  in  Figure  4.1  are  not  artifacts  of  special  selection 
criteria,  we  checked  the  backgrounds  from  two  sources:  one  is  the  combination  of  a 
fake  A+  candidate  taken  A+  sidebands  and  7r+7r~,  the  other  is  the  combination  of 
a  A+  with  two  same-sign  pions  (7r±7r±).  As  Figure  4.2  shows,  these  mass-difference 
distributions  are  smooth  and  similar  to  the  A+7T+7T-  background.  These  plots  can 
be  fitted  very  well  using  the  fixed  background  parameterization  as  used  in  Figure 
4.1. 

As  shown  in  Figure  4.1,  it  is  obvious  that  the  A+^2593)  peak  is  not  fitted  well 
by  a  Breit-Wigner  convoluted  with  a  Gaussian.  Fitting  uncertainties  on  both  the 
mass  and  the  width  are  big.  We  will  discuss  this  further  in  the  next  chapter  after 
presenting  the  Sc  mass  measurements.* 

The  extracted  Breit-Wigner  width  of  the  A+j(2625)  with  a  1.62  MeV/c2  Gaussian 
resolution  is  0.7±0.3  MeV/c2,  close  to  zero.  But  the  question  is,  have  we  measured 
a  finite  width  of  the  A+^2625)? 

As  we  mentioned  in  the  last  section,  for  CLEO-II  the  actual  detector  resolution 

may  be  as  much  as  15%  larger  than  the  Monte  Carlo  predicted  value.  If  we  fit  the 

A^i(2625)  with  a  single  Gaussian  with  a  floating  width,  the  fitted  width  is  1.76±0.12 

MeV/c2,  less  than  10%  larger  than  the  Monte  Carlo  predicted  value.  So  taking 

statistical  errors  into  account,  we  cannot  conclude  that  we  have  observed  a  finite 

natural  width  of  the  A+(2625).  Therefore,  in  the  following  study,  we  calculated  the 

*We  also  remeasured  the  A+(2593)  signal  by  a  new  technique.  This  second  method  is  only 
possible  after  a  measurement  of  Ec  width  and  mass.  The  measurements  give  different  mass  and 
width  values. 
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Figure  4.2:  The  top  plot  shows  the  mass  difference  M{K+tt~ti-)  -  M(A+)  and  the 
middle  plot  shows  M (A+7r+7r+)-M(A+),  and  Ac+'s  candidates  are  taken  with  masses 
with  1.6  aM^  from  2.2867  MeV/c2.  The  bottom  plot  shows  M(A+7r+7r-)  -  M("A+ 
side"),  where  a  "Ac+  side"  is  from  A+  sidebands  which  are  taken  from  3  4  to  5  0 
(Jms+  away  from  2.2867  MeV/c2. 
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production  rate  of  the  A+^2625)  by  using  the  yield  of  the  signal  fitted  by  a  single 
Gaussian,  which  is  roughly  10%  smaller  than  the  yield  obtained  from  the  fit  using  a 
Breit-Wigner  convoluted  with  a  Gaussian.  Corresponding  systematic  uncertainties 
were  added  in  the  final  results  of  branching  ratios  and  production  ratios. 

Table  4.1  shows  the  yields,  mass  differences,  and  widths  obtained  by  using  dif- 
ferent fits  to  Figure  4.1.  Fit  Type  1  uses  a  Breit-Wigner  convoluted  with  a  Gaussian 
for  the  A+^2593)  peak  and  a  single  Gaussian  for  the  A+j(2625)  peak.  Type  2  fits 
the  region  right  to  294  MeV/c2,  excluding  two  entries  from  the  plot.  Type  3  uses 
a  third-order  polynomial  to  fit  the  background  instead  of  the  second-order  one. 
Higher-order  polynomials  were  also  used,  and  the  results  are  very  similar.  Types  4, 
5,  and  6  only  fit  the  A+^2593)  peak.  Type  4  excludes  the  region  from  325  MeV/c2  to 
360  MeV/c2;  Type  5  excludes  from  335  MeV/c2  to  350  MeV/c2;  Type  6  only  fits  the 
region  from  290  MeV/c2  to  325  MeV/c2.  Types  7  and  8  only  fit  the  A+(2625)  peak. 
Type  7  fits  only  from  335  MeV/c2  to  390  MeV/c2  and  Type  8  excludes  from  300 
MeV/c2  to  325  MeV/c2.  Type  9  uses  a  Breit-Wigner  convoluted  with  a  Gaussian  to 
fit  the  A+ (2625)  peak  instead  of  a  single  Gaussian.  The  errors  from  Type  2  to  Type 
8  are  not  shown  since  they  are  roughly  the  same  as  in  Type  1.  We  have  performed 
more  variations  of  fits  other  than  listed  in  Table  4.1.  Variations  of  the  results  are 
all  within  the  range  Table  4.1  shows. 

Different  values  of  xp  and  cut  around  the  A+  mass  may  vary  the  results  of  the 
mass  and  the  width.  As  we  checked,  the  variation  of  the  widths  on  both  particles 
and  the  mass  difference  of  the  A+(2625)  is  less  than  the  one  shown  in  Table  4.1,  but 
the  mass  of  the  A+(2593)  varies  from  308.2  MeV/c2to  309.1  MeV/c2. 
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Table  4.1:  The  measured  yield,  mass,  and  extracted  width  of  A^-^A+tf+jr-  by 
various  fits. 


Fit  Tvnp 
ail  J-  v ^/c 

YiplH 

I  lclU 

iviass^ivie  v  c  I 

Widtn(MeV/c  ) 

A+  ( 25931 

1 

1 

i  qi  +31 
iyl-28 

dUo.o  ±  U.4 

A  -7+1.3 

9 
£t 

904 

QOQ  1 

oUo. 

5.2 

Q 
O 

1  (  D 

oUo.D 

4.2 

zuy 

oUo.y 

5.3 

5 

194 

308.7 

4.9 

6 

181 

308.7 

4.4 

Ac+i(2625) 

1 

431  ±  23 

341.80  ±0.10 

2 

428 

341.78 

3 

428 

341.81 

4 

433 

341.84 

7 

429 

341.77 

8 

431 

341.86 

9 

475l2289 

341.84  ±0.11 

0.7  ±0.3 

4.3.2  A+ttV 

To  search  for  A+  -)•  A+7r°7r°  signals,  we  combined  each  A+  candidate  with  two 
7r°'s  and  measured  the  mass  difference  with  Ac+.  tt0  selection  criteria  are  described 
in  the  last  chapter.  However,  to  suppress  large  background  we  need  an  additional 
cut  on  tt0  momentum.  We  chose  the  cut  based  on  the  optimization  of  the  Monte 
Carlo  signal  versus  the  background  from  the  data.  As  Table  4.2  shows,  the  Monte 
Carlo  predicts  that  150  MeV/c  gives  best  signal  to  background  ratio  for  different 
momentum  cuts.  By  fitting  the  signals  in  the  M(A+7r°7r°)  -  M(AC+)  spectrum  of  the 
data  using  different  P(n°)  cuts,  we  have  found  the  results  are  consistent  with  the 
Monte  Carlo  efficiencies. 

The  spectrum  and  its  fit  of  the  mass  difference  of  Ac+7r°7r°  combinations  are 
shown  in  Figure  4.3.  As  we  expected,  two  A+  signals  present  at  309  MeV/c2  and 
342  MeV/c2.  Unlike  the  A+tt+tt"  case  in  which  we  see  a  much  larger  A+(2625) 
signal  than  A+1(2593)'s,  the  two  signals  are  roughly  the  same  size.  The  Ac+X(2593) 
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Table  4.2:  Signal  versus  background  of  A+7r°7r°  combinations  for  different  7r°  mo- 
mentum cuts.  Si  and  S2  in  the  table  are  the  fitted  areas  of  the  A+^2593)  and  the 
A+^2625),  respectively,  decaying  into  A+7r°7r°  from  Monte  Carlo  studies;  while  B  is 
the  background  from  plotting  A+7r°7r°  combinations  of  real  data. 


P{ir°)  (GeV/c2) 

Si 

s2 

B 

S2/y/B 

0.50 

747 

661 

6480 

9.3 

8.2 

0.10 

697 

575 

4380 

10.5 

8.7 

0.15 

518 

424 

2290 

10.8 

8.9 

0.20 

279 

277 

1078 

8.5 

7.2 

0.25 

69 

109 

414 

3.4 

5.4 

signal  is  three  times  as  large  as  our  estimation  from  Monte  Carlo  efficiencies  and 
isospin  symmetry.  We  will  provide  an  explanation  in  the  next  chapter. 

To  fit  the  plot,  we  have  used  the  same  functions  as  in  the  A+7T7T  mass-difference 
spectrum.  Also  shown  in  the  plot  as  a  histogram  is  the  mass-difference  spectrum 
obtained  by  using  A+  sidebands  instead  of  A+  candidates.  Sidebands  are  picked 
from  3.4  to  5.0  oM/^+  on  both  side  of  A+  mass.  No  evidence  of  peaks  can  be  seen* 
in  this  A+  sideband  mass-difference  distribution  and  its  shape  is  consistent  with  the 
M(A+7r°7r°)  -  M(A+)  background. 

We  have  fitted  Figure  4.3  in  many  different  ways.  Table  4.3  shows  the  results 
from  some  different  fitting  procedures. 

Fit  Type  1  uses  a  Breit-Wigner  convoluted  with  a  Gaussian  with  a  fixed  width 
of  2.4  MeV/c2  to  fit  the  A+1(2593)  signal  and  a  single  Gaussian  with  a  fixed  width 
of  2.9  MeV/c2  to  fit  the  A+(2625)  signal.  Both  Gaussian  widths  are  determined 
by  Monte  Carlo  studies.  Type  2  fixes  the  width  of  the  A+x(2593)  signal  using  the 
measured  value  from  the  A+tttt  decay,  with  the  Gaussian  width  floated.  Types  3 
and  4  only  fit  the  A+^2593)  signal.  Type  3  fits  the  region  from  290  MeV/c2  to  320 
MeV/c2,  and  Type  4  excludes  the  region  from  230  MeV/c2  to  355  MeV/c2.  Type  6 
-Q^tTwMVI       fiXed  maSS6S        WidthS'  WG  find  yidds  °f  23  ±  14  at  307  MeV/c2  and 
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Figure  4.3:  Mass  difference  M(A+7r°7r0)  -  M(A+)  together  with  the  rescaled  his- 
togram which  is  the  mass  difference  M(A+7r°7r°)  -  M("A+  side"),  where  "A+  side" 
is  the  same  as  the  one  in  Figure  4.2. 
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Table  4.3:  The  yield,  mass  difference(A(A/)),  Breit-Wigner  width(r),  and  Gaussian 
width(cr)  of  signal  events  from  various  fits  to  \^  ->  A+7r°7r°.  The  numbers  labeled 
by  *  are  fixed  values  in  the  fit. 


rit  lype 

\  leld 

A(M)(MeV/cr) 

r(MeV/c2) 

a(MeV/c2) 

A£(2593) 

1 

99+^ 

—  ZD 

306.3  ±  0.7 

2  9+on 

^••3— 2.0 

2.4* 

2 

113+H 

306.4  ±  0.7 

4.7* 

2  3+11 
2.4* 

3 

129^ 

306.6  ±  0.8 

5.8±1;I 

4 

109ta! 

306.2  ±  0.7 

q  0+2.8 

2.4* 

5 

102^ 

306.2  ±0.7 

o  -,  +2.7 
0\1_21 

2.4* 

Ac+x(2625) 

1 

105  ±22 

342.9  ±0.9 

2.9* 

5 

342.2  ±  1.1 

4.61};? 

6 

132  ±27 

342.7+?;jj 

1.5* 

2.9* 

7 

106  ±  24 

343.0  ±  0.9 

2.9* 

8 

107  ±23 

342.8  ±0.9 

2.9* 

and  further  only  fit  the  A^i(2625)  signal.  Type  6  fits  the  signal  with  a  Breit-Wigner 
with  a  fixed  width  of  1.5  MeV/c2  convoluted  with  a  Gaussian  with  a  fixed  with  of 
2.9  MeV/c2.*  Type  7  fits  the  region  to  the  right  of  320  MeV/c2,  and  Type  8  excludes 
the  region  from  300  MeV/c2  to  315  MeV/c2. 

Many  other  variations  of  fits  like  these  are  also  performed.  We  find  that  using 
higher-order  polynomials  to  fit  the  background  gives  almost  the  same  results.  Also, 
results  are  consistent  with  the  ones  in  Table  4.3  if  we  fit  only  one  peak  by  excluding 
the  region  of  the  other  peak  in  the  plot. 

These  fitting  results  show  very  strong  evidence  that  the  decay  mode  A+x  ->• 
A+ttV  does  exist  for  both  the  A+(2593)  and  the  A+(2625).  The  next  point  of 
interest  is  to  calculate  the  branching  ratio  B(A^  ->•  A+7r°7r°)/S(A+1  ->  A+tt+tt-). 
To  do  this  we  need  to  estimate  the  errors  on  the  yields  of  A+ttV  signals. 

We  not  only  have  large  statistical  errors  because  of  the  poor  signal,  but  also  large 
systematic  errors,  mostly  due  to  the  uncertainties  of  the  resolutions  and  the  natural 

*The  reason  we  do  this  fit  is  that  we  do  not  trust  the  result  on  fit  Type  5  since  the  width  is  too 
much  wjder  than  the  Monte  Carlo  predicted  resolution.  Here  we  use  this  fit  trying  to  estimate  the 
upper  limit  of  the  yield 
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Table  4.4:  Comparison  of  productions  of  A+7r°7r°  for  different  momentum  cuts. 
The  A^(2593)  yields  are  from  the  fits  using  a  Breit-Wigner  function  with  T  =  4.7 
MeV/c2  convoluted  with  a  Gaussian  with  fixed  width  from  Monte  Carlo  studies. 
The  A+j(2625)  yields  are  from  fits  using  a  single  Gaussian  with  fixed  width  from  the 
Monte  Carlo.  The  statistical  errors  on  efficiencies  are  tiny. 


P,o(GeV) 

Yield 

Efficiency(%) 

Yield/Efficiency(xl00) 

A+(2593) 

0.05 

133  ±  30 

2.35 

57  ±13 

0.10 

145  ±  28 

2.19 

66  ±  13 

0.15 

115  ±20 

1.63 

70  ±  12 

0.20 

59  ±  13 

0.88 

67  ±  15 

0.25 

11  ±9 

0.22 

50  ±41 

Ac+i(2625) 

0.05 

107  ±34 

2.46 

43  ±  14 

0.10 

95  ±30 

2.31 

41  ±  13 

0.15 

105  ±  22 

1.78 

59  ±  12 

0.20 

65  ±  17 

0.95 

68  ±  18 

0.25 

38  ±  11 

0.37 

103  ±  30 

width  of  the  particle.  Thus  we  have  fitted  the  signals  by  varying  the  resolution 
by  ±20%  of  the  Monte  Carlo  predicted  values,  varying  the  width  of  the  A;!i(2593) 
from  2  MeV/c2  to  7  MeV/c2,  and  varying  the  width  of  the  A+j(2625)  from  0  to  1.5 
MeV/c2.  We  then  estimate  that  the  fitting  uncertainties  of  the  yields  are  +g  for  the 
A+(2593)  and  tU  for  the  A+(2625). 

Since  a  150  MeV/c  momentum  cut  on  both  independent  7r°'s  is  a  very  hard  cut,  a 
consistency  check  has  been  made  by  comparing  the  efficiency-corrected  productions 
for  different  tt°  momentum  cuts.  As  shown  in  Table  4.4,  for  the  A+j(2593)  decay, 
the  Monte  Carlo  and  the  data  agree  very  well,  but  for  the  A+(2625),  we  see  some 
disgreement.  Thus  we  need  to  intoduce  a  corresponding  systematic  error  when  we 
calculate  the  decay  ratios. 


66 


4.4  Substructures 

Studying  A^-^A+tttt  decay  substructures  through  intermediate  Ec's  is  very  use- 
ful to  identify  Jp  of  the  A^'s  and  gain  more  understanding  of  these  two  states.  In 
this  section,  we  describe  the  measurements  on  the  number  of  A+j's  which  decay  to 
Ec's.  The  ratio  is  defined  to  be 


hc  =  B(A+  ->  Ectt)/£(A+  ->  A+tttt). 

4.4.1    A+  -»  S++7T-  and  A+  -»  E°tt+ 

Figure  4.4  plots  the  mass  difference  M(A+ti-)  -  M(A+)  against  M(A+tt+)  - 
M(A+)  for  the  same  A+tt+tt-  combinations  in  Figure  4.1.  By  looking  directly  at  the 
figure,  it  is  easy  to  tell  that  there  are  two  very  prominent  clusters  which  correspond 
to  the  decays  through  intermediate  E++  and  E°at  the  two  ends  of  the  diagonal  band 
of  A+1(2593)'s  in  which  all  entries  are  shown  as  inverted  triangles;  in  the  diagonal 
band  of  A+1(2625)'s  in  which  all  entries  are  shown  as  triangles,  no  obvious  clusters 
can  be  seen.  Therefore,  by  simply  inspecting  this  scatter  plot,  we  can  tell  that 
the  Ectt  mode  dominates  the  A+(2593)  decay,  and  the  non-resonant  A+tttt  mode 
dominates  the  A+^2625)  decay. 

To  determine  the  fraction  /E++  and  /Eo  of  the  A+(2593),  we  first  measured 
the  A+(2593)  yield  as  a  function  of  A+TT+  submass.  Figure  4.5  shows  the  mass 
enhancement  of  E?  and  Ec++  from  A+tt+tt"  combinations  with  the  mass  difference 
M(A+tt+tt-)  -  M(A+)  between  304  MeV/c2  and  314  MeV/c2. 

We  fit  each  plot  in  Figure  4.5  with  a  bifurcated  Gaussian*  (left)  and  a  single 
Gaussian  (right)  together  with  the  same  background  function.  The  same  combina- 
tions in  the  two  plots  are  constrainted  to  have  the  same  area. 

*A  bifurcated  Gaussian  is  A  Gaussian  with  different  slopes  on  either  side  of  the  mean. 
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Figure  4.4:  Scatter  plot  of  the  mass  difference  M(A+7i-)-Af  (A+)  vs.  the  mass  differ- 
ence M (A+7T+)  -  M (A+).  Shown  in  the  plot  are  (1)  the  vertical  band  at  A+tt+  mass 
difference  of  about  168  MeV/c2,  corresponding  to  E++  ->  A+7T+;  (2)  the  horizontal 
band  at  A+tt  mass  difference  of  about  168  MeV/c2,  corresponding  to  E°  ->  A+tt^- 
(3)  the  diagonal  band  extending  from  (0.14,  0.20)  to  (0.20,  0.14),  corresponding  to 
combinations  from  the  A+(2625)  signal  region  in  Figure  4.1  (triangles);  and  (4)  the 
diagonal  band  extending  from  (0.14,  0.168)  to  (0.168,  0.14),  corresponding  to  com- 
binations from  the  A+(2593)  signal  region  in  4.1  (inverted  triangles),  within  which 
there  are  two  lobes  at  the  kinematic  limits:  (a)  the  enhancement  in  the  £++  band 
just  above  the  threshold  for  A+7T+,  and  (b)  the  enhancement  in  the  S°  band  just 
above  the  threshold  for  A+7r~. 
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Figure  4.5:  Events  and  fits  of  £°'s  and  E++'s  with  the  mass  difference  of 
M (A+tt+tt  )  -  M (A+)  between  304  MeV/c2  and  314  MeV/c2  (the  A+(2593)  signal). 
Plot  (a)  shows  the  mass  difference  M{^~)  -  Af  (A+).  The  right  peak  centered  at 
167  MeV/c2  is  the  £°  signal;  and  the  left  peak  is  the  reflection  of  the  £++  signal 
on  the  A;?!  (2593)  band.  Plot  (b)  shows  the  mass  difference  M(A+?r+)  -  M (A+) 
The  right  peak  centered  at  167  MeV/c2  is  the  £++  signal;  and  the'left  peak  is  the 
reflection  of  the  £°  signal  on  the  A+(2593)  band.  Since  plots  (a)  and  (b)  represent 
the  same  events  projected  in  two  different  directions,  they  are  fitted  simultaneously 
by  constraining  the  area  of  the  left/right  peak  in  (a)  to  be  the  same  as  the  area 
of  the  right/left  peak  in  (b)  The  left  peaks  in  both  (a)  and  (b)  are  fitted  by  bifur- 
cated Gaussians,  and  the  right  peaks  are  fitted  by  single  Gaussians.  They  are  all 
constrainted  to  have  the  same  widths. 
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Shapes  of  the  bifurcated  Gaussian  and  the  single  Gaussian  are  determined  by 
Monte  Carlo  studies  of  A^Ecn  and  Ec->-A+7r  decays.  The  background  is  expected 
to  be  a  mixture  of  two  parts:  random  A+7T+7T-  combinations  and  non-resonant 
A^(2593)^A+7r+7r~  events.  Alhough  it  is  very  difficult  to  predict  what  its  exact 
shape  is,  we  know  this  shape  has  cut-offs  at  phase  space  lower  limit  139  MeV/c2 
and  upper  limit  about  172  MeV/c2,*.  We  used  a  "double-threshold"  function  y(x) 
with  upper  and  lower  thresholds  at  a  and  b  (in  this  particular  case  a  =  139  MeV/c2 
and  b  =  172  MeV/c2), 


y(x)  =  N(x  -  a)(b  -  x)eMx-a)+S2(x-af es3(b-x)+s4(b-x)^ 


(4.1) 


to  fit  the  background  of  Figure  4.5.  Here  N  is  a  normalization  factor  and  si5  where  i 
=  1,  2,  3,  and  4,  determines  the  shape  of  the  function.  The  fitted  background  shape 
obtained  by  fitting  the  plot  in  Figure  4.5  with  floating  s^s  is  very  consistent  with  the 
shape  from  the  Monte  Carlo  generated  A+tt+tt'  non-resonance  decay.  As  we  tested, 
the  fitted  results  are  not  sensitive  to  the  variation  of  the  background  shape  as  far  as  a 
and  b  are  fixed.  Reasonably  varying  the  parameters  of  the  double-threshold  function 
(reasonably  means  still  giving  a  good  fit)  makes  the  yields  of  signals  change  within 
5%  of  the  mean  value  of  fitted  areas.  We  believe  that  statistical  errors  dominate  the 
uncertainties  of  this  measurement. 

The  areas  of  E°  and  E++  signal  we  measured  from  Figure  4.5  are  77±7  and  72±7, 
respectively.  These  values  are  obtained  by  fitting  Figure  4.5  by  fixing  the  back- 
ground, with  shape  determined  by  the  Monte  Carlo  studies  of  A+1(2593)^A+7r+7r- 
non-resonance  decay.  To  obtain  the  number  of  A+(2593)'s  which  decay  into  Ectt, 
these  numbers  need  to  be  subtracted  by  the  number  of  "independent"  Ec's  in  the 

'The  lower  limit  is  really  the  pion  mass,  and  the  upper  limit,  as  we  can  tell  in  FiKure  4  4 
depends  on  the  width  of  the  A+  (2593)  band. 
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sideband  of  the  A^(2593)  signal.  Due  to  the  kinematic  limit,  only  the  higher-mass 
A^(2593)  sideband  contains  Ec's.  Therefore  we  estimated  the  background  as  being 
half. 

The  number  of  Ec's  in  the  A^(2593)  sideband  were  measured  similarly  as  we 
did  with  the  number  of  Ec's  in  the  A^(2593)  peak  except  that  the  mass-difference 
range  is  wider  and  the  kinematic  reflections  of  the  peaks  are  in  different  mass- 
difference  ranges.  The  sidebands  were  taken  from  the  range  of  the  mass  difference 
M(A+7r+7r-)  -  M(A+)  between  320  MeV/c2  and  330  MeV/c2.  Figure  4.6  shows  the 
histograms  and  the  fits.  The  fitted  yields  are  19±4  of  E++'s  and  15±4  of  E°'s. 

After  subtracting  one-half  of  the  yield  from  the  sideband,  we  obtained  the  num- 
ber of  A;!j(2593)'s  which  decay  to  Ec's:  67±8  for  the  E++  and  65±8  for  the  E° 
Compared  with  the  total  number  of  A+1(2593)'s  in  the  mass-difference  range  be- 
tween 304  MeV/c2  and  314  MeV/c2,  which  is  137±22  measured  in  Figure  4.1,  we 
see  that  resonance  via  the  Ec  dominates  these  transitions. 

The  same  methods  were  applied  to  measure  /Ec's  of  the  A+^2625)  except  that 
in  this  case  there  are  two  A+1(2625)  sidebands  other  than  the  one  in  the  A+^2593) 
case.  The  numbers  of  E++'s  and  E°'s  inside  the  A+^2625)  signal  are  measured  as 
shown  in  Figure  4.7,  and  the  fitted  yields  are  46±9  for  E++  and  36±9  for  E°  The 
lower  A+!(2625)  sideband  is  chosen  to  be  between  320  MeV/c2  and  330  MeV/c2, 
which  is  the  same  as  the  one  shown  in  Figure  4.6,  and  there  are  19±4  E++'s  and 
15±4  E°'s.  The  higher  sideband  is  therefore  in  the  range  between  354  MeV/c2 
and  364  MeV/c2,  and  the  number  of  E++'s  and  E°'s  is  measured  to  be  13±3  and 
6±4,  respectively.  The  fit  for  the  higher  sideband  is  shown  in  Figure  4.8.  After 
the  sideband  subtraction,  the  results  are  30±10  A^'s  decay  to  E++'s  and  25±10 
decay  to  E°'s.  These  numbers  indicate  that  resonance  through  Ec  is  very  small 
in  the  decay  of  the  A+1(2625),  since  the  total  number  of  A+(2625)'s  between  mass 
difference  354  MeV/c2  and  364  MeV/c2,  450±29,  is  much  larger. 
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Figure  4.6:  Events  and  the  fits  of  £°'s  and  £++'s  in  the  A^(2593)  sideband.  Plot 
(a)  shows  the  mass  difference  M(A+tt~)  -  M(A+),  and  plot  (b)  shows  M(A+tt+)  - 
M(A+).  The  peaks  in  the  right  at  about  168  MeV/c2  on  both  plots  are  the  signals 
of  the  £°  and  the  £++.  In  each  plot  a  bump  is  centered  at  about  158  MeV/c2, 
representing  reflections  of  the  £°(£++)  peak  seen  in  the  direction  of  £++(£°).  Each 
plot  is  fitted  by  two  single  Gaussians  with  fixed  widths  (3.4  MeV/c2  of  the  left  bump 
and  1.1  MeV/c2  of  the  right  peak)  from  £°  ->  A+tT  (£++  ->  A+7T+)  Monte  Carlo 
studies  and  a  double-threshold  function  the  same  as  the  one  used  to  fit  Figure  4.5, 
with  the  lower  threshold  140  MeV/c2  and  the  higher  threshold  185  MeV/c2.  The 
areas  representing  the  same  combinations  are  constrainted  to  be  the  same,  and  the 
two  plots  are  constrained  to  have  an  identical  background. 


72 


0.15  0.16  0.17  0.18  0.19 

GeV/c2 

Figure  4.7:  The  mass  difference  of  M(A+7r+7r_)  -  M(A+)  between  337  MeV/c2  and 
347  MeV/c2  (the  A+1(2625)  signal).  The  meanings  of  the  signals  are  the  analog  to 
those  in  Figure  4.5.  The  reflection  of  the  peak  of  one  combination  is  overlaid  with 
the  peak  of  the  other  combination.  These  shapes  are  predicted  by  A+1(2625)->Ec7r 
Monte  Carlo  studies.  Each  plot  is  fitted  by  two  single  Gaussians  and  a  second-order 
Chebyshev  polynomial.  The  widths  of  the  Gaussians  are  fixed  by  the  values  from 
the  Monte  Carlo  (1.2  MeV/c2  of  the  peak  in  the  left  and  3.4  MeV/c2  of  the  bump 
in  the  right).  The  areas  representing  the  same  combinations  are  constrained  to  be 
the  same,  and  the  two  plots  are  constrained  to  have  the  same  background. 
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Figure  4.8:  The  mass  difference  of  M(A+7r+7r-)  -  M(A+)  between  354  MeV/c2  and 
364  MeV/c2  (the  A+(2625)  sideband),  fitted  in  the  same  way  as  fitting  Figure  4.7. 
The  widths  of  the  Gaussians  are  fixed  by  the  values  from  the  £c->-A+7r  Monte  Carlo 
studies  (1.3  MeV/c2  in  the  left  and  4  MeV/c2  in  the  right). 


After  further  considerations  of  systematic  errors,  we  listed  our  results  of 
and  /Eo,  compared  with  the  measurements  with  other  experiments,  in  Table  4.13 
and  Table  4.14. 

4.4.2    A+  ->  £+7T° 

The  technique  of  measuring 

/Sc+  =  B(A+  ->  J£**)/B(AS  ->  ACW), 

is  similar  to  that  of  measuring  and  /Eo,  and  much  simpler.  Since  A+7r°7r°  decay 
only  has  one  substructure,  which  is  A+7r°,  instead  of  two  (£++7r~  and  £°tt+),  we 
do  not  need  to  do  complicated  fits  to  correlate  £++  and  £°  signals  as  we  did  in  the 
A+7r+7r~  case. 

The  scatter  plot  of  M(A+7r°7r°)  -  M(A+)  versus  Af(A^"7r°)  —  M(A^)  shown  in 
Figure  4.9  shows  the  evidence  of  A+1(2593)^-E+7r°  substructure.  In  the  plot,  we 
can  see  two  concentrated  clusters  of  events  in  the  A^i(2593)  mass-difference  band 
around  306  MeV/c2:  one  at  about  167  MeVc2  which  is  consistent  with  £+  mass 
difference,  and  another  which  is  the  kinematic  reflection  of  the  first,  at  about  140 
MeV/ c2  just  above  the  kinematic  threshold.  We  cannot  see  this  phenomenon  for 
the  A+1(2625)  signal. 

Low  statistics  and  high  background  do  not  allow  us  to  report  the  A+1(2625)->£c7r° 
measurement  (the  upper  limit  would  be  1).  Thus  only  a  /E+(A+1(2593))  measure- 
ment was  performed. 

Figure  4.10  shows  the  histogram  and  the  fit  of  the  signal  and  its  kinematic  re- 
flection at  the  A+!(2593)  mass-difference  band  taken  from  the  M(A+7r°7r°)  -  M(A+) 
mass-difference  range  between  304  and  314  MeV/c2.  The  plot  was  fitted  by  two 
single  Gaussians  with  fixed  widths  (2.9  MeV/c2  for  the  left  signal  and  2.4  MeV/c2 


75 


0.27  I — ^ — 1 — 1 — 1 — 1 — ' — i — i — I — i — i — i  i  I  ■    ■    ■    .  I 

0.130  0.155  0.180  0.205  0  230 

M(ACV)  -  M(AC+)  (GeV/c2) 


Figure  4.9:  Scatter  plot  of  the  mass  difference  M(A+7r°7r°)  -  M(A+)  vs.  the  mass 
difference  M (A+tt°)  -  M(A+).  Very  visible  in  the  figure  is  a  horizontal  band  at  about 
306  MeV/c2  of  M(A+7r°7r°)  -M(A+),  corresponding  to  A+(2593)  -+  A+ttV.  It  has 
two  lobes  at  the  kinematic  limits:  one  at  about  167  MeV/c2  of  M(A+tt°)  -  M(A+), 
which  is  consistent  with  £+  mass  difference;  and  the  other  just  above  the  kinematic 
threshold  which  is  tt0  mass.  We  assume  there  are  three  other  bands  in  the  plot 
(visible  but  not  obvious)  which  are  not  shown  clearly:  (1)  a  vertical  band  at  about 
167  MeV/c2  of  M(A+tt°)  -  M(A+),  corresponding  to  E+  ->  A+tt0;  (2)  a  diagonal 
band  extending  from  (0.135,  0.308)  to  (0.2,  0.37),  which  is  the  kinematic  reflection 
from  band  (1),  corresponding  to  the  single  7r°'s  from  the  random  E+tt°  combinations; 
(3)  a  horizontal  band  at  about  342  MeV/c2  of  M(A+7r°7r°)  -  M(A+),  in  which  no 
E+  clustering  can  be  seen,  corresponding  to  A+j(2625)  ->  A+7r°7r°. 
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Figure  4.10:  The  mass  difference  M(A+7r°7r°)  -  M(A+n°)  between  304  MeV/c2  and 
314  MeV/c2  (the  A+j(2593)  signal).  The  two  peaks  correspond  to  E+'s  (left)  and 
independent  7r°'s  (right)  in  the  A+j(2593)  ->  E+tt°  decay. 
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Figure  4.11:  The  mass  difference  M(A+7r°7r°)  -  M(A+7r°)  between  320  MeV/c2  and 
330  MeV/c2  (the  A+^2593)  sideband).  The  positions  and  shapes  of  the  two  bumps 
are  determined  by  the  E+->A+7r°  Monte  Carlo.  The  signal  at  right  corresponds  to 
E+^A+tt0,  and  the  left  bump  reflects  single  tt0's  from  the  E+tt°  combinations. 

for  the  right  signal)  determined  by  the  A+1(2593)^S+tt°  Monte  Carlo,  and  a  double- 
threshold  function  which  is  defined  in  equation  4.1.  The  shape  of  the  double- 
threshold  function  with  the  two  thresholds  135  MeV/c2  and  176  MeV/c2  is  de- 
termined by  the  non-resonance  A+1(2593)-^A+7r°7r°  Monte  Carlo.  The  areas  of  the 
the  two  Gaussians  are  constrained  to  be  the  same.  Finally,  the  fitted  yield  is  89±11. 
Consistent  fitting  results  were  obtained  when  we  fitted  the  background  with  a  float- 
ing shape  of  the  double-threshold  function. 

Figure  4.11  shows  the  M(A+tt0)  -  M(A+)  distribution  at  the  A+(2593)  side 
band  and  the  fit  to  it.  The  sideband  was  taken  from  the  M(A+7r°7r°)  -  M(A+) 
range  between  320  and  330  MeV/c2.  The  plot  is  fitted  by  two  single  Gaussians  with 


fixed  mean  values  and  shapes  determined  by  the  £+— >A+7T°  Monte  Carlo,  and  a 
double-threshold  function  with  fixed  thresholds  and  other  floating  parameters.  The 
areas  of  the  these  two  Gaussians  are  constrained  to  be  the  same.  The  fitted  yield  is 
Mill. 

We  finally  calculated  that  81±16  of  A+1(2593)  decays  are  to  £+'s  from  a  total 
number  of  93±20  A^i(2593)'s  in  the  mass-difference  range  between  304  MeV/c2  and 
314  MeV/c2.  This  is  consistent  with  what  we  have  seen  in  the  A^(2593)->A+7r7r 
decay,  for  which  the  Ec7r  substructure  dominates.  The  result  of  /E+(A+1(2593))  is 
presented  in  Table  4.12  after  proper  systematic  errors  are  evaluated. 

4.5    Search  for  Other  Decays 

Since  isospin  is  conserved  in  strong  decays,  the  A+x  ->•  A+7T0  is  forbidden.  If 
any  of  the  two  A+7r+7r~  states  we  observed  were  excited  Ec  baryons,  A+7T  decay 
should  be  not  only  allowed  but  more  favorable  than  A+7T7T  decays  due  to  the  larger 
phase  space.  A^'s  are  allowed  to  decay  electromagnetically  to  A+7,  and  if  any 
A^i's  intrinsic  width  is  sufficiently  narrow,  its  decay  may  be  competitive  with  A+7T7T 
decays. 

To  search  for  these  decays,  we  reconstructed  A+7r°  and  A+7  combination  using 
our  previously  defined  candidates  of  A+,  tt°  and  7.  Based  on  the  optimization  of 
Monte  Carlo  signals  versus  the  background  from  data,  a  200  MeV/c  n°  momentum 
cut  for  A+7r°  and  a  300  MeV/c  photon  momentum  cut  for  A+7  are  used. 

Figure  4.12  shows  the  mass-difference  spectrum  of  M(A+tt0)  -  M(A+),  and  no 
signal  can  be  observed.  To  measure  the  upper  limit  of  decay  ratios,  three  functions 
are  used  to  fit  the  distribution.  A  Breit-Wigner,  with  fixed  4.7  MeV/c2  A+j(2593) 
width,  convoluted  with  a  Gaussian  with  fixed  5.7  MeV/c2  Monte  Carlo  predicted 
resolution,  is  used  to  fit  at  a  fixed  mean  value  of  307  MeV/c2;  a  single  Gaussian 
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Figure  4.12:  Mass  difference  of  M(A+tt°)  -  M(A+). 
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Figure  4.13:  Mass  difference  of  M(A+7)  -  M(A+). 
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Table  4.5:  Fitting  results  on  M(A+tt°)-M(A+)  and  M(A+7)-M(A+)  distributions. 


Decay 

Yield 

Yield  Upper  Limit  (90%  C.  L.) 

A+(2593)->  A+tt° 
A+(2625)->  A+tt° 
A+(2593)^Ac+7 
Ac+1(2625)^A+7 

26  ±28 
-8  ±36 
96  ±71 
-42  ±  52 

160 
165 
381 
120 

with  6.7  MeV/c2  Monte  Carlo  predicted  resolution  is  used  to  fit  at  a  fixed  mean 
value  of  342  MeV/c2;  and  the  background  is  fitted  with  a  second-order  polynomial. 

Also,  no  signal  is  found  at  mass-difference  spectrum  of  M(A+7)  —  M(A+), 
as  shown  in  Figure  4.13.  We  used  the  exact-same  fitting  procedures  as  in  the 
M(A+7r0)  — M(A+)  case,  except  different  Gaussian  resolutions,  which  are  7.4  MeV/c2 
at  mass  difference  307  MeV/c2  and  7.9  MeV/c2  at  mass  difference  342  MeV/c2. 

The  fitted  yields  and  their  projected  upper  limits  are  shown  in  Table  4.5.  We 
conclude  that  we  have  observed  no  evidence  of  A+7r°  or  A+7  decays  from  A^. 

4.6    Masses  and  Widths 

4.6.1    Systematic  Uncertainties 

One  major  part  of  the  systematic  error  of  mass  measurements  is  the  fitting 
error,  mostly  from  the  uncertainty  of  the  background  shape  or  sometimes  even  the 
signal  shape.  To  evaluate  this  error,  we  use  many  different  fit  types  and  check 
the  variations  of  the  mean  value  of  the  signal.  For  a  particle  with  a  considerably 
wide  natural  width,  such  as  the  or  a  particle  with  an  uncertain  shape  of  the 

signal,  such  as  the  A+x(2593),  the  fitting  error  generally  dominates  the  systematic 
uncertainty. 
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Since  we  measure  the  mass  difference  between  the  A+7T7T  and  the  A+  instead 
of  measuring  the  directly,  most  of  the  experimental  uncertainties  from  the  re- 
construction of  the  A+  cancel.  As  we  checked,  3  MeV/c2  change  of  the  A+  mass 
shifts  the  M(A+7r+7r~)  -  M(A+)  by  no  more  than  0.01  MeV/c2.  So  the  systematic 
uncertainties  other  than  the  fitting  errors  mostly  come  from  the  measurement  of  the 
momenta  and  decay  angles  of  the  two  extra  pions.  For  charged  pions,  these  errors 
are  dominated  by  the  uncertainties  in  the  magnetic  field  normalization  and  the  en- 
ergy loss  correction  which  is  applied  to  the  tracks  traversing  the  beam  pipe  and  the 
drift  chambers.  For  7r°'s,  the  errors  come  from  the  shower  energy  measurement  by 
the  crystal  calorimeters. 

To  evaluate  the  systematic  errors  of  mass-difference  measurements  caused  by  the 
independent  charged  pions  with  low  momentum,  we  measured  the  mass  difference 
M(D*+)  -  M(D°)  in  decay  D*+  ->  D°ir+  since  the  mass-difference  measuring  tech- 
nique is  the  same  as  M(A+7r+7r_)  -  M(A+)  and  the  independent  pion  is  soft.  Our 
fitted  result  is  145.39±0.01  MeV/c2,  in  excellent  agreement  with  the  PDG  value, 
which  is  145.397±0.030  MeV/c2.  We  found  that  0.2%  momentum  change  of  all 
charged  tracks  will  shift  the  M(D*+)  -  M(D°)  by  about  0.05  MeV/c2,  larger  than 
the  error  of  the  PDG  average.  Therefore  we  believe  0.2%  change  in  momentum  is 
a  reasonable  amount  to  study  the  momentum  scale  in  mass-different  measurement 
with  charged  soft  pions. 

For  the  systematic  errors  caused  by  the  measurement  of  the  soft  7r°'s,  similar  to 
the  charged  ones,  we  check  the  photon  energy  scale  by  looking  at  the  mass  difference 
of  M(D*°)~M(D°),  and  our  fitted  value  agrees  very  well  with  the  PDG  value,  which 
is  142.12±0.07  MeV/c2.  We  think  1%  energy  change  is  a  reasonable  amount  used 
for  an  energy-scale  study  since  it  shifts  M(D*°)  -  M(D°)  by  about  0.10  MeV/c2. 

In  the  study  of  M{D*+)  -  M(D°)  measurement,  we  find  that  mass  difference 
is  more  sensitive  to  the  systematic  shifts  of  the  position  and  angle  parameters  of 
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tracks.  However,  the  charged  tracks  are  measured  much  more  accurately.  Our  ver- 
tex constraint  program  would  shift  the  track  positions  and  angles  by  fair  amounts. 
We  find  that  calculating  the  mass  difference  by  constraining  the  D°  and  the  pion 
to  be  from  the  event  vertex,  gives  up  to  0.02±0.01  MeV/c2  mean- value  difference, 
compared  with  the  mean  value  of  mass  difference  reconstructed  without  vertex  con- 
straint. This  difference  reasonably  reflects  the  systematic  error  from  the  decay  angle 
measurement. 

The  systematic  error  on  the  width  measurement  comes  from  two  sources:  the 
fitting  error  and  the  mass  resolution  of  the  detectors  determined  by  Monte  Carlo 
studies.  The  fitting  error  dominates  the  total  error  when  the  natural  width  of 
a  particle  is  wide,  and  it  can  be  evaluated  by  varying  parameterizations  of  the 
background  and  the  signal  shapes.  The  detector  resolution  becomes  important  when 
it  is  close  to,  or  bigger  than  the  natural  width,  and  its  correctness  can  be  checked 
by  fitting  the  signals  of  other  particles  with  well-known  widths  and  comparing  the 
results  with  corresponding  Monte  Carlo  signals.  As  we  tested,  the  CLEO  Monte 
Carlo  program  does  a  very  good  job  with  tracks  with  medium-to-high  momentum 
spectra,  but  for  tracks  with  low-momentum  spectra,  it  tends  to  overestimate  the 
detector  performance  and  therefore  predicts  a  mass  resolution  with  a  lower  value 
than  the  real  resolution.  But  by  the  time  this  research  work  was  performed,  we  are 
unable  to  make  corrections,  since  it  is  very  difficult  to  detect  what  the  exact  reason 
is  and  by  exactly  how  much  the  Monte  Carlo  is  wrong.  So  we  have  to  be  very 
careful  and  conservative  with  our  estimates  of  the  uncertainties  on  narrow-width 
measurements. 
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4.6.2    A+j(2625)  Mass  and  Width  Limit 

The  result  presented  on  the  A^j(2625)  mass  difference  with  the  A+,  341.80±0.10 
MeV/c2  is  our  measured  mean  value  by  fitting  the  A+7r+7r~  signal  with  a  width- 
float  single  Gaussian.  By  fitting  the  A+7r°7r°  signal,  we  obtained  a  value  of  342.9  ± 
0.9  MeV/c2.  We  have  also  checked  the  fits  to  the  A+7r+7r~  signal  using  a  double 
Gaussian,  a  bifurcated  Gaussian,  and  a  Breit-Wigner  convoluted  with  a  Gaussian, 
with  variations  of  different  shapes.  All  fitted  means  agree  within  0.13  MeV/c2. 

In  the  momentum-scale  study,  we  find  that  a  0.2%  momentum  change  for  all 
charged  tracks  shifts  M(A+1(2625))  -  M(A+)  by  about  0.22  MeV/c2.  We  find  that 
the  mass  difference  M(A+1(2625))  -  M(A+)  shifts  for  about  0.2±0.1  MeV/c2  if  we 
do  not  constrain  the  A+  and  the  two  pions  to  be  from  a  single  vertex  compared 
with  constraining  them  into  the  main  vertex.  This  proved  that  the  uncertainty  of 
the  pion  decay  angle  does  not  bring  large  systematic  errors.  We  also  compared 
M(A+7r+7r-)  -  M(A+)  with  M(A+7r+7r-)  -  M(A+),  and  their  mean  values  agree 
within  statistical  errors. 

According  to  Table  4.1,  the  fitting  error  due  to  the  background  parameteriza- 
tion with  the  signal  fitted  by  a  single  Gaussian  is  about  ±0.04  MeV/c2,  which  is 
tiny  compared  with  the  uncertainty  of  the  soft  pion  momentum  measurement.  Con- 
sidering all  these  sources  of  errors,  we  conservatively  estimate  that  the  systematic 
uncertainty  is  0.35  MeV/c2. 

With  the  width  of  0.7±0.3  MeV/c2  extracted  from  the  Breit-Wigner  convoluted 
with  a  Gaussian  fit,  considering  that  our  Monte  Carlo  predicted  resolution  should 
not  be  larger  than  the  actual  one,  we  calculated  the  upper  limit  of  the  A+1(2625) 
width  to  be  1.4  MeV/c2  at  90%  confidence  level. 
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Table  4.6:  Measurements  of  the      (2625)  mass  and  width  limit. 


Experiments 

[M(Ac+1)-M(A+)](MeV/c:i) 

rA+  (MeV/c2) 

E687[29] 
CLEO-II[30] 
ARGUS[32] 
CLEO-II(this) 

340.4±0.6±0.3 
342.2±0.2±0.5 
342.1±0.5±0.5 
341.80±0.10±0.35 

<  1.9 

<  3.2 

<  1.4 

Table  4.7:  Measurements  of  A^(2593)  mass  and  width  together  with  the  mass  res- 
olutions in  the  corresponding  experiments.  All  quantities  are  in  units  of  MeV/c2. 
Note  that  E687[31]  fitted  their  signal  using  a  single  Gaussian  and  the  fitted  Gaussian 
width  is  consistent  with  that  A+^2593)  has  a  zero  intrinsic  width. 


Experiments 

Resolution 

[M(A+)  -  M(A+)](MeV/c2) 

TA+]  (MeV/c2) 

CLEO-II[30] 

E687[31] 
ARGUS[32] 
CLEO-II(this) 

2.0 
1.8 
1.8 
1.62 

307.5±0.4±1.0 
309.2±0.7±0.3 
309.7±0.9±0.4 
308.8±0.4±0.4 

O  Q+1.4+2.0 

°-y-1.2-1.0 

o  q+2.9+1.8 
z-y-2. 1-1.4 
A  7+1.3+1.4 
4-  '-1.1-0.9 

The  results  of  mass  difference  with  the  A+  and  the  upper  limit  of  the  width  are 
listed  in  Table  4.6.  A  comparison  with  the  previous  experiment  results  is  given  in 
the  table. 


4.6.3    A+(2593)  Mass  and  Width 

In  Figure  4.1,  we  fit  the  A+1(2593)^A+7r+7r-  signal  in  the  same  manner  as  the 
previous  experiments:  using  a  Breit-Wigner  function  convoluted  with  a  Gaussian. 
The  results  of  the  mass  and  width  are  listed  in  Table  4.7,  together  with  the  mea- 
surements of  previous  experiments.  Here  the  systematic  errors  are  mostly  fitting 
errors,  evaluated  based  on  Table  4.1. 

Although  our  measurements  agree  with  previous  experiments,  we  doubt  the  re- 
sults, and  the  most  obvious  reasons  comes  from  the  following  facts:  (1)  A+x(2593) 
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has  a  nearly  100%  dominating  decay  channel,  £c7r;  (2)  the  combination  of  H°7r+  (or 
£^+7r-)  has  a  kinematic  limit  at  the  mass  difference  of  about  307  MeV/c2;  (3)  the 
fit  to  the  A^(2593)  signal  in  Figure  4.1  is  poor.  Because  of  these  and  some  other 
thoughts  which  will  be  stated  later,  we  feel  that  using  a  full  Breit-Wigner  distribu- 
tion to  parameterize  the  S°7r+  (or  Ec7r")  decay  very  close  to  the  kinematic  threshold 
may  not  be  correct.  We  will  describe  a  new  way  to  parameterize  the  signal  in  the 
next  chapter,  since  the  new  method  uses  the  £++  and  £°masses  and  widths. 


4.7    Fragmentation  Functions 

As  stated  in  Chapter  1,  all  quarks  "fragment"  into  hadrons.  The  initially  pro- 
duced pair  of  hadron  and  anti-hadron  should  be  expected  to  have  more  energy  than 
the  subsequent  hadrons  produced  further  in  the  decay  chains.  Each  hadron,  as  it 
is  produced,  carry  away  a  fraction,  z,  of  total  energy  of  the  corresponding  original 
quark.  A  fragmentation  function  (~,  where  N  is  the  number  of  hadron  produced 
from  a  type  of  specific  physical  process)  proportional  to  the  probability  a  hadron 
being  produced  in  the  interval  of  [z,  z  +  dz). 

In  e+e~  annihilation  experiments,  measurement  of  the  fragmentation  functions  of 
heavy  quarks  provides  information  about  non-perturbative  particle  production  in  a 
variety  of  experimental  environments.  Many  forms  of  functions  have  been  suggested 
to  describe  the  normalized  momentum  spectra  for  heavy  quarks.  Among  those,  the 
functional  form  given  by  Peterson  et  al.[33]  is  the  most  widely  used.  The  Peterson 
function  has  form 

dN  _  1 

~dz~~  z[l-(l/*)-eP/(l -*)]*' 

where  z  =  (E  +  p\\)  hadron /{E  +  pQ),  and  p\\  is  the  longitudinal  momentum  and  pQ 
is  the  total  momentum  of  the  heavy  quark.  The  quantity  z  is  not  experimentally 
accessible,  and  in  CLEO  experiments  a  close  approximation  is  made  by  using  the 
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Table  4.8:  Efficiencies  and  productions  of  A^'s  for  different  xp's.  The  statistical 
errors  on  efficiencies  are  negligible. 


Xp 

Efficiency  (%) 

Yield 

Normalized  Yield 

A+j(2593) 

0.5-0.6 

5.45 

34  ±  16 

0.13  ±0.06 

0  6-0  7 

6  44 

1  '-t  HZ  io 

U.iO  X  u.uo 

0.7-0.8 

7.07 

91  ±  18 

0.27  ±0.05 

0.8-0.9 

7.11 

71  ±  15 

0.21  ±0.05 

0.9-1.0 

5.47 

35  ±  11 

0.14  ±0.05 

Ac+!  (2625) 

0.5-0.6 

5.63 

114  ±  13 

0.18  ±0.02 

0.6-0.7 

6.47 

217  ±  16 

0.30  ±0.02 

0.7-0.8 

7.08 

214  ±  16 

0.27  ±0.02 

0.8-0.9 

7.78 

161  ±  14 

0.19  ±0.02 

0.9-1.0 

7.48 

48  ±8 

0.06  ±0.01 

scaling  variable  xp  =  p/pmax,  as  denned  for  A+  in  the  last  chapter.  eP  physically 
represents  the  "hardness"  of  the  heavy  quark  fragmentation.  The  smaller  the  value 
of  eP  is,  the  larger  the  portion  of  the  total  momentum  the  heavy  quark  carries. 

Although  the  Peterson  function  agrees  pretty  well  with  the  fragmentation  of  D° 
and  D*+  measured  in  the  experiments  [34]  [35],  there  is  no  reason  why  any  charmed 
particle  momentum  spectrum  should  follow  an  analytical  curve  of  this  type,  as  many 
particles  are  the  decay  results  of  excited  states.  However,  the  measurement  of  eP  of 
Peterson  function  remains  useful  for  comparing  the  spectra  of  different  particles. 

To  measure  the  Peterson  function  of  a  particle,  we  fit  the  distribution  of  efficiency- 
corrected  cross  section  at  a  different  xp  range  and  extract  the  parameter  eP.  Table 
4.8  shows  the  yields  and  Monte  Carlo  derived  efficiencies  at  different  xp  ranges.  Also 
listed  are  the  efficiency-corrected  yields  normalized  to  xp  >  0.5. 

Figure  4.14  shows  the  ^  (N  is  the  efficiency-corrected  number  of  A&'s  for 
xp  >  0.5)  distribution  for  the  A+^2593)  and  the  A+(2625)  based  on  the  efficiency- 
corrected  yields  listed  in  Table  4.8.  Only  statistical  errors  are  shown  in  Table  4.8. 
Systematic  errors  are  dominated  by  fitting  procedures  and  are  found  to  be  a  little 
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Figure  4.14:  Normalized  spectra  of  scaled  momentum  for  A+x  samples.  The  data  for 
the  A+t(2593)  (diamonds)  and  the  A+^2625)  (triangles)  are  fitted  to  the  Peterson 
function  overlaid  as  the  solid  and  dashed  curves,  respectively. 

smaller  than  the  statistical  errors.  The  error  bars  shown  in  Figure  4.14  are  the  total 
errors.  Also  shown  in  the  figure  are  the  curves  obtained  by  fitting  the  distributions 
by  Peterson  function.  The  measurement  results  are  listed  in  Table  4.9  compared 
to  the  previous  measurements.  The  measurements  show  that  A^  fragmentation 
at  10  GeV/c2  is  close  to  the  fragmentation  of  D*+,  which  was  measured  to  have 
eP  =  0.078  ±0.008(36]. 
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Table  4.9:  Measurements  of  the  fragmentation  of  A^'s  using  Peterson  functions. 


Experiments 

ep 

Ac+1(2593) 

CLEO-II[30] 
ARGUS[32] 
CLEO-II(this) 

0.057±0.023±0.016 
0.069±g-git  ±  0.040 
0.064l8fJ 

Ac+!(2625) 

ARGUS[28] 
CLEO-II[30] 
CLEO-II  (this) 

0.044±0.018 
0.065±0.016±0.013 
0.1081q;oo9 

4.8    Decay  Ratios  and  Production  Ratios 


4.8.1    Systematic  Errors  on  Yields  and  Efficiencies 

To  calculate  the  branching  ratios  and  productions,  the  yields  on  the  signals  and 
the  Monte  Carlo  derived  efficiencies  are  used.  Both  bring  systematic  errors  to  the 
result. 

The  systematic  errors  on  the  yield  of  a  signal  come  from  the  fitting  proce- 
dures, including  the  signal  and  background  parameterizations  and  the  mass  res- 
olution determined  by  the  Monte  Carlo.  As  we  explained  in  the  last  section,  for  the 
£*(2593)->>A+7r+7r~  signal,  the  systematic  error  is  dominated  by  the  uncertainty  of 
the  signal  shape,  and  background  parameterization  is  also  a  factor  since  the  signal 
is  just  at  the  kinematic  threshold.  Based  on  Table  4.1  and  other  fit  results,  the 
fitting  errors  of  the  £«(2593)-fA+7r+7r  yield  are  ±15%  from  the  peak  parameteri- 
zation and  ±8%  from  the  background  parameterization.  As  we  tested,  varying  the 
Gaussian  resolution  by  a  reasonable  amount  only  changes  the  A+^2593)  yield  by 
2%.  We  therefore  estimate  the  fitting  error  on  the  A+1(2593)^A+7r+7r-  yield  to  be 
±17%.  For  A+1(2625)->A+7r+7r-,  we  have  a  clean  signal  at  very  low  background, 
so  the  systematic  error  comes  mostly  from  the  Monte  Carlo  predicted  resolution  or 
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equivalently,  the  uncertain  natural  width  of  this  state.  We  finally  estimate  that  the 
fitting  error  of  the  A+1(2625)->A+7r+7r-  yield  is  t\°%. 

As  the  A^— ^A+7T07r°  signals  are  poor,  both  statistical  and  fitting  errors  are  very 
large.  The  fitting  errors  are  mostly  due  to  the  uncertainty  of  the  mass  resolutions 
and  the  natural  widths  of  the  states.  We  expect  the  A+1(2593)->A+7r°7r°  signal  to  be 
nearly  a  full  Breit-Wigner  convoluted  with  a  Gaussian,  since  the  kinematic  thresh- 
old is  at  301  MeV/c2.  The  fitting  errors  are  evaluated  by  reasonably  varying  the 
background  parameterizations,  the  mass  resolutions,  and  the  widths  of  the  states. 
From  the  results  listed  in  Table  4.3  and  some  further  tests,  we  estimate  the  fitting 
errors  on  A+^A+tt0^0  signals  are  ±13%  for  A+^2593)  and  tf  %. 

Since  we  measure  the  branching  fractions  by  ratios,  the  systematic  errors  on  the 
efficiencies  of  A+  reconstruction  will  cancel.  Based  on  comprehensive  studies  by 
the  CLEO  collaboration,  the  systematic  errors  on  efficiencies  are  1%  on  finding  one 
charged  track  and  5%  on  finding  a  7r°. 

4.8.2    A+ttV,  A+tt0,  and  A+7  Decays 

Since  the  yield  of  A^-^A+tt+tt-  was  the  best  measured  for  AJ,  the  branching 
ratios  we  measured  are  all  given  relative  to  this  decay  mode. 

This  is  the  first  time  that  the  decay  A+j-^A+ttV  is  reported.  Of  interest  is 
the  branching  ratio  #(A+7r07r°)/#(A+7r+7r-).  If  only  isospin  symmetry  is  applied  on 
these  states,  this  ratio  should  be  1/2. 

Based  on  the  discussion  in  the  last  subsection,  we  assign  the  systematic  errors  on 
S(Ac+7r°7r0)/fi(Ac+7r+7r-)  to  be  ±26%  for  the  A+(2593)  and  ±™%  for  the  A+(2625). 
Table  4.10  gives  the  result  calculated  from  the  yield  and  the  Monte  Carlo  predicted 
efficiency  of  each  decay  mode. 
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Table  4.10:  Branching  ratios  of  A+7r°7r°  relative  to  A+7r+7r  . 


Particle 

5(A+7r°7r°)/B(A+7r+7r-) 

(2593) 
A+(2625) 

2.6±0.6±0.7 
1.0±0.3±S;| 

Table  4.11:  Limits  on  branching  ratios  of  A+x— >-A+7r0  and  A^— >-A+7  relative  to 
A^-^A+Tr+Tr-.  All  values  are  evaluated  at  90%  C.L. 


Experiments 

5(A+7r°)/5(A+7r+7r-) 

£(A+7)/£(A+7r+7r-) 

(2593) 
A+(2625) 

CLEO-II[30] 
CLEO-II(this) 

CLEO-II[30] 
CLEO-II(this) 

<  3.53 

<  1.1 

<  0.91 

<  0.45 

<  0.98 

<  3.0 

<  0.52 

<  0.36 

The  result  on  the  A^(2625)  is  larger  than  0.5  but  consistent  with  it,  but  the 
number  on  the  A^(2593)  is  significantly  larger  than  0.5.  This  can  be  explained  by 
the  phase-space  suppression.  The  A^(2593)  mostly  decays  into  Ec7r  first  then  goes 
to  the  final  state  A+7T7T.  Since  the  mass  of  the  A^(2593)  is  just  at  the  threshold 
to  decay  into  a  Ec  plus  a  charged  pion,  so  the  decay  A+x(2593)  — >■  £++7r~  and 
A^(2593)  — >  E^7r+  will  be  suppressed.  On  the  other  hand,  isospin  symmetry  is  not 
exact,  so  the  n°  mass  is  about  5  MeV/c2  lower  than  charged  ones,  so  the  decay 
A^(2593)  ->  £+7T°  is  not  suppressed  by  phase  space.  Thus  for  the  A^(2593),  we 
expect  B(A^  -4  A+7r°7r°)/5(A+1  ->  A+tt+tt")  to  be  higher  than  1/2.  We  will  discuss 
this  more  in  the  next  chapter. 

According  to  Table  4.5,  taking  the  systematic  errors  on  A+7r+7r~  yields  into 
account,  we  estimated  the  upper  limits  of  A+7r°  and  A+7  decay  ratios  relative  to 

A+7T+7T". 
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Table  4.12:  Branching  fractions  of  substructure  E+  in  A+j— »A+7r+7r  decays.  Note 
the  A^i  yields  are  smaller  than  the  fitted  results  in  Figure  4.1  since  we  measure 
the  Ec  yields  from  304  MeV/c2  to  314  MeV/c2  of  the  A+^2593)  peak  and  from  337 
MeV/c2  to  347  MeV/c2  of  the  A+j(2625)  peak. 


Decay 

A+  Yield 

Ec  Yield 

A+(2593) 

E++7T- 

E°7T+ 
E+7T0 

137  ±22  ±22 
137  ±  22  ±  22 
93  ±  20  ±  22 

67  ±  7  ±  5 
65  ±  7  ±  5 
74  ±11  ±10 

0.45  ±0.12  ±0.11 
0.47±  0.12  ±0.11 
0.73  ±0.20  ±0.20 

A+x(2625) 

S++7T" 
E°7T+ 

450  ±  29  ±  9 
450  ±  29  ±  9 

24  ±  10  ±  7 
32  ±  10  ±  7 

0.05  ±0.02  ±0.02 
0.07  ±0.02  ±0.02 

Table  4.13:  Comparison  of  measurements  of  branching  ratios  of  E++  and  E°  sub- 
structures in  A^(2593)^A^7r+7r"  decays  from  different  experiments. 


Experiment 

/ej> 

E687[29] 
CLEO-II[30] 
ARGUS[32] 
CLEO-II(this) 

0.36  ±  0.09  ±  0.09 
0.37  ±0.12  ±0.13 
0.45  ±0.12  ±0.11 

0.42  ±0.09  ±0.09 
0.29  ±0.10  ±0.11 
0.47  ±0.12  ±0.11 

>  0.51  90%  C.L. 
0.66±8;ii  ±  0.07 

4.8.3    Ec7r  Decays 

The  A(1405),  the  spin  1/2  strange  baryon,  decays  100%  to  E7r;  and  the  spin 
3/2  state  A(1520)  decays  to  both  E7r  and  Ann.  We  have  observed  the  similar 
phenomenon  for  the  A;!;  (2593)  and  the  A+1(2625).  Measurements  of  fEc  =  B{A^  -> 
Ectt)/B(A+1  ->  A+tttt)  for  E++,  E+,  and  E°  are  tabulated  in  Table  4.12.  The 
systematic  errors  are  dominated  by  the  fitting  errors  on  Ec  yields  on  A£  signals  and 
backgrounds. 

Our  observation  is  close  to  the  previous  CLEO  measurements  and  consistent 
with  the  prediction  that  the  A+x(2593)  and  A+^2625)  have  Jp  1/2"  and  3/2",  re- 
spectively. The  A+j(2593)  can  decay  through  both  S- wave  Ecn  and  P-wave  non- 
resonance  A+7r+7r-  channels,  but  the  5-wave  is  preferred.  However,  angular  mo- 
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Table  4.14:  Comparison  of  measurements  of  branching  ratios  of  E^"+  and  £°  sub- 
structures in  A^(2625)— >-A+7r+7r~  decays  from  different  experiments. 


Experiment 

/s++ 

ARGUS[28] 

E687[31] 
CLEO-II[30] 
CLEO-II(this) 

<  0.08  90%  C.L. 
0.05  ±0.02  ±0.02 

<  0.07  90%  C.L. 
0.07  ±0.02  ±0.02 

0.46  ±0.14 
<  0.36  90%  C.L. 

Table  4.15:  Measurements  of  A^— ^A+7r+7r  production  ratio  relative  to  the  total 
inclusive  A+  production. 


Experiment 

A+(2593) 

CLEO[30] 
ARGUS  [32] 
CLEO(this) 

1.44  ±0.24  ±0.30 
2.1±f;$  ±1.1 

1  /.fi+0.23+0.29 
1-4O-0.22-0.29 

A+(2625) 

ARGUS[28] 
CLEO[30] 
CLEO(this) 

4.1  ±  1.0  ±0.8 

3.51  ±0.34  ±0.28 
3.47  ±0.1718$ 

mentum  and  parity  conservations  require  Al1(2625)— >T,cir  to  be  jD-wave,  so  for  the 
A^(2625),  the  P-wave  non-resonance  Al7r+7r~  is  preferred  and  expected  to  have  a 
smaller  width. 


4.8.4         Production  Ratios 

Since  we  do  not  have  good  statistics  with  Al7r°7r°  decay  channels,  we  only  present 
the  A+^A+n+n'  production  as  a  fraction  of  A+  baryons  in  the  10  GeV/c2  e+e" 
annihilations.  We  have  only  measured  the  number  of  A+'s  and  A^'s  at  xp(A^)  >  0.5, 
so  we  need  to  extrapolate  the  Peterson  function  to  zero  momentum.  Fortunately, 
the  two  A^'s  are  fragmented  very  hard  (only  about  20%  are  fragmented  below 
xp  =  0.5,)  so  the  systematic  errors  from  extrapolation  are  not  very  large.  Based  on 
our  estimation,  the  sources  and  their  values  of  systematic  errors  of  the  A^(2593) 
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are  15%  for  the  peak  parameterization,  8%  for  the  background  parameterization, 
and  10%  for  the  xp  =  0  extrapolation.  We  estimate  the  total  systematic  error  of 
the  A^(2593)  production  ratio  measurement  to  be  20%.  For  the  A+1(2625),  the 
fitting  errors  are  dominated  by  the  uncertainty  of  the  detector  resolution.  The  total 
fitting  error  is  estimated  to  be  t\0%,  and  the  xp  =  0  extrapolation  uncertainty  is 
estimated  to  be  7%.  Thus  the  total  systematic  error  of  the  A^(2625)  production 
ratio  measurement  is  estimated  to  be  lg2%-  Our  results  are  consistent  with  the 
previous  experiments. 


CHAPTER  5 
STUDIES  OF  Ec  AND  E*  BARYONS 

5.1  Introduction 

The  isospin  triplet  (E++,  £+,  £[!)  is  the  charm  counterpart  of  the  E  triplet  with 
the  two  light  quarks  antisymmetric  under  exchange.  The  quark-model-predicted 
spin-parity  Jp  for  the  triplet  is  l/2+  with  the  two  light  quarks  in  a  Jp  =  1+  con- 
figuration. Although  all  three  states  have  been  observed  by  many  experiments[37, 
38,  39,  40,  41,  42,  43,  44],  only  the  mass  splittings  and  production  ratios  of  the  E++ 
and  the  E°  are  reported,  and  no  attempt  has  been  made  to  measure  their  intrinsic 
widths.  With  a  much  larger  data  sample  and  the  improved  resolution  of  the  CLEO- 
II  detector,  we  are  now  able  to  measure  the  widths  and  production  ratios  of  all  three 
states. 

The  3/2+  excitation  of  the  Ec  is  the  E*  triplet.  CLEO  has  recently  made  the 
first  observation  of  the  E*++  and  the  E*°[21].  Here  we  present  an  update  of  the 
measurements  of  these  two  particles  and  give  evidence  of  the  observation  of  the 
previously  undiscovered  E*+. 

In  this  chapter  we  first  describe  the  observations  and  the  fits  of  Ec  and  E* 
signals  in  section  2.2,  then  in  later  sections  we  present  the  results  of  masses,  widths, 
fragmentation  functions,  and  production  ratios  of  these  particles.  The  physical 
significance  of  these  measurements  is  also  discussed.  Finally,  we  describe  a  tentative 
measurement  on  the  mass  and  width  of  the  A+j(2593)  in  a  new  way,  different  from 
the  methods  used  in  previous  experiments.  We  present  this  last  because  it  requires 
the  knowledge  of  the  Ec  mass  and  width. 
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5.2  Signals 

All  Ec's  and  E*'s  have  one  decay  mode:  A+7T.  To  obtain  the  signals,  we  combine 
each  A+  candidate  with  a  pion,  then  measure  the  mass  difference  M(A+n)  —  M(A+). 
The  A+  candidates  are  selected  in  the  same  way  as  we  did  in  observing  A^'s.  To 
improve  the  signal  to  background  ratio,  a  cut  of  xp  >  0.5  for  A+n*  and  a  harder 
cut  of  xp  >  0.6  for  A+7T0  are  used,  since  the  A+7r°  combination  has  much  larger 
background  and  is  more  difficult  to  observe.  Here  xp  is  defined  in  the  same  way  as 
the  xp  of  the  A+x. 

All  charged  tracks  are  used  to  be  the  independent  candidates.  In  order  to 
suppress  the  large  background  from  these  low-momentum  charged  pions,  a  cut  on 
cos(9dec)  is  used  for  A+n*  modes.  9dec  is  defined  to  be  the  angle  between  the  the 
pion  momentum  measured  in  the  rest  frame  of  the  A+tt*,  and  the  direction  of  the 
A+7r±  in  the  laboratory  frame.  According  to  Monte  Carlo  studies,  9dec  cut  is  more 
efficient  than  simply  cutting  on  7r±  momentum  to  improve  the  signal  to  background 
ratio  of  the  E*++  and  the  E*°.  We  used  cos(9dec)  >  -0.4  which  is  the  optimization 
from  the  E*-»Ec7rMonte  Carlo  studies. 

When  reconstructing  A+tt*  invariant  mass,  we  constrain  the  A+  and  the  inde- 
pendent pion  to  be  from  the  event  vertex,  whose  finding  procedure  was  described  in 
Chapter  3.  This  improves  the  mass  resolution  by  roughly  5%  for  both  the  E++  and 
the  E°. 

7T°  selection  is  described  in  Chapter  3.  To  suppress  the  large  background,  a  7r° 
momentum  cut  is  needed.  The  Monte  Carlo  shows  that  a  direct  tt0  momentum  cut 
is  slightly  more  efficient  than  9dec  cut  for  the  E*+  observation  and  predicted  that  a 
150  MeV/c  cut  gives  largest  signal  to  background  ratio. 

Sufficient  Monte  Carlo  events  have  been  generated.  The  mass  differences  we  used 
are  167.87  MeV/c2  for  the  E++,  168.70  MeV/c2  for  the  E+,  167.30  MeV/c2  for  the 
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E°,  and  233  MeV/c2  for  all  three  E*'s.  We  found  that  slightly  varying  the  generated 
mass  does  not  affect  the  Monte  Carlo  prediction  for  either  width  or  efficiency.  We 
first  generated  events  with  zero  width  to  get  the  resolution  of  the  signal;  then  after 
we  extracted  the  width  from  real  data,  we  generated  events  with  the  found  width. 
Efficiencies  are  calculated  using  the  latter  Monte  Carlo  events.  Using  the  same 
signal-fitting  functions,  the  same  natural  width  can  be  extracted  from  the  Monte 
Carlo  signal.  All  reconstructed  mass  differences  in  Monte  Carlo  shift  by  less  than 
0.03  MeV/c2  relative  to  the  generated  ones. 

5.2.1  A+tt* 

Figure  5.1  shows  the  mass-difference  spectrum  of  M(A+7r±)  -  M(A+).  The  clear 
peaks  near  167  MeV/c2  are  Ec's,  and  wide  bumps  near  233  MeV/c2  are  E*'s.  In  addi- 
tion, we  should  expect  broad  enhancements  below  205  MeV/c2  which  are  A+7r  com- 
binations from  A+1(2625)^A+7r+7r_  decays.  It  is  because  of  this  background  that 
we  did  the  A^  analysis  before  the  E^'s.  It  is  calculated  from  A+1(2625)->A+7r+7r- 
signal  and  Monte  Carlo  predicted  efficiencies  that  the  excess  from  this  source  yields 
728  events  between  150  MeV/c2  and  205  MeV/c2. 

In  the  same  plots  shown  as  solid  histograms  are  the  mass-difference  spectra  using 
the  sidebands  of  the  A+signal.  The  sidebands  are  taken  from  the  A+  combinations 
in  the  mass-difference  region  3  -  4.6  aM  away  from  the  known  A+  mass,  2286.7 
MeV/c2.  No  enhancement  can  be  observed  from  these  histograms,  and  they  can  be 
very  well  fitted  by  second-order  polynomials. 

In  Figure  5.1,  Ec  and  E*  signals  are  fitted  with  Breit-Wigners  convoluted  with 
a  Gaussian  with  a  fixed  resolution  width  of  1.04  MeV/c2  for  the  Ec  or  1.53  MeV/c2 
for  the  E*.  The  Gaussian  resolutions  are  determined  by  the  Monte  Carlo  studies. 
The  broad  excesses  in  the  region  below  205  MeV/c2  are  fitted  by  a  double-threshold 
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Figure  5.1:  Mass  difference  for  (a)  M(A+tt+)  -M(A+),  and  (b)  M(A+tt-)  -M(A+). 
The  solid  histogams  are  the  spectra  of  M(A+7r±)  -M(A+  "side"),  where  a  "A+  side" 
is  from  Aj"  sidebands. 
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function  with  a  fixed  shape  determined  from  A+7r  shapes  in  A+1(2625)— >A+7r+7r~ 
Monte  Carlo  (with  thesholds  at  139  MeV/c2  and  205  MeV/c2),  and  a  fixed  area, 
728,  determined  by  the  efficiency-corrected  A+1(2625)  production.  The  rest  of  the 
background  is  fitted  by  a  second-order  polynomial  with  all  parameters  floated. 

The  results  are  sensitive  to  different  fitting  procedures.  The  extracted  parameters 
from  some  variations  of  fits  are  listed  in  Table  5.1.  Fit  Type  1,  which  has  the 
fitted  curves  shown  in  Figure  5.1,  provides  the  presented  results.  Type  2  fixes  the 
convoluted  Gaussian  width  of  the  Ec's  to  be  1.25  MeV/c2,  widening  the  Monte 
Carlo  predicted  resolution  of  the  Ec's  by  20%.  Type  3  reduces  the  area  of  the 
double-threshold  function  by  50%,  and  Type  4  enlarges  it  by  50%.  Type  5  uses  a 
third  order  of  polynomial  instead  of  the  2nd  order  to  fit  the  background.  Type  6 
fits  the  region  in  the  plot  below  200  MeV/cc  and  Type  7  fits  the  region  above  200 
MeV/c2.  Type  8  uses  4  MeV/c2  bin  size  instead  of  1  MeV/c2.  The  statistical  errors 
given  in  the  table  are  only  for  fit  Type  1  since  the  errors  for  the  other  types  are 
roughly  the  same. 

Besides  the  fit  types  listed  in  Table  5.1,  many  other  fitting  procedures  are  used, 
such  as  varying  and  floating  the  parameters  of  the  double-threshold  function  (we 
found  that  our  results  are  much  more  sensitive  to  the  area  than  the  actual  shape 
of  the  double-threshold  function),  using  higher-order  Chebychev  polynomials  to  fit 
the  entire  background,  using  a  single  Gaussian,  a  double-Gaussian,  or  a  bifurcated 
Gaussian  to  fit  the  peaks  to  check  the  variation  of  the  means,  fitting  the  two  peaks 
separately  in  different  mass-difference  ranges  and  with  different  backgrounds.  The 
measured  parameters  from  all  reasonable  fits  are  within  the  range  listed  in  Table 
5.1.  Also,  we  measured  the  same  parameters  with  different  oM  cuts  around  the  A+ 
mass  and  found  that  the  signals  we  have  are  consistently  associated  with  the  A+. 

To  study  the  fragmentation  function  of  these  particles,  we  divide  the  0.5  -  1.0  xp 
range  into  4  bins  and  add  an  extra  bin  from  0.4  to  0.5,  then  fit  each  mass-difference 
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Table  5.1:  The  yield  (Y),  mass  difference  AM,  and  extracted  width  T  measured 
from  various  fits  to  A+7r±. 


Fit  Type 

y 

AM(MeV/c2) 

r(MeV/c2) 

1 

880±™ 

167.34  ±0.10 

2.1  ±0.3 

2 

845 

167.35 

1.7 

3 

915 

167.35 

2.3 

4 

851 

167.33 

2.1 

5 

875 

167.34 

2.1 

6 

908 

167.36 

2.2 

1 

895t& 

167.24  ±0.11 

o  q+0.4 
^•°-0.3 

2 

870 

167.23 

1.9 

3 

942 

167.24 

2.1 

4 

925 

167.22 

2.0 

5 

879 

167.22 

2.3 

c 
0 

oyo 

10/  .2,1 

2.3 

1 

231.4  ±0.8 

14  ±3 

2 

825 

231.4 

15 

3 

649 

231.6 

12 

4 

931 

231.2 

16 

5 

829 

231.4 

15 

7 

715 

231.4 

13 

8 

821 

231.7 

15 

Ef 

1 

10151^ 

231.9  ±  1.2 

203 

2 

769 

232.8 

16 

3 

811 

233.1 

16 

4 

1249 

231.6 

24 

5 

822 

232.7 

17 

7 

1193 

231.9 

23 

8 

1116 

231.7 

22 

100 


Table  5.2:  The  yield  (Y),  efficiency,  and  normalized  production  of  A+tt*  events  for 
different  xv  cuts. 


Xp 

Y 

Efficiency  (%) 

Normalized  Production 

0.4  -  1.0 

1094  ±  54 

6.64 

16988  ±  813 

0.4  -  0.5 

211  ±29 

6.57 

3212  ±  441 

0.5  -  0.6 

246  ±  27 

6.61 

3722  ±  408 

0.6  -  0.7 

295  ±  26 

6.69 

4410  ±  389 

0.7-0.8 

232  ±  23 

6.64 

3493  ±  364 

0.8  -  1.0 

137  ±  16 

6.42 

2134  ±249 

0.4  -  1.0 

1109  ±55 

6.64 

16702  ±  828 

0.4  -  0.5 

214  ±28 

6.57 

3257  ±426 

n  ^   n  fi 

v.o  u.u 

941  -1-  97 

O.Ol 

o040  ±  olo 

0.6  -  0.7 

289  ±  27 

6.69 

4319  ±  403 

0.7-0.8 

209  ±  23 

6.64 

3148  ±  346 

0.8  -  1.0 

120  ±  14 

6.42 

1869  ±218 

0.4  -  1.0 

1097 ± 121 

7.06 

15538  ±  1714 

0.4  -  0.5 

244  ±  67 

7.07 

3451  ±  948 

0.5  -  0.6 

321  ±  55 

7.23 

4440  ±  761 

0.6  -  0.7 

282  ±  46 

7.09 

3977  ±  649 

0.7-0.8 

155  ±  33 

6.70 

2313  ±  493 

0.8  -  1.0 

70  ±  15 

6.88 

1017  ±218 

0.4-  1.0 

1348  ±  189 

7.06 

19093  ±  2677 

0.4  -  0.5 

311  ±79 

7.07 

4399 ±  1117 

0.5  -  0.6 

384  ±  71 

7.23 

5311  ±  982 

0.6  -  0.7 

277  ±  60 

7.09 

3907  ±  846 

0.7  -  0.8 

227  ±47 

6.70 

3388  ±  701 

0.8  -  1.0 

134  ±  34 

6.88 

1948  ±  494 
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spectrum  individually.  Each  Gaussian  width  is  fixed  by  the  predicted  resolution 
from  corresponding  Monte  Carlo  events.  If  we  float  the  masses  and  the  Breit- Wigner 
widths  of  all  signals,  both  mass  and  width  of  the  same  particle  vary  considerably 
due  to  poor  statistics.  However,  if  we  fix  all  the  masses  to  the  measured  values 
in  the  0.5  -  1.0  xp  range,  all  widths  vary  around  our  presented  values,  measured 
in  the  0.5  -  1.0  xp  range,  within  statistical  errors.  If  we  fix  all  the  widths  to  be 
our  presented  values,  Ec  masses  vary  within  0.16  MeV  around  our  presented  masses 
with  statistical  errors  about  ±0.20  MeV/c2,  and  E*  masses  vary  within  tlx  MeV/c2 
around  our  presented  masses  with  statistical  errors  of  about  ±2  MeV/c2.  Table  5.2 
listed  the  yields  and  efficiency-corrected  productions  with  all  widths  fixed.  Each 
efficiency  is  determined  by  the  Monte  Carlo,  and  its  error  is  within  0.01%. 

We  also  fitted  the  mass-difference  spectra  with  variations  of  the  9dec  cut  and, 
almost  equivalently,  of  the  n4  momentum  cuts  to  see  if  there  are  any  systematic 
changes  in  masses  and  widths.  We  find  the  signals  of  the  data  are  consistent  with 
Monte  Carlo  studies,  and  all  extracted  parameters  agree  within  statistical  errors. 

5.2.2  A+tt° 

For  A+7T0  combinations,  7r°  momentum  cut  is  very  important  to  reduce  the  back- 
ground. The  150  MeV/c2  cut  we  used  is  the  prediction  from  the  Monte  Carlo.  Table 
5.3  displays  signal-to-background  ratios  of  A+7r°  combinations  for  different  7r°  mo- 
mentum cuts.  In  the  table  S"s  are  the  fitted  areas  of  the  Ec  or  the  E*  from  the 
Monte  Carlo  plots;  £Ts  are  the  size  of  the  background  from  the  Ec  or  E*  sidebands 
from  A+7r°  combinations  of  real  data. 

The  plot  of  mass  difference  M(A+tt°)  -  M(A+)  is  shown  in  Figure  5.2.  The 
two  signals  and  the  background  are  fitted  by  the  same  four  functions  used  in  Fig- 
ure 5.1.   The  convoluted  Gaussian  widths  are  determined  by  Monte  Carlo  stud- 
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Table  5.3:  Signal  versus  background  of  A+7r°  combinations  for  different  7r°  momen- 
tum cuts.  In  the  table  S  is  the  fitted  area  of  the  Ec  or  the  E*  signal  from  the  Monte 
Carlo  events;  B  is  the  area  of  the  background  (taken  from  the  Ec  or  E*  sidebands) 
from  A+7T0  combinations  of  real  data. 


P(7T°)  (GeV) 

5E+ 

0.05 

2877 

2610 

5448 

9977 

39.0 

26.1 

0.10 

2873 

2526 

5376 

7760 

39.3 

28.7 

0.15 

2734 

2398 

4685 

6068 

40.0 

30.8 

0.20 

2363 

2010 

3676 

4685 

39.0 

30.7 

0.25 

1929 

1694 

2632 

3654 

37.6 

28.0 

ies:  1.92  MeV/c2  for  the  E+  at  167  MeV/c2  and  3.71  MeV/c2  for  the  E*+  at  229 
MeV/c2.  The  shape  of  the  double-threshold  function  is  determined  by  A+7r°  shape 
in  the  A+1(2625)->A+7r07r0  Monte  Carlo,  and  the  two  thresholds  are  fixed  to  be  at 
135  MeV/c2  and  216  MeV/c2.  The  area  of  the  double-threshold  function  is  the 
number  of  A+1(2625)->A+7r°7r°  feed-down,  and  it  is  calculated  (from  the  efficiency- 
corrected  A+1(2625)->A+7r°7r°)  to  be  216.  The  solid  histogram  in  the  plot  shows  the 
spectrum  of  the  combination  of  7r°'s  and  the  A+  sidebands,  taken  from  the  same 
mass-difference  region  (3  -  4.6  aM)  as  in  the  A+ir*  case. 

Fitting  procedures  affect  the  extracted  parameters.  Many  different  fits  were 
performed  to  investigate  the  stabilities  of  the  signals  using  different  fitting  methods. 
Table  5.4  shows  some  results  from  different  fit  types.  Fit  Type  1  is  the  procedure 
just  described  whose  fitted  curve  is  shown  in  Figure  5.2,  and  the  fitted  results  are 
our  presented  values.  Type  2  fixes  the  width  of  the  convoluted  Gaussian  of  the  E+ 
peak  to  be  2.20  MeV/cc,  15%  larger  than  the  Monte  Carlo  predicted  resolution;  and 
Type  3  fixes  it  to  be  2.50  MeV/c2,  30%  larger  than  the  Monte  Carlo  value.  We  found 
that  the  fitted  T  of  the  E*+  is  insensitive  to  the  change  of  the  convoluted  Gaussian 
width.  In  fit  Type  4  the  double-threshold  function  representing  the  feed-down  of 
A^i(2625)  is  taken  away;  and  in  Type  5  its  area  is  enlarged  by  100%.  Type  6  uses  a 
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Figure  5.2:  Mass  difference  M(A+tt°)  -  M(A+). 

Table  5.4:  The  yield  Y,  mass  difference  AM,  and  extracted  width  T  measured  from 
various  fits  to  the  A+7r°  signal. 


Fit  Type 

Y 

AM(MeV/c2) 

r(MeV/c2) 

1 

367±^ 

166.9  ±0.5 

4.5±& 

2 

356 

166.9 

3.9 

3 

343 

167.0 

3.2 

4 

415 

167.0 

5.4 

5 

329 

166.8 

3.9 

6 

337 

166.8 

4.0 

7 

376 

166.9 

4.7 

9 

352 

166.9 

3.7 

10 

355 

167.1 

4.5 

1 

236+^9 

229.2  ±  1.3 

r  o+().4 
°-°-4.6 

2 

224 

229.2 

5.2 

3 

231 

229.2 

5.1 

4 

190 

229.2 

3.3 

5 

297 

229.0 

7.7 

6 

203 

229.4 

3.7 

8 

215 

229.5 

4.2 

9 

238 

229.3 

5.5 

10 

217 

229.4 

4.1 
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third-order  polynomial  instead  of  the  second  order  in  Type  1  to  fit  the  background. 
Fit  Type  7  and  8  try  to  eliminate  the  influence  of  one  signal  to  the  other.  Type  7 
excludes  the  region  from  200  MeV/c2  to  250  MeV/c2  so  only  the  E+  signal  is  fitted. 
Type  8  fits  the  region  above  200  MeV/c2  so  the  E+  signal  is  excluded.  All  the  above 
fit  types  fit  the  spectrum  with  2  MeV/c2  bin  size.  Types  9  and  10  are  the  fits  to 
plots  with  the  bin  size  1  MeV/c2  and  4  MeV/c2,  respectively.  The  statistical  errors 
for  all  fit  types  are  roughly  the  same  as  the  ones  given  in  Type  1. 

As  in  the  A+7r±  cases,  we  have  used  more  variations  of  fits,  and  the  variations 
on  all  measured  parameters  do  not  exceed  the  ranges  of  those  listed  values  in  Table 
5.4. 

It  is  worth  mentioning  that  although  the  yield  of  the  E*+,  which  is  236l7°9, 
does  not  seem  statistically  significant,  if  we  fit  the  signal  by  a  single  Gaussian  with 
floated  width,  a  162±™  signal  is  obtained  with  a  mean  value  of  229.4  ±1.4  MeV/c2. 
We  also  notice  that  if  we  fit  the  mass-difference  region  from  190  MeV/c2  to  290 
MeV/c2  (E+  signal  only)  in  Figure  5.2,  the  x2  obtained  in  our  best  fit  is  31.5  for 
44  d.o.f.  (degrees  of  freedom),  whereas  in  a  fit  with  no  signal  allowed,  the  \2  rises 
to  49.5.  The  fit  strongly  supports  the  existence  of  a  signal,  and  this  represents  the 
first  evidence  for  the  particle  E*+. 

When  we  vary  7r°  momentum  cuts  from  50  MeV/c  to  250  MeV/c,  we  found  that 
the  mass  difference  shifts  from  166.1  to  166.9  MeV/c2  for  the  E+  and  from  233.5  to 
229.7  MeV/c  for  the  E*+;  and  the  extracted  width  varies  from  3.4  to  4.8  MeV/c2 
for  the  Ec  and  from  5.0  to  8.4  MeV/c2.  Also,  for  the  yield  in  the  E+  signal,  the 
data  agrees  fairly  well  with  the  Monte  Carlo.  Although  we  have  obtained  the  best 
statistical  significance  for  the  E*+  signal  at  200  MeV/c7r°  momentum  cut,  the  Monte 
Carlo  optimized  150  MeV/c  cut  is  actually  used.  To  prove  the  consistency  of  our 
signal  with  A+  samples,  we  looked  at  the  signal  with  different  oM  cuts  around  the 
A+  mass,  and  the  measured  masses,  widths,  and  yields  agree  very  well. 
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Table  5.5:  The  yield  (Y),  mass  difference  (AM),  efficiency,  and  normalized  produc- 
tion of  the  E+  for  different  xp  cuts.  Both  M  and  T  are  in  units  of  MeV/c2. 


Xp 

AM 

Y 

Efficiency(%) 

} '/Efficiency 

0.5  -  1.0 

167.1  ±0.4 

559  ±  60 

2.77 

20181  ±2166 

0.5  -  0.6 

167.4  ±0.7 

208  ±  38 

2.45 

8490  ±  1551 

0.6  -  0.7 

167.2  ±0.8 

168  ±  34 

2.82 

5957  ±  1206 

0.7-0.8 

166.4  ±0.7 

182  ±  28 

2.88 

6319  ±  972 

0.8  -  1.0 

168.2  ±  1.3 

72  ±  18 

2.85 

2526  ±  632 

We  tabulated  the  masses  and  the  yields  of  the  E+  for  different  xp  ranges  in 
Table  5.5.  For  the  E*  signal,  lacking  statistics  makes  it  very  difficult  and  not  quite 
meaningful  to  fit  the  mass-difference  spectra  at  different  xp  bins.  The  numbers  in  the 
table  are  all  obtained  by  fixing  the  Breit- Wigner  width  to  be  4.5  MeV/c2  and  fixing 
the  width  of  the  convoluted  Gaussian  to  be  the  number  obtained  by  Monte  Carlo, 
ranging  from  1.80  MeV/c2  to  2.32  MeV/c2.  All  statistical  errors  on  efficiencies  are 
no  larger  than  0.01%. 

5.3  Masses 

The  masses  we  present  are  the  mean  values  of  the  Breit- Wigner  convoluted  with 
a  Gaussian.  We  have  checked  that  fitting  with  a  single  Gaussian,  a  double  Gaussian, 
or  a  bifurcated  Gaussian  (constraining  the  two  width  difference  no  larger  than  15%) 
gives  the  mean  value  within  1  a  of  the  statistical  error. 

As  discussed  in  the  last  chapter,  the  uncertainties  in  the  magnetic  field  nor- 
malization and  the  energy-loss  correction  applied  to  the  charged  soft  pions,  or  the 
uncertainty  of  the  shower  energy  calibration  applied  to  the  neutral  pions,  bring  sys- 
tematic errors  on  mass  measurement  for  M(A+7r)  -  M(A+).  We  find  that  0.2% 
change  of  charged  track  momentum  shifts  M(£c)  -  M(A+)  by  about  0.10  MeV/c2 
and  M(E*)  -  M(A+)  by  about  0.30  MeV/c2,  1.0%  change  of  photon  energy  shifts 
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Table  5.6:  Measurements  of  mass  differences  between  Ec's  and  the  A+.  All  numbers 
are  in  units  of  MeV/c2. 


Experiment 

M(E++)  -  M(A+) 

M(E+)  -  M(A+) 

M(E»)  -  M(A+) 

ARGUS[41] 
CLEO-II[42,  43] 
E687[45] 
E791[44] 
CLEO-II(This) 

168.2  ±  0.5  ±  1.6 

168.2  ±0.3  ±0.2 
167.6  ±0.6  ±0.6 
167.8  ±0.3  ±0.2 

167.3  ±0.1  ±0.2 

168.5  ±0.4  ±0.2 
166.9  ±0.5  ±0.3 

167.0  ±  0.5  ±  1.6 

167.1  ±0.3  ±0.2 
166.6  ±  0.5  ±  0.6 
167.4  ±0.3  ±0.2 

167.2  ±0.1  ±0.2 

M(E+)  -  M(A+)  by  about  0.13  MeV/c2  and  Af(E*+)  -  M(A+)  by  about  0.38 
MeV/c2. 

Reconstructing  mass  M(A+7r)  using  vertex  constraint  on  the  A+  and  the  in- 
dependent pion  shifts  M(EC)  -  M(A+)  by  about  0.05±0.10  MeV/c2.  For  the  Ec, 
M(A+tt+)  -  M(A+)  and  M(A+tt+)  -  M(A+)  are  about  0.2±0.10  MeV/c2  higher 
than  M(A+7r")  -  M(A+)  and  M(A+7r~)  -  M(A+),  respectively. 

For  Ec's  which  have  clean  signals  and  narrow  widths,  fitting  errors  are  less  impor- 
tant than  the  errors  from  other  sources,  but  for  E*'s  which  have  very  wide  widths, 
fitting  errors  dominate  the  systematic  uncertainties. 

Based  on  Table  5.1,  Table  5.4,  and  more  other  fitting  results  of  M(A+7r)-M(A+), 
we  conclude  that  the  systematic  errors  due  to  fitting  procedures  are  ±0.03  MeV/c2 
for  both  M(E++)-A/(A+)  and  M(E°)-M(A+),  ±0.2  MeV/c2  for  M(E+)-M(A+), 
±8:1  MeV/c2  for  M(E:++)-M(A+),  ±0.5  MeV/c2  for  M(E*+)-M(Ac+),  MeV/c2 
for  M(E;++)  -  Af (A+). 

To  be  conservative,  we  finally  assign  the  total  systematic  errors  to  be  0.2  MeV/c2 
for  M(E++)  -  M(A+)  and  M(EC°)  -  M(A+),  0.3  MeV/c2  for  M(E+)  -  M(A+),  ^ 
MeV/c2  for  M(EC*++)  -  M(A+),  ±0.7  MeV/c2  for  M(E^)  -  M(A+),  ±JJ  MeV/c2 
for  M(E:++)  -  M(A+). 
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Table  5.7:  Measurements  of  mass  splitting  of  the  Ec  triplet  and  comparison  with 
theoretical  calculations.  All  numbers  are  in  units  of  MeV/c2. 


M{Zt  )  -  M(s") 

IT/       -4-  \             H   I"  /       0  \ 

M(£+)  -  M(E°) 

Theory 

Chan[46j 

0.4 

-0.7 

Hwanef47l 

3.0 

-0  5 

Wright  [48] 

-1.4 

-2.0 

Deshpande[49] 

-3.3 

-2.5 

Sinha[50] 

1.5 

-0.3 

Capstick[51] 

1.4 

-0.2 

Experiments 

ARGUS[41] 

1.2  ±0.7  ±0.3 

E791[44] 

0.38  ±0.40  ±0.15 

CLEO-II[42,  43] 

1.1  ±0.4  ±0.1 

1.4  ±0.5  ±0.3 

CLEO-II(this) 

0.1  ±0.1  ±0.1 

-0.3  ±0.5  ±0.1 

Table  5.6  summarizes  the  measurements  of  the  Ec  masses  from  different  experi- 
ments. Note  that  we  have  larger  errors  on  M(E+)  -M(A+)  than  the  previous  CLEO 
measurement.  This  is  because  we  have  fitted  the  mass  difference  with  a  floated  nat- 
ural width  (Breit-Wigner  function)  instead  of  the  assumed  negligible  natural  width 
in  the  previous  works.  For  E++  and  E°  mass  differences,  we  have  better  statistical 
errors  due  to  the  significant  improvement  of  the  resolution  of  soft  charged  pions. 

The  experimental  measurements  of  mass  splitting  among  the  three  Ec  states 
together  with  some  theoretical  predictions  from  these  models  are  listed  in  Table  5.7. 
Our  M(E°)  result  is  in  vary  good  agreement  with  previous  results.  Our  M (E++) 
result  is  lower  than  the  average  of  the  previous  results,  but  nearly  equal  to  M(E°). 
Only  CLEO  has  observed  E+,  and  our  measured  mass  is  obviously  lower  than  the 
previous  result. 

Theoretical  predictions  about  Ec  isospin  mass  splitting  range  from  -6.5  MeV/c2 
to  ±18.0  MeV/c2.  Table  5.7  lists  the  predictions  which  are  close  to  our  measure- 
ments. These  models  have  considered  the  following  effect:  (1)  the  intrinsic  mass 
difference  between  the  u  and  d  quarks;  (2)  electromagnetic  interactions  consisting 
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Table  5.8:  Measurements  of  mass  difference  between  E*'s  and  the  A+,  in  unit  of 
MeV/c2. 


Experiment 

M(E;++)  ~  M(AC+) 

M(Z*+)  -  M(A+) 

M(E*°)  -  M(A+) 

CLEO-II[21] 
CLEO-II(this) 

234.5  ±1.1  ±0.8 
231.4  ±0.8^ 

229.2  ±  1.3  ±0.5 

232.6  ±  1.0  ±0.8 
231.9  ±  1.2j£? 

of  Coulomb  interactions  between  every  pair  of  quarks  and  spin-spin  interactions;  (3) 
contributions  of  QCD  effects  (in  some  of  models). 

CLEO  is  the  only  experiment  which  has  reported  the  measurements  on  E*'s. 
Table  5.8  compared  the  mass  differences  measured  in  this  work  with  the  recent 
reported  values.  Both  measurements  used  the  same  datasets,  but  the  latter  used 
the  recompressed  data  in  which  the  charged  tracks  are  processed  differently  and 
supposed  to  be  measured  better.  In  this  work  we  have  roughly  the  same  errors  since 
we  could  not  take  much  advantage  of  better  charged  track  resolution  due  to  the 
large  natural  width  of  these  states.  We  have  strongly  confirmed  the  existence  of  the 
E*++  and  the  E*°  and  shown  the  first  evidence  of  the  existence  of  the  E*+. 

5.4  Widths 

As  discussed  in  Chapter  3,  the  measurement  of  a  particle's  natural  width  depends 
on  the  Monte  Carlo  simulated  resolution  of  the  detector,  especially  when  the  width 
is  narrow.  Also  as  stated  in  Chapter  3,  we  know  that  the  real  detector  resolution  of 
the  M(D*+)  -  M(D°)  measurement  is,  in  the  worst  case,  15%  larger  than  the  Monte 
Carlo  predicted  one.  We  also  found  that  the  A;ti(2625)  actual  resolution  is  likely 
bigger,  and  can  be  up  to  10%  larger  than  the  Monte  Carlo  predicted  resolution 
if  A^i  (2625)  intrinsic  width  is  negligible.  Since  the  momentum  spectrum  on  the 
soft  charged  pion  for  the  E++  or  the  E°  is  similar  to  the  two  independent  pions  in 
the  A+1(2625)->A+7r+7r-  decay,  we  should  expect  that  the  Monte  Carlo  predicted 
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Table  5.9:  Measurements  of  widths  of  Ec's. 


r(E++)(MeV/c2) 

r(E+)(MeV/c2) 

r(E°)(MeV/c2) 

Fitted  Value 
Limit 

2.1  ±0.3t£ij 
>  0.6  (90%  C.L.) 

A  r-KU  +  l.U 

4-°-1.7-1.5 

0  0+0.4+0.3 
^■^-O.S-O.b 

>  0.7  (90%  C.L.) 

resolution  on  M(A+7r±)  -  M(A+)  may  be  lower  than  the  actual  one  by  roughly  the 
same  amount  and  take  this  into  account  when  we  evaluate  the  systematic  errors. 

The  Monte  Carlo  resolution  on  both  the  E++  and  the  E°  is  1.04  MeV/c2.  To  make 
sure  that  we  have  seen  the  finite  widths  of  these  states,  we  fit  the  M(A+7r+)  -  M(A+) 
signal,  which  is  slightly  narrower  than  M(A+7r~)  -M(A+)  signal,  by  a  width-floated 
single  Gaussian.  The  fitted  widths  range  from  1.81  to  1.95  MeV/c2,  depending  on 
how  the  background  is  parameterized.  We  get  obviously  bad  fits  if  we  fit  the  signal 
by  a  single  Gaussian  with  a  fixed  width  lower  than  1.6  MeV/c2.  We  have  also  tried 
to  fit  the  signal  by  a  double-Gaussian,  with  a  fixed  shape  determined  by  Monte 
Carlo,  and  could  only  get  reasonable  fits  by  enlarging  the  Monte  Carlo  width  by  at 
least  50%.  Therefore  we  believe  that  our  detector  resolution,  which  should  not  be 
poorer  than  1.35  MeV/c2,  has  enabled  us  to  extract  finite  widths  of  Ec's  from  the 
fits. 

Even  conservatively,  we  are  very  certain  that  the  detector  resolution  must  be 
better  than  1.35  MeV/c2,  which  is  30%  worse  than  the  Monte  Carlo  predicted  res- 
olution. Then  we  looked  at  the  x2  distribution  of  fits  to  the  mass  difference  using 
different  Breit-Wigner  widths  convoluted  with  a  single  Gaussian  with  fixed  1.35 
MeV/c2  width.  It  can  be  concluded  that  at  90%  confidence  level,  T(E++)  >  0.6 
MeV/c2  and  r(E°)  >  0.7  MeV/c2. 

Table  5.9  gives  the  results  on  Ec  widths  measurements.  The  fitted  values  and 
their  errors  are  based  on  Table  5.1  and  Table  5.4.  Uncertainties  of  the  detector 
resolution  dominate  the  systematic  errors.  No  other  experiments  has  reported  Ec 
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Table  5.10:  Measurements  of  widths  of  E*'s. 


Experiment 

t(e;++) 

t(e;+) 

t(e:°) 

CLEO-II[21] 
CLEO-II(this) 

17.9^  ±  4.0 
14  ±  3  ±  3 

°-5-2 

13.0ljj;S  ±  4-0 
20±l  ±  5 

Table  5.11:  CLEO  Measurements  of  fragmentation  functions  using  the  Peterson 
function  for  the  Ec  and  the  E*. 


Particles 

E++ 
E° 

E++  +  E° 
EcAverage 

o.2ai*|B 

0  9<S+0-18 
u-zo-0.09 

n  9Q+0-05 
u-zy-0.04 

0.28  ±0.03 
0.28  ±0.03 

e:° 

K++  +  K° 

E*++  E*°  Average 

0.36^;{i 
0.32±g?? 
0.34l°;ii 
o  n+011 

widths  measurements,  and  all  previously  published  Ec  observations  have  the  detector 
resolution  no  better  than  2  MeV/c2. 

Our  results  of  the  E*  widths,  together  with  the  CLEO  previously  measured 
values,  are  tabulated  in  Table  5.10.  The  new  results  are  consistent  with  the  previous 
ones.  As  we  explained  in  the  last  section,  this  measurement  cannot  do  much  better 
than  the  previous  work  with  E*.  The  E*+  width  we  first  measured  is  not  very 
meaningful,  since  the  signal  needs  to  be  further  confirmed  and  the  error  is  too  large. 

5.5    Fragmentation  Functions 

Figure  5.3  and  Figure  5.4  show  the  distributions  for  the  Ec  and  the  £*.  The 
yields  in  the  figures  are  the  numbers  in  Table  5.2  and  Table  5.5  with  the  correct 
systematic  errors  added,  and  they  are  normalized  to  the  total  yield  in  the  0.4  -  1.0 
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Figure  5.3:  Normalized  production  spectra  of  scaled  momentum  for  Ec  samples.  The 
data  for  the  £++  (triangles),  the  E+ (circles),  and  the  £°  (inverted  triangles)  are  fitted 
to  the  Peterson  functions  overlaid  as  solid,  dotted,  and  dashed  curves,  respectively. 
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Figure  5.4:  Normalized  production  spectra  of  scaled  momentum  for  £*  samples. 
The  data  for  the  £*++ (triangles),  and  the  £*°  (inverted  triangles)  are  fitted  to  the 
Peterson  functions  overlaid  as  solid  and  dashed  curves,  respectively. 
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xp  range.  The  distributions  are  fitted  using  the  Peterson  function.  Table  5.11  gives 
the  measured  values  of  e  for  Ec's  and  E*'s. 

Our  value  of  ep(E++  +  E°),  0.28±0.03,  is  in  agreement  with  the  previous  result, 
0.29  ±  0.06,  done  by  ARGUS[41]  also  in  the  10  GeV/c2  e+e~  experiment.  Our  E* 
values  are  also  agree  with  the  previous  value,  ep(E*++  +  E*°)  =  0.301°;^ . 

Based  on  the  measurements,  it  seems  that  E^'s  are  fragmented  slightly  harder 
than  A+,  which  has  eP  =  0.25±0.03[1]  from  10GeV/c2  e+e~  annihilation.  A^'s  are 
obviously  fragmented  much  harder.  This  is  possibly  because  many  A^'s  are  directly 
created  from  the  annihilation,  whereas  a  much  larger  fraction  of  E^'s  are  the  decay 
products  from  higher  excited  states. 

5.6    Production  Ratios 

To  calculate  the  ratio  of  the  A+  production  from  Ec  and  E*  decays  relative  to 
the  total  inclusive  production,  the  reconstruction  efficiencies  of  the  A+,  the  Ec,  and 
the  E*  from  the  Monte  Carlo  simulation  are  used.  Our  calculation  covers  the  xp 
region  of  0.4  -  1.0  for  SW++  and  E^0,  and  0.5  -  1.0  for  the  E+.  Since  we  measure 
the  productions  by  ratio,  most  of  the  systematic  errors  on  the  A+  reconstruction 
cancel.  We  estimate  that  the  efficiency  errors  are  3%  for  E(*)++  or  E(*)0  and  6%  for 
the  E+.  The  statistical  errors  on  the  Monte  Carlo  efficiency  calculation  are  tiny. 

The  uncertainty  of  the  detector  mass  resolution  brings  large  errors  in  the  calcu- 
lation of  E++  and  E°  productions.  We  estimate  this  uncertainty  by  reducing  and 
increasing  the  Monte  Carlo  resolution  by  25%  and  comparing  the  areas  of  the  Breit- 
Wigner  distributions.  We  estimate  this  error  to  be  t\Vo.  Other  fitting  errors  are 
mostly  due  to  different  background  shapes,  and  they  are  tested  to  be  no  more  than 
3%. 
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Table  5.12:  Measurements  of  the  Ec  production  relative  to  total  inclusive  A+  pro- 
duction from  different  experiments.  All  numbers  are  in  units  of  %. 


Experiment 

Beam 

E691 

145  GeV  7 

5±3±2 

13±4±2 

E687 

220  GeV  7 

6.7  ±  1.9 

7.8  ±2.1 

NA32 

230  GeV  7T" 

<  5.2(90%C.L.) 

<  5.3(90%C.L.) 

E791 

500  GeV  tt- 

6.7±  1.0  ±  1.0 

7.9  ±  1.1  ±  1.0 

CLEO-I.5 

10  GeV  e+e~ 

6±3 

6±3 

CLEO-II(this) 

10  GeV  e+e- 

7.2  ±0.5  ±0.7 

7.3  ±0.5  ±0.7 

Fitting  procedures  dominate  the  systematic  errors  for  the  E+  and  E*'s,  and  the 
uncertainties  on  the  resolution  are  much  less  important.  We  estimate  the  fitting 
errors  are  t\2%  for  the  E+,  ±25%  for  the  two  E*'s. 

We  extrapolate  the  production  distribution  down  to  xp  =  0  of  the  Peterson 
function.  The  uncertainties  on  fragmentation  function  measurements  bring  large 
errors.  This  is  not  only  because  we  lack  the  statistics  to  measure  the  fragmentation 
function  accurately,  but  also  because  there  is  no  reason  to  expect  the  Peterson 
function  to  perfectly  represent  the  data.  We  estimate  the  errors  in  the  extrapolations 
are  6%  for  the  E++  and  the  E°,  22%  for  the  E+  and  18%  for  E*'s. 

We  finally  assign  the  total  systematic  errors  on  production  ratio  measurements 
to  be  tl°%  for  the  E++and  the  Sc°,  ±™%  for  the  E+  and  ±31%  for  E*'s. 

The  background  is  not  simply  formed  like  the  phase  space  of  arbitrary  particle 
combinations.  We  already  know  that  there  is  wide  enhancement  around  the  region 
between  140  MeV/c2  to  210  MeV/2  due  to  a  feed-down  of  A+!(2625)->-A+7r7r  decay. 
Very  likely  there  are  some  other  unknown  physical  processes  making  non-phase  space 
contributions  elsewhere  in  the  mass-difference  spectra,  for  instance,  feed-down  from 
more  higher  states  of  the  A+  and  the  Ec  which  decay  into  A+'s. 

Table  5.12  lists  the  measurements  of  <tEc/<ta+  for  the  E++  and  E°.  The  E+ 
relative  production  ratio  has  not  been  measured  by  any  previous  experiments,  and 
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Table  5.13:  Measurements  of  E*  production  ratio  relative  to  total  inclusive  A+ 
production  in  units  of  %. 


Experiment 

(7e;++/(7a+ 

crS^°/f7A+ 

CLEO-II  1996 
CLEO-II  (this) 

5.9lJi  ±  1.8 

7.8+l;i  ±  2.4 

12.8^3  ±3.2 
13.71?:?  ±  3.0 

for  the  first  time  we  have  calculated  that 

(TEc+  /cta+  =  10.3  ±  1.6  ±  2.7%. 

Note  if  we  fit  the  signal  with  fixed  natural  width  the  same  as  the  E°  signal,  the 
result  would  be  7.5  ±  1.1  ±  2.0.  Also,  if  we  assume  that  Ec's  do  not  have  natural 
widths  and  fit  the  signals  with  width-floated  single  Gaussians,  the  results  on  the 
El+  and  the  E°  are  15%  smaller  and  the  result  on  the  E+  is  20%  smaller. 

Table  5.13  lists  the  cr^-/aA+  measurements  from  CLEO.  The  new  result,  which 
agrees  with  the  previous  one,  shows  that  the  E*  and  the  Ec  are  nearly  equally 
produced. 

According  to  their  spins,  which  are  1/2  for  the  Ec  and  3/2  for  the  E*,  the  number 
of  E*'s  created  should  be  two  times  as  large  as  the  number  of  Ec's  from  the  e+e~~ 
annihilation.  Then  why  is  the  number  of  Ec's  in  CLEO  experiments  even  larger 
than  1/2  of  E*'s?  As  the  E*  has  higher  mass,  its  creation  should  be  slightly  more 
suppressed  than  the  creation  of  a  Ec;  however,  we  expect  this  to  be  a  small  effect. 
We  think  the  real  reason  is  very  likely  that  there  are  more  higher  states  which  decay 
to  the  Ec  than  decay  to  the  E*.  For  instance,  we  know  the  decay  Al1(2593)->Ec7r 
contributes  about  20%  of  the  Ec  production. 


116 


5.7    New  Measurements  of  A+1(2593)^A+7r+7r 


According  to  our  measurements,  A+1(2593)  decays  almost  100%  to  Ec7r,  and  using 
our  average  £++  and  £°  mass  value,  we  can  calculate  that  the  kinematic  threshold 
of  +7r~  and  S°7r+  is  at  a  mass  difference  of  306.86  MeV/c2.  Since  this  threshold 
is  very  close  to  the  central  mass  of  the  resonance,  the  shape  of  the  signal  will  not 
be  an  "ordinary"  Breit- Wigner.  We  believe  this  is  the  reason  why  a  Breit- Wigner 
convoluted  with  a  Gaussian  resolution  does  not  fit  the  A^(2593)— >A+7r+7r"  signal 
well  as  shown  in  Figure  5.5.  The  fit  gives  a  ln(Likelihood)*  value  of  48.0  for  30 
d.o.f. 

If  we  assume  that  A^(2593)  decay  100%  to  Ec7r,  a  reasonable  way  to  fit  the 
A+1(2593)— >A+7r+7r~  signal  is  to  use  a  relativistic  Breit- Wigner  form[l,  52]  with  the 
amplitude 


where  M  is  the  invariant  mass  of  the  Ec7r  in  the  final  state  that  results  from  the 
decay  of  the  A+!(2593),  Mr  is  the  mass  of  the  A+j(2593),  and  rr(7r°)  and  rr(?r±)  are 
the  mass-dependent  widths  of  A+x(2593)  corresponding  to  the  decay  A+7r°7r°  and 
A+7r+7r",  respectively.  Note  both  rr(7r°)  and  T^n*)  are  in  the  denominator  but 
only  rr(7r±)  is  in  the  numerator.  This  is  because  we  must  put  the  "full"  width  in 
the  denominator,  and  A^(2593)  has  both  A+7r°7r°  and  A+7r+7r_  decay  channels,  but 
here  we  are  only  measuring  the  signal  of  the  A+7r+7r~  decay.  We  need  a  factor  \JM  in 
the  numerator  of  the  Breit- Wigner,  because  the  phase  space  factor  is  a  factor  times 
dM2  =  MdM,  and  we  are  fitting  the  data  as  a  function  of  M.  The  factor  would 
not  be  present  if  we  fit  the  data  as  a  function  of  M2.  For  A;|i(2593),  M  »  Tr,  so 
this  y/M  factor  virtually  does  not  affect  the  Breit- Wigner  shape.  The  normalization 
factor  in  equation  5.1  is  omitted,  and  it  is  evaluated  numerically  in  the  fit  procedure. 
*ln(Likelihood)  is  the  same  as  x2  for  large  statistics. 


(5.1) 


M2  -  M2  -  iMr{rr(ir°)  +  rr(7r±)]' 
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Figure  5.5:  The  fit  to  the  A+1(2593)->A+7r+7r-  signal  using  a  simple  Breit-Wigner 
function  convoluted  with  the  detector  resolution.  The  mass-difference  spectrum  and 
the  fit  are  the  same  as  Figure  4.1,  but  here  we  only  fitted  A+x(2593)  in  order  check 
the  likelihood  of  the  fit,  and  1-MeV/c2  bin  size  is  shown  instead  of  1-MeV/c2  in 
order  to  get  a  clearer  view  near  the  threshold. 
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For  this  relativistic  Breit-Wigner  function,  the  width  parameter  of  the  resonance 
is  mass  dependent.  For  a  decay  by  a  partial  wave  of  orbital  angular  momentum  L, 
the  width  can  be  described  as 

Tr  =  (p*/M)2Lp*a,  (5.2) 

where  p*  is  the  momentum  in  the  rest  frame  of  the  resonance  of  either  one  of  the 
daughters,  and  a  is  a  mass-independent  parameter.  It  has  been  tested  in  experiments 
that  the  width  parameterization  of  equation  5.2  can  reasonably  fit  to  the  mass- 
dependent  width  of  P- wave  decay  of  the  p[53,  54]. 
A^(2593)— )-E°7r~  is  two-body  decay,  so  we  have 

p*^-)  =  [(M2  -  (M£o  +  A/„-)2)(M2  -  (Mso  -  A/,-)2)]2  (5  3) 

p*(£+7T°)  and  p*(E++7r+)  can  be  expressed  similarly.  In  our  fit  procedure  we  as- 
sume that  Ms++  =  Mso;  thus  we  have  p*^)  =  p*(£°7r_)  =  p*(E++7r+).  Since 
Acl(2593)— >Ec7r  is  5-wave  decay,  according  to  equation  5.2, 

rr(*±)  =  P^Mtt*),  (5.4) 
rr(7r°)   =  pV)a(A  (5.5) 

From  isospin  symmetry,  we  can  assume  that  a(7r°)  =  fax*).  Now  we  have  three 
parameters  to  fit  the  data:  Mr,  a(7r±),  and  the  area  of  the  signal. 

To  fit  the  A+1(2593)-)'A+7r+7r_  signal,  we  used  the  relativistic  Breit-Wigner  func- 
tion as  in  the  equation  5.1,  smeared  with  a  single  Gaussian  and  a  regular  Breit- 
Wigner  function.  The  Gaussian  is  used  to  parameterize  the  detector  resolution  of 
A+1(2593)^A+7r+7r-,  1.28  MeV/c2,  determined  by  Monte  Carlo  studies;  and  the 
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Figure  5.6:  New  fit  to  the  A^i(2593)-)-A+7r+7r_  signal  using  the  relativistic  Breit- 
Wigner  function  described  in  the  text  convoluted  with  the  detector  resolution  and 
the  Ec  width. 

Breit-Wigner  represents  the  natural  width  of  E++  or  E°,  with  an  average  value  of 
2.2  MeV/c2.  To  input  a  Ec  mass,  we  used  the  mass  difference,  167.3  MeV/c2,  which 
is  our  measured  average  of  E++  and  E°.  The  background  is  fitted  by  a  second  or- 
der polynomial  with  the  area  and  shape  floated.  Figure  5.6  shows  such  a  fit.  The 
ln(Likelihood)  value  for  the  fit  is  30.4  for  30  d.o.f,  much  better  compared  with 
48.0  from  the  fit  using  the  regular  Breit-Wigner  function. 

The  results  of  the  fit  shown  in  Figure  5.6,  which  are  used  as  our  presented  values, 
are  listed  in  the  first  row  of  the  values  in  Table  5.14.  In  the  following  rows  of  the 
table,  we  tabulated  the  results  obtained  by  different  types  of  signal  parameterization. 
The  first  three  numbers  in  each  row  are  input  values  of  the  fit  and  the  last  three 
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Table  5.14:  Fit  results  of  A+1(2593)— ►A+7r+7r_  signal  in  Figure  5.6.  The  Ec  mass 
difference  (with  A+)  AMSc,  the  Ec  width  rEc,  and  the  A^1(2593)-)'A+7r+7r_  mass 
resolution  erA+  are  fixed  as  the  input  to  the  fitting  function.  The  mass  difference 

(with  A+)  of  the  A^(2593),  a(7r±)  and  the  yield  of  the  signal  are  the  measured 
parameters.  All  mass  differences  and  widths  in  the  table  are  shown  in  units  of 
MeV/c2. 


Fixed  parameter 

Fitted  parameter 

A.Uv, 

rv, 

AMA+ 

a(7r±) 

Yield 

167.3 

2.2 

1.28 

305.3  ±0.4 

83™ 

195i£ 

167.0 

2.2 

1.28 

305.2 

96 

194 

167.6 

2.2 

1.28 

305.4 

73 

198 

167.3 

0.7 

1.28 

305.9 

16 

164 

167.3 

1.5 

1.28 

305.0 

68 

187 

167.3 

2.8 

1.28 

305.4 

98 

195 

167.3 

2.2 

1.00 

304.8 

90 

195 

167.3 

2.2 

1.70 

305.7 

82 

183 

numbers  are  the  extracted  parameters.  In  each  row  except  the  first,  there  is  a  number 
shown  in  boldface  in  order  to  emphasis  that  this  is  the  number  we  reasonably  altered 
from  the  central  value.  Note  that  the  mass,  but  not  the  mass  difference,  is  used  in 
the  fit.  So  A+  mass  is  actually  used  here.  We  proved  that  varying  mass  of  A+  by 
a  few  MeV/c2  virtually  does  not  changed  the  measured  parameters.  Although  the 
A+j  (2593)  mass  is  measured,  we  still  present  the  mass  difference  since  it  is  a  more 
accurate  value. 

From  Table  5.14,  we  see  that  the  measured  mass  slightly  varies  when  other 
input  parameters  change.  So  all  input  parameters,  within  reasonable  ranges,  do 
not  influence  the  value  of  the  mass  very  much.  Also  note  that  the  yield  measured 
here  is  almost  the  same  as  the  measurement  described  in  Chapter  4  with  similar 
statistical  error  and  fitting  error,  thus  all  the  A+!(2593)  related  branching  ratios  and 
production  ratios  measured  in  Chapter  4  still  hold. 
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a(7r±)  is  very  sensitive  to  changes  of  any  input  parameter,  but  we  also  notice  that 
the  shape  of  the  resonance  varies  only  slightly  when  a(7r±)  changes.  It  seems  that 
the  a  parameter  is  not  proportional  to  the  "width"  of  the  signal,  and  the  value  of  a 
is  very  sensitive  to  any  change  of  parameterization  near  the  threshold.  We  think  the 
measurement  of  the  a  parameter  in  our  fit  is  not  very  reliable,  and  the  systematic 
error  we  quote  is  large. 

Using  the  measured  value  of  a(7r±)  and  assuming  a(7r±)  =  2a(7r°),  we  can  predict 
the  branching  ratio  r(A+7r°7r°)/r(A+7r+7r-)  using  the  Asw^^)  formula  in  equation 
5.1  and  the  similar  formula  for  £+7r°  decay, 

0  yjMrr(ir°) 
ABw{n  ]  ~  M?-M*-iMr[rr{*°)  +  rr{«±)Y  (5'6) 

We  found  that 


r(Ac+7rV)/r(Ac+7r+7r-)  = 


2.2   if  a(7r±)  =  100, 
3.1   ifa(7r±)  =  84, 
4.0   ifa(7r±)  =  60. 


This  proves  that  the  range  of  a(ir±)  in  Table  5.14  is  consistent  with  our  our  measured 

r(A+7r°7r°)/r(A+7r+7r-)  in  Chapter  4,  which  is  2.6  ±  0.6  ±  0.7. 

Similar  to  the  A+1(2593)->A+7r+7r_  case,  we  can  also  fit  the  A,?j(2593)->A+7r07r0 

signal  using  a  Breit-Wigner  function  with  the  amplitude  Abw(tt°)  in  equation  5.6. 

Note  that  Abw(-k°)  is  very  close  to  a  regular  Breit-Wigner  amplitude  since  the  £+7r° 

threshold  is  at  about  302  MeV/c2,  4  MeV/c2  below  the  center  mass,  and  therefore 

fitting  the  A+1(2593)->A+7r°7r°  signal  using  Abw(-kq)*  obtains  almost  the  same  mass 

"Another  way  to  do  the  measurement  is  to  fit  the  A+7r+7r~  and  A+7r°7r°  signals  simultaneously 
by  constraining  them  to  have  the  same  a's,  but  since  A+7r°7r°  signal  is  much  poor,  the  measured 
parameters  would  be  almost  determined  by  the  A+7r+7r_  signal.  We  think  it  is  more  reasonable 
to  fit  them  separately. 
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Figure  5.7:  New  fit  to  the  A+1(2593)^A+7r°7r°  signal  using  the  relativistic  Breit- 
Wigner  function  described  in  the  text  convoluted  with  the  detector  resolution  and 
the  Ec  width. 


and  yield  as  the  fit  in  Figure  4.3.  Figure  5.7  shows  such  a  fit.  In  the  fit,  we  used 
a  convoluting  Gaussian  with  a  width  of  2.4  MeV/c2,  representing  the  A+7r07r°  mass 
resolution,  and  the  same  convoluting  Breit-Wigner  as  we  used  to  fit  the  A+7r+7r~ 
signal*.  However,  due  to  low  statistics,  this  fit  measures  the  a  parameter  with 
enormous  statistical  error,  so  we  cannot  get  more  information  about  ABw(n°)  and 
Abw^)  by  directly  fitting  them  with  the  A+7r°7r°  signal.  Figure  5.7  shows  a  very 
good  fit  by  fixing  Abw^)  to  be  83. 


'Since  the  S+  width  is  very  poorly  measured,  we  assume  that  it  is  the  same  as  the  E++  and 
E°.  The  fit  does  not  change  much  even  if  we  use  =  4.5  MeV/c2. 


123 

Finally,  using  a  new  method  never  used  on  A+1(2593)  in  any  previous  experi- 
ments, we  have  measured  A^(2593)  with  the  results 

M(A+(2593)  -  M(A+)  =  305.3  ±  0.4  ±  0.6  MeV/c2, 

a(7r±)  =  831J0  ±  40. 

The  fitting  error,  0.5  MeV/c2evaluated  from  Table  5.14,  dominates  the  total  sys- 
tematic error  of  A^(2593)  measurements.  The  error  estimated  from  the  momentum 
scale  study  is  0.2  MeV/c2.  The  value  of  mass  is  significantly  lower  than  the  pre- 
viously measured  values  using  a  normal  Breit- Wigner  as  the  fit  function  listed  in 
Table  4.7.  But  the  value  is  consistent  with  our  mean  value  of  A^(2593)— >A+7r07r0, 
306.3  ±0.7  MeV/c2. 

The  fit  using  the  lower  limit  of  the  i\c,  0.7  MeV/c2,  gives  a  very  small  a. 
Although  the  fit  is  acceptable,  it  is  not  used  to  evaluated  the  systematic  error  of 
the  a(7r±)  value  because  the  0.7  MeV/c2  limit  is  far  beyond  one  unit  of  standard 
deviation.  More  careful  studies  are  needed  about  the  a  parameter  of  A^(2593)  in 
both  theoretical  models  and  experiments. 


CHAPTER  6 
SUMMARY  AND  DISCUSSION 


In  this  research,  we  have  made  a  complete  study  of  eight  charmed  baryons:  the 
A^j  doublet,  which  decays  to  A+7T7T,  and  the  Ec  triplet  and  the  E*  triplet  which 
decay  to  A+7r. 

With  a  large  data  sample  and  improved  detector  resolution,  we  measured  the  Ec 
masses  more  accurately  than  previous  experiments,  and  the  results  are  167.3  ±  0.1  ± 
0.2MeV/c2  for  E++,  166.9  ±  0.5  ±  0.3MeV/c2  for  E+,  and  167.2  ±  0.1  ±  0.2MeV/c2. 
Our  results  do  not  indicate  any  noticeable  isospin  mass  splittings  among  this  triplet. 

Studying  the  mass  splitting  among  baryons  containing  different-valence  quarks 
is  a  good  way  to  understand  hadronic  structure  and  the  forces  which  determine  the 
structure.  The  mass  splittings  in  the  charm  baryon  system  are  especially  interesting 
since  the  large  mass  of  the  c  quark  simplifies  the  component  structure  in  theoretical 
models  so  the  experimental  tests  are  more  meaningful.  The  Ec  is  the  first  baryon 
isospin  multiplet  for  which  the  intrinsic  light-quark  mass  difference  is  not  the  domi- 
nating factor  of  the  isospin  mass  splitting.  However,  the  theoretical  predictions  vary 
largely  for  different  models.  Our  results  support  that  the  splittings  are  very  small, 
and  the  numbers  are  close  to  the  prediction  by  Chan [46]. 

In  this  work,  the  Ec  widths  were  measured  for  the  first  time,  and  the  results  are 
2.1  ±  0.3i°j  MeV/c2  for  £++,  4.5^};°  MeV/c2  for  £+  and  2.3±g:|±g;|  MeV/c2 
for  E°.  These  values  are  in  very  good  agreement  with  the  prediction  of  heavy-quark 
chiral  perturbation  theory[16].  We  measured  the  E*  widths  r(E*++)  =  14  ±  3  ± 
3  MeV/c2  and  r(E*°)  =  20±f  ±  5  MeV/c2.  These  are  consistent  with  the  previous 
measurement  by  CLEO[21].  Note  that  based  on  equations  1.3  and  1.5,  the  Ec  and 
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E*  widths  are  correlated  with  each  other  through  the  coupling  constant  g2-  Our 
results  agree  with  the  g2  value  predicted  by  many  theoretical  models. 

We  have  confirmed  the  previous  E*++  and  E*°  observation  made  by  CLEO, 
with  masses  measured  to  be  231.4  ±  0.8^|  MeV/c2  for  the  £*++  and  231.9  ± 
1.2tJ;|  MeV/c2  for  the  E*°,  and  the  widths  measured  to  be  14  ±  3  ±  3  MeV/c2  for 
the  £*++  and  20^  ±  5  MeV/c2  for  the  E*°,  all  consistent  with  previous  results.  In 
addition,  for  the  first  time  we  observed  the  E*+  evidence,  with  the  mass  measured 
to  be  229.2  ±  1.3  ±  0.5  MeV/c2  and  the  width  measured  to  be  5t^  MeV/c2. 

We  have  measured  the  ratios  of  the  A+  productions  from  relative  to  the 
total  inclusive  A+  production  from  the  10GeV/c2  e+e~  production,  and  the  results 
are  7.2  ±  0.5  ±  0.7%  from  E++,  10.3  ±  1.6  ±  2.7%  from  E+,  7.3  ±  0.5  ±  0.7%  from 
E°,  5.91^  ±  1.8%  from  E*++,  and  7.8±i;§  ±  1.8%  from  E^0.  We  are  unable  to 
make  a  E*+  measurement  since  the  signal  is  too  poor.  This  is  the  first  time  that 
aE+/aA+  's  measured,  and  all  the  other  ratios  measured  are  in  agreement  with 
previous  experiments.  Based  on  our  measurements,  we  can  conclude  that  25  ±  4% 
of  A+'s  are  from  Ec  decays.  If  we  assume  the  three  E*'s  are  equally  produced,  then 
about  20%  of  A+'s  are  produced  from  E*  decays. 

The  following  evidences  provide  very  strong  proof  that  these  two  states  are  indeed 
A+!  baryons  with  light  degrees  of  freedom  having  L  =  1:  (1)  the  A+7r°  signals  are  not 
seen;  (2)  the  Ec7r  channel  transition  dominates  the  A+(2593)  decay  and  the  same 
channel  is  measured  small  for  A+j(2625);  and  (3)  the  A+7r°7r°  channel  is  observed 
for  both  A^i's. 

This  is  the  first  time  the  A+7r°7r°  decay  channel  has  been  observed,  and  the 
relative  decay  ratio  B(A+7r°7r°)/,B(A+7r+7r-)  is  measured  to  be  2.6  ±  0.6  ±  07  for 
A+(2593)  and  1.0  ±  0.3l<^  for  A+(2625).  The  value  for  A+(2593)  does  not  agree 
with  isospin  symmetry  but  can  be  explained  by  phase-space  suppression. 
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The  A+i(2625)  mass  was  measured  to  be  341.80  ±  0.10  ±  0.35  MeV/c2.  We 
measured  the  upper  limit  of  the  A+x(2625)  width  to  be  1.4  MeV/c2  at  90%  confidence 
level.  These  values  agree  with  previous  experiments. 

By  fitting  the  A^(2593)— >A+7r+7r~  signal  the  same  way  as  the  previous  ex- 
periments, it  is  measured  to  have  a  mass  308.8  ±  0.4  ±  0.4  MeV/c2  and  a  width 
4.7l};itJ;9  MeV/c2  which  are  consistent  with  previous  results  but  not  consistent 
with  our  A+1(2593)— »A+7r°7r°  measurements,  in  which  we  measured  the  mass  to  be 
306.3  ±  0.7  MeV/c2.  But  these  measurements  are  based  on  the  signal  parameter- 
ization using  an  ordinary  Breit-Wigner  function.  Because  A^j(2593)  resonance  is 
just  at  the  £c7r±  kinematic  threshold,  a  relativistic  Breit-Wigner  function  should  be 
used.  In  this  function,  the  width  of  a  resonance  is  mass  dependent  and  can  only  be 
parameterized  by  the  a  parameter.  Using  the  new  method,  we  obtained  very  good 
fit  and  measured  the  A^(2593)  to  have  a  mass  305.3  ±0.4  ±0.6  MeV/c2.  Assuming 
that  g(7t±)  =  2a(n°)  based  on  isospin  symmetry,  we  measured  a(7r±)  =  83lg°  ±  40. 
These  new  measured  values  are  consistent  with  our  r(A+7r°7r°)/r(A+7r+7r_)  result, 
2.6  ±0.6  ±0.7. 

The  non-relativistic  quark  potential  model [6]  fails  in  its  prediction  of  the  A£ 
masses.  The  relativized  model  is  close,  but  their  10  MeV/c2  spin-orbital  mass  split- 
ting is  too  small. 

The  mass  splitting  between  the  two  A^  states  is  36.5  ±  0.4  ±  0.6  MeV/c2  from 
our  measurements,  about  4  MeV/c2  larger  than  32.4  ±  1.0  ±  0.7  MeV/c2  measured 
by  ARGUS[32].  From  heavy-quark  symmetry,  we  should  expect  that 

M(A+(2625))  -  M(A+(2593))  =  M(JD2*(2460)°)  -  M(A(2420)°). 

The  PDG[1]  average  value  of  M(Z^(2460)°)  -  M(A(2420)°)  is  37 ±3  MeV/c2.  We 
see  that  HQET  works  quantitatively  very  well  in  this  case. 
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To  date,  few  theoretical  models  directly  predicted  the  A^j  widths.  The  most  re- 
cent calculation  was  made  by  Tawfiq  et  al.[55].  They  used  a  light-front  constituent 
quark  model  which  utilized  the  SU(2iV/)xO(3)  diquark  symmetry,  where  Nf  rep- 
resents the  number  of  light  flavors.  Assuming  the  mass  of  a  light  quark  to  be  340 
MeV/c2,  they  predicted  that  r(A+!(2593)  ->  E  c7r)  is  around  1.6  ~  2.5  MeV/c2. 
They  also  predicted  the  D-wave  transition  rate  r(A+1(2625)  — >  Sc7r)  to  be  about 
0.7  MeV/c2.  But  our  measurements  contradict  these  values.  No  model  has  directly 
predicted  the  widths  of  non-resonance  two-pion  transitions.  Better  Monte  Carlo 
simulation  and  detector  resolution  are  needed  to  perform  further  measurement  of 
the  A+!(2625)  width. 

Theoretical  computations[56]  give  small  values  of  of  decay  rates  of  the  radiative 
decays:  0.016  MeV/c2for  A+^2593)  and  0.021  MeV/c2for  A+^2625).  Ac- 
cording to  these  values  we  should  not  be  able  to  see  A^(2593)->A+7  smce  its  width 
is  obviously  large,  and  we  might  be  able  to  see  A;Jj(2593)->A+7  if  T(A+1(2625)) 
A+7T+7T-  is  very  narrow.  We  did  not  find  any  evidence  of  A^-^A+7  decays.  The 
upper  limit  we  measured  is  too  large  to  provide  any  indications. 
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